Karlsruher Institut für Technologie 


Schriftenreihe 
Kontinuumsmechanik im Maschinenbau 


Julian Karl Bauer 


Fiber Orientation Tensors 
and Mean Field Homogenization 


Application to Sheet Molding Compound 


NIT zusishing 


Julian Karl Bauer 


Fiber Orientation Tensors and 
Mean Field Homogenization 


Application to Sheet Molding Compound 


Schriftenreihe 
Kontinuumsmechanik im Maschinenbau 
Band 23 


Karlsruher Institut fur Technologie (KIT) 
Institut fur Technische Mechanik 
Bereich Kontinuumsmechanik 


Hrsg. Prof. Dr.-Ing. habil. Thomas Bóhlke 


Eine Übersicht aller bisher in dieser Schriftenreihe erschienenen Bànde 
finden Sie am Ende des Buchs. 


Fiber Orientation Tensors and 
Mean Field Homogenization 


Application to Sheet Molding Compound 


by 
Julian Karl Bauer 


SICHT Fikisha 


Karlsruher Institut fur Technologie 
Institut fur Technische Mechanik 
Bereich Kontinuumsmechanik 


Fiber Orientation Tensors and Mean Field Homogenization: 
Application to Sheet Molding Compound 


Zur Erlangung des akademischen Grades eines Doktor-Ingenieurs 
von der KIT-Fakultät für Bauingenieur-, Geo- und Umweltwissen- 
schaften des Karlsruher Instituts fur Technologie (KIT) genehmigte 
Dissertation 


von Julian Karl Bauer, M.Sc. 


Tag der mündlichen Prüfung: 29. Juli 2022 
Hauptreferent: Prof. Dr.-Ing. habil. Thomas Seelig 
Korreferent: Prof. Dr.-Ing. habil. Thomas Bóhlke 
Korreferent: Prof. Dr.-Ing. Andrew Hrymak 


Impressum 
NC Scientific 
Publishing 
Karlsruher Institut für Technologie (KIT) 
KIT Scientific Publishing 


StraBe am Forum 2 
D-76131 Karlsruhe 


KIT Scientific Publishing is a registered trademark 
of Karlsruhe Institute of Technology. 
Reprint using the book cover is not allowed. 


www.ksp.kit.edu 


© OO This document - excluding parts marked otherwise, the cover, pictures and graphs — 
baa is licensed under a Creative Commons Attribution-Share Alike 4.0 International License 
(CC BY-SA 4.0): https://creativecommons.org/licenses/by-sa/4.0/deed.en 


© OXO) The cover page is licensed under a Creative Commons 
Rare Attribution-No Derivatives 4.0 International License (CC BY-ND 4.0): 
https://creativecommons.org/licenses/by-nd/4.0/deed.en 


Print on Demand 2023 - Gedruckt auf FSC-zertifiziertem Papier 


ISSN 2192-693X 
ISBN 978-3-7315-1262-2 
DOI 10.5445/KSP/1000152989 


Zusammenfassung 


Die effektiven mechanischen Eigenschaften von faserverstárkten Ver- 
bundwerkstoffen hángen stark von der Mikrostruktur, d. h. von der 
Ausrichtung der Fasern im Verbundwerkstoff, ab. In dieser Arbeit wird 
der Einfluss von Faserorientierungsverteilungen auf die effektiven Ei- 
genschaften von langfaserverstärkten Verbundwerkstoffen untersucht. 


Die vorliegende Arbeit wird von der Deutschen Forschungsgemein- 
schaft (DFG) im Rahmen des interdisziplinären Graduiertenkollegs 
GRK 2078 finanziert. Daher liegt der Fokus auf sogenannten „Sheet 
Molding Compound“ (SMC) Materialien. SMC ist ein Verbundwerk- 
stoff aus langen Glasfasern und duroplastischer Polymermatrix, für 
Leichtbauanwendungen mit flachen Bauteilformen, beispielsweise in 
der Automobilindustrie. 


Gemittelte Information einer Faserorientierungsverteilung kann durch 
Faserorientierungstensoren quantifiziert werden. Ein vollständiges Bild 
des Einflusses dieses Richtungsmaßes auf die effektiven mechanischen 
Eigenschaften wird präsentiert, basierend auf einer neu abgeleiteten 
Vielfalt von Faserorientierungstensoren vierter Stufe. Diese Vielfalt wird 
identifiziert durch lineare invariante Zerlegung von Faserorientierungs- 
tensoren vierter Stufe, was zu einfachen und flexiblen Parametrisie- 
rungen führt. Beschränkungen durch materielle Symmetrie werden in 
diese Parametrisierungen einbezogen und zulässige Parameterberei- 
che werden bestimmt, indem positive Semi-Definitheit eingefordert 
wird. In dieser Arbeit werden Schließungsapproximationen durch die 
Vielfalt der Faserorientierungstensoren vierter Stufe kontrastiert. Mi- 
krostrukturen von SMC werden durch planare Orientierungszustände 
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approximiert. Die Beziehung zwischen planaren Faserverteilungen und 
planaren Faserorientierungstensoren wird im Detail diskutiert. Dazu 
werden tensorielle Fourierreihen in 2D- und 3D-Darstellungen sowie 
die Maximum-Entropie-Rekonstruktion von Faserorientierungsvertei- 
lungen durch führende Faserorientierungstensoren untersucht. Eine 
detaillierte Analyse des Advani-Tucker-Orientierungsmittelwertes gibt 
Aufschluss über die Struktur verschiedener orientierungsgemittelter 
Meanfield-Homogenisierungen. Ein vollstandiges Bild der Richtungs- 
abhängigkeit der effektiven linear-elastischen Steifigkeiten, die durch 
diese Homogenisierungsverfahren vorhergesagt werden, wird beschrie- 
ben. Dazu wird eine innovative Visualisierungsmethode basierend auf 
Übersichtsplots und Polarplot-Ensembles entwickelt. Diese Methode 
ist anwendbar auf jede richtungsabhángige skalare Gróf?e, die eine 
Funktion verschiedener planarer Faserorientierungstensoren ist. 

Diese Arbeit motiviert die Verwendung von Faserorientierungsinfor- 
mationen vierter Stufe und hat somit weitreichende Implikationen für 
virtuelle Prozessketten für planare langfaserverstárkte Bauteile. 
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Summary 


Effective mechanical properties of fiber reinforced composites strongly 
depend on the microstructure, i.e., the fibers’ alignment within the 
composite. In this work, the influence of fiber orientation distribution on 
the effective properties of long-fiber-reinforced composites is studied. 
The present work is funded by the German Research Foundation (DFG) 
within the interdisciplinary research training group GRK 2078. Within 
this research group, the material class sheet molding compound (SMC) 
is of primary interest. SMC is a composite of long glass fibers and a 
thermoset polymer matrix, serving lightweight applications with flat 
component designs, for instance in the automotive industry. 

Averaged information of a fiber orientation distribution can be quantified 
by fiber orientation tensors. This work gives a complete picture of the 
influence of fiber orientation tensors up to fourth order on effective me- 
chanical properties, based on the newly derived variety of fourth-order 
fiber orientation tensors. This variety is identified by linear invariant 
decomposition of fourth-order fiber orientation tensors leading to simple 
and flexible parameterizations. Constraints of material symmetry are 
incorporated into these parameterizations and admissible parameter 
ranges are identified by demanding positive semi-definiteness. Through- 
out this work, closure approximations are contrasted by the variety 
of fourth-order fiber orientation tensors. Microstructures of SMC are 
approximated by planar orientation states. The correspondence of planar 
fiber distributions and planar fiber orientation tensors are presented in 
detail, deploying tensorial Fourier series in 2D- and 3D-frameworks 
as well as maximum entropy reconstruction of fiber orientation distri- 
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Summary 


butions by leading fiber orientation tensors. Detailed analysis of the 
Advani-Tucker orientation average yields insights on the structure of 
several orientation-averaging mean field homogenization approxima- 
tions. A complete picture of the directional dependence of effective 
linear elastic stiffnesses predicted by those mean field homogenizations 
is drawn. Therefore, an innovative visualization method based on 
overview plots and polar plot ensembles is derived. This method 
is applicable to any direction-dependent scalar quantity, which is a 
function of distinct planar fiber orientation tensors. 

This work emphasizes the importance of fourth-order fiber orientation 
information and therefore has far-reaching implications for virtual pro- 
cess chains for planar long-fiber-reinforced components. 
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Chapter 1 


Introduction 


1.1 Motivation 


Fiber-reinforced composites solve lightweight construction tasks in a 
wide variety of application areas. A major reason for the frequent 
use of these composites is their excellent mechanical performance at 
relatively low density. In the case of epoxy resins reinforced with 
glass fibers, competitive mechanical properties of the composite are 
achieved by the combination of the high strength and stiffness of the 
fibers with the flexibility of the matrix material. The microstructure 
of the composite, i.e., the amount, arrangement and orientation of the 
fibers, largely determines the effective properties. The realization of 
beneficial microstructures and, in particular, the suitable alignment of 
fibers in highly stressed areas of the component is therefore a crucial 
task in the design and manufacturing process. Optimal component 
design and optimized manufacturing requires a profound understand- 
ing of the relationship between fiber orientations and effective local 
mechanical material properties. The local mechanical properties can 
be determined by material models within the framework of continuum 
mechanics, on scales relevant to engineering. A practically relevant class 
of such material models are orientation-averaging linear elastic mean 
field homogenization methods, which are the focus of this work. These 
homogenization methods combine information on the microstructure 
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with mechanical properties of the constituents, fiber and matrix, to 
identify the effective material behavior. Fiber orientation tensors are 
well-established descriptors for the orientation of fibers on a selected 
discretization scale. They can be predicted numerically by process simu- 
lations as well as determined experimentally by computer tomography. 
The variety of second-order fiber orientation tensors is known and 
is represented graphically by the so-called fiber orientation triangle. 
Knowledge of the variety of all possible second-order fiber orientation 
tensors offers far-reaching simplification options for individual steps 
within a virtual process chain and facilitates the assessment and classifi- 
cation of microstructures. However, since the effective material behavior 
of solids depends on both second- and fourth-order fiber orientation 
tensors, the knowledge of the variety of fourth-order fiber orientation 
tensors is of great interest to obtain a complete picture of the realizable 
effective properties. The present work has been funded by the German 
Research Foundation (DFG) interdisciplinary research training group 
GRK 2078. This research training group focuses on the material class 
sheet molding compound (SMC), a composite of long glass fibers and 
thermoset polymer matrix, which is suitable for the mass production of 
flat components. In this work the variety of fourth-order fiber orientation 
tensors is studied and for the first time the complete picture of the 
realizable effective properties using microstructure properties in the 
form of planar fourth-order fiber orientation tensors is outlined. The 
subgroup of planar fourth-order fiber orientation tensors is of special 
interest, as the typical fiber lengths in SMC are in the range of two to 
ten times of the component thickness. Due to this length ratio and low 
tendency towards fiber bending, SMC microstructures may be described 
by a low-dimensional planar approximation. 
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The main objective of this thesis is the connection between local fiber ori- 
entation states in SMC components and effective mechanical properties. 
A compact and nevertheless complete representation of this connections 
is aimed at, which is in contrast to limited examinations of concrete 
examples, present in the literature. The identification of the variety of 
fiber orientation tensors of fourth-order represents the starting point for 
the investigations. This identification is done for general orientation 
states. In a subsequent step, subspaces of these general orientation 
states, which state reasonable approximations to SMC microstructures 
are of interest, possibly leading to compact representations. For the 
special case of planar orientation states, correlations between orientation 
tensors and orientation distributions are of concern. Therefore, methods 
reconstructing fiber orientation distributions by a set of leading fiber 
orientation tensors are given. These methods lead to a flexible numerical 
scheme for the calculation of orientation averages. Consequences of 
orientation states on effective mechanical properties are to be drawn 
for different orientation-averaging homogenizations. Visual represen- 
tations for all realizable planar orientation states are of special value 
giving a simple picture of the dependencies between microstructure and 
effective mechanical properties. Emphasis is to be taken on linearity of 
effective stiffnesses in fiber orientation tensors within the Advani-Tucker 
orientation average. 


The diversity of fiber orientation tensors is in direct contradiction to 
commonly used closure algorithms. Therefore, representations which 
contrast fourth-order information obtained by closure approximations 
and the variety of possible fourth-order information, are of interest. 
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The major novelties of this work, which is based on Bauer and Bóhlke 

(2022c;a;b), are summarized. 

* Aninvariant framework for parameterizations of fourth-order fiber 
orientation tensors based on deviators from the isotropic orientation 
state is developed. The most general case is obtained by an orthotropic 
deviator of second order and a triclinic deviator of fourth order. 
Within this framework, deviators of stronger material symmetries 
are obtained as special cases. 

* Explicit parameterizations and admissible parameter ranges for or- 
thotropic as well as planar fiber orientation tensors of fourth order are 
identified and visualized in Figures 4.5 to 4.8 and 5.1b. 

* Boundaries of the set of admissible fiber orientation tensors of fourth 
order are identified to correspond to localized fiber orientation states 
which can be represented by empirical fiber orientation distribution 
functions with a limited number of Dirac-distributions. 

* Redundancies and degeneration in the planar set of admissible fourth- 
order fiber orientation tensors are discussed in detail. 

* A generic visual representation of the dependence of any scalar quan- 
tity on planar fourth-order fiber orientation tensors is derived and 
applied to fiber orientation distribution functions and effective me- 
chanical properties. 

* Drawbacks of fiber orientation distribution reconstruction based on 
truncated Fourier series expansion within two- and three-dimensional 
frameworks are demonstrated, highlighting the role of the isotropic 
reference states. 

* For the planar case and within the two dimensional framework it is 
shown that interference of second- and fourth-order contributions 
leads to the variety of reconstructed FODF based on truncated Fourier 


series. 
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The maximum entropy reconstruction of Müller and Bóhlke (2016) is 
recast for the planar case in a two-dimensional framework, leading to 
a low-dimensional optimization problem. 

Fiber orientations of sheet molding compound are identified to be 
approximately planar, leading to a reduction of the independent 
components of fourth-order fiber orientation tensors from eleven 
to three in the orientation coordinate system. 

The orientation average of an elasticity tensor based on a fiber ori- 
entation tensor following Advani and Tucker III (1987) is linear in 
both the transversely isotropic elasticity tensor and the second- and 
fourth-order fiber orientation tensors. 

A new numerical formulation of the Advani-Tucker orientation av- 
erage in fiber orientation tensors Advani and Tucker III (1987) for 
the special class of planar fiber orientations is proposed based on 
a maximum entropy reconstruction of fiber orientation distribution 
functions. 

Effective stiffnesses obtained by homogenizations which perform 
the orientation average in the stiffness domain and those which 
perform the average in the compliance domain are demonstrated 
to differ significantly with respect to the orientation dependence of 
the generalized bulk modulus. 

Possible directional dependence of the elastic response for Advani- 
Tucker averaged two-phase materials of isotropic constituents and 
planar orientation measures are comprehensively presented. This 
presentation is complete and can be used to express the orientation 
dependence of any quantity which is a function of planar fiber orien- 
tation tensors up to fourth order. The dependence is restricted by the 
limited averaged information given by fourth-order fiber orientation 
tensors and due to the constraints of linear elasticity. 


1 Introduction 


1.4 State of the art’ 


1.4.1 Fiber orientation tensors 


Local fiber orientation distributions are completely defined by a fiber 
orientation distribution function (FODF). However, in applications the 
FODF is usually approximated by a limited number of fiber orientation 
tensors (FOT), as tensor representations fit into the continuum mechanics 
framework. Multiple kinds of orientation tensors exist, originating 
from Kanatani (1984). Applications related to fiber orientations include 
both identification of the fiber orientation information by experimental 
methods (Pinter et al., 2018) or process simulation (Meyer et al., 2020a) 
and usage of the orientation information in fullfield (Górthofer et al., 
2020), mean field (Müller and Bóhlke, 2016; Kehrer et al., 2020), or 
damage simulations (Schemmann et al., 2018b). Both groups of applica- 
tions, identification and usage of orientation information, benefit from 
a well known variety of fiber orientation tensors. If the variety of fiber 
orientation tensors is well-defined, the identified orientation tensors can 
be assessed and the space of input parameters of methods predicting 
effective mechanical properties is set. The variety of second-order fiber 
orientation tensors is well known (Cintra Jr and Tucker III, 1995; Chung 
and Kwon, 2002; Cowin, 1985) and used, e.g., in Kóbler et al. (2018); 
Górthofer et al. (2020); Kóbler et al. (2021) to generate clear pictures of 
the dependence of mechanical properties on the second-order directional 
measures. However, concise descriptions and parameterizations of 
fourth-order fiber orientation tensors are rare. The variety of transversely 
isotropic fourth-order orientation tensors is identified by Nomura et al. 
(1970) and used in Müller and Bóhlke (2016); Jack and Smith (2008). For 
a given second-order orientation tensor, closure approximations (Advani 
and Tucker III, 1990; Han and Im, 1999; Cintra Jr and Tucker III, 1995; 


1 Most of the content of this section is taken directly from Bauer and Böhlke (2022c;a;b). 
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Chung and Kwon, 2002; Montgomery-Smith et al., 2011a) identify a 
corresponding fourth-order orientation tensor based on assumptions. 
Although literature on closure approximations involves much infor- 
mation on fourth-order fiber orientation tensors, most closures state 
simplifying assumptions on the material symmetry of fourth-order FOT. 
In consequence, generic, i.e., triclinic FOT are not covered in the literature 


on closure approximations. 


1.4.2 Fiber orientation distribution reconstruction 


Fiber orientation distribution functions represent exact microstructure 
descriptors of the orientation of axisymmetric fibers within a specified 
volume of a fiber-reinforced composite (Kanatani, 1984; Advani and 
Tucker III, 1987). However, in practice (Bohlke et al., 2019; Górthofer 
et al., 2019) the exact distribution of fibers is commonly approximated 
by fiber orientation tensors which represent averaged directional mea- 
sures and can be directly determined either by non-destructive analysis 
methods, such as computer tomography (Pinter et al., 2018; Schöttl et al., 
2020) or by numerical flow simulations (Meyer et al., 2020a; Bóhlke 
et al., 2019; Karl et al., 2021a). Some applications, e.g., damage model- 
ing (Schemmann et al., 2018b) or averaging schemes (Hessman et al., 
2021; Brylka, 2017), require FODF information and therefore motivate 
reconstruction of FODF from a given set of leading fiber orientation 
tensors. Due to the averaged character of fiber orientation tensors, no 
one-to-one correspondence between a set of leading FOT and a FODF 
exists. Truncated Fourier series are used by Schóttl et al. (2020); Jack 
and Smith (2004); Eik et al. (2016) to identify an FODF based on leading 
FOT. However, the identified functions do not meet the non-negativity 
requirement of a FODF. Non-negativity of the reconstructed FODF 
based on transversely isotropic leading FOT of order two and four, 
is enforced by Müller and Bóhlke (2016) using a maximum entropy 
method. Maximum entropy methods are deployed, e.g., in Bóhlke 
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(2005; 2006). Closure approximations (Advani and Tucker III, 1990; 
Chung and Kwon, 2001; 2002; Cintra Jr and Tucker III, 1995; Han and 
Im, 1999; Jack et al., 2010; Kuzmin, 2018; Montgomery-Smith et al., 
2011a;b; Karl et al., 2021c; Dray et al., 2007) are commonly used to 
estimate fourth-order FOT based on second-order FOT and therefore 
have reconstruction character. The closures in Montgomery-Smith et al. 
(2011a;b) are based on assumptions on an FODF associated with the 
second- and fourth-order FOT of the closure and therefore state a FODF 
reconstruction scheme. Breuer et al. (2019) assess the reconstruction of 
FODFs based on closure approximations and minimum entropy method 
used by Bóhlke (2005); Müller and Bóhlke (2016). 


1.4.3 Orientation-averaging mean field homogenization 


Approximation of effective mechanical properties of heterogeneous 
materials by mean field homogenization dates back to Voigt (1889); 
Reuß (1929); Hill (1952) with major contributions by Hashin and Shtrik- 
man (1962); Walpole (1966a;b); Willis (1977; 1981); Walpole (1969); Mori 
and Tanaka (1973) with a strong dependence on Eshelby (1957). The 
transformation of volume averages into averages over orientations (Bóh- 
Ike, 2001) leads to orientation-averaging homogenizations (Benveniste, 
1987; Duschlbauer et al., 2006; Hessman et al., 2021; Müller and Bóhlke, 
2016; Müller, 2016; Brylka, 2017; Schemmann et al., 2018b; Kehrer et al., 
2020). Hessman et al. (2021) cast several popular orientation-averaging 
homogenizations (Hill, 1965b; Budiansky, 1965; Mori and Tanaka, 1973; 
Benveniste, 1987; Castafieda and Willis, 1995; Zheng and Du, 2001; 
Pierard et al., 2004) into a unified framework. Effective mechanical 
properties of SMC are identified by Trauth et al. (2021); Schemmann et al. 
(2018b); Kehrer (2019); Kehrer et al. (2020) using orientation-averaging 
mean field homogenization, whereas Górthofer et al. (2020) deploy a fast 
Fourier transform method (Schneider, 2021; Moulinec and Suquet, 1994; 
1998). Most mean field homogenizations are based on the averaging 
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scheme of Advani and Tucker III (1987) which is directly formulated 
in FOT. However, some authors apply direct numerical integration 
based on FODF, e.g., Pettermann et al. (1997). Müller and Bohlke 
(2016) showed that the identification of effective mechanical properties 
based on second-order orientation tensors yields insufficient accuracy 
(Breuer et al., 2019). Closure approximations are used by, e.g., Jack 
and Smith (2007; 2008); Buck et al. (2015); Müller and Bóhlke (2016); 
Goldberg et al. (2017) to model mechanical properties. However, as 
closure approximations are based upon assumptions and do not generate 
information, the assessment of Müller and Bóhlke (2016) applies to those 
studies. 


1.5 Outline 


The present work starts with a compact introduction to selected topics 
of continuum mechanics, i.e., kinematics, stresses, balance equations 
and the role of material modeling, followed by tensor algebra with focus 
on linear invariant decompositions, irreducible tensors and material 
symmetry in Chapter 3. 

The variety of fiber orientation tensors of order two and four is investi- 
gated in Chapter 4, which is based on Bauer and Bóhlke (2022c). Positive 
semi-definiteness is discovered to define the set of admissible parameter 
ranges for parameterizations of fourth-order fiber orientation tensors. 
Constraints induced by material symmetry are utilized explicitly to 
derive admissible orientation tensors for the transversely isotropic and 
planar case. The orthotropic case is presented visually. Notes on closure 
approximations and their limitations close this chapter. 

In Chapter 5, which is based on Bauer and Böhlke (20222), the in- 
sights on admissible planar fiber orientation tensors are refined and 
the correspondence with fiber distributions is investigated. For the 
first time, a complete picture of planar fiber orientation distributions, 
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reconstructed from the planar subspace of fiber orientation tensors, is 
given. This picture is based on a visualization method which slices 
the space of admissible fiber orientation tensors. Reconstruction of 
fiber orientation distribution functions by truncated tensorial Fourier 
series demonstrates advantage of a two-dimensional framework for 
planar orientations. The variety of orientation distributions is explicitly 
explained by interference of second- and fourth-order contributions 
on fiber orientation distributions. A maximum entropy reconstruction 
method for planar fiber distributions is given in compact and explicit 
notation in a two-dimensional framework. An explicit representation of 
the exact closure for the planar case closes this chapter. 

The transition from microstructure descriptors to mechanical properties 
is investigated in Chapter 6 based on Bauer and Bóhlke (2022b). Local 
fiber orientation states of SMC components are identified as approx- 
imately planar. Based on this planarity assumption, the variety of 
associated fiber orientation tensors and distributions is known. In 
consequence, the complete picture of the effective mechanical prop- 
erties obtained by orientation-averaging mean field homogenizations is 
presented. Characteristics of homogenization methods averaging in the 
stiffness or compliance domain are discussed. The orientation average of 
Advani-Tucker Advani and Tucker III (1987) is recast as linear invariant 
decomposition and a new numerical formulation based on reconstructed 
fiber distributions is given. 

Finally, in Chapter 7 the results are summarized and an outlook on 
future research is given. 
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1.6 Notation, frequently used acronyms, sym- 
bols, and operators 


Symbolic tensor notation is preferred in this work. Tensors of first 
order are denoted by bold lowercase letters such as q, n, v, e. Ten- 
sors of second-order are denoted by bold uppercase letters like N or 


Q and fourth-order tensors are denoted by, e.g., N or D. Tensors in 


representations for varying tensor order are represented by, e.g., Dik), 


where k defines the tensor order. Tensorial quantities are defined in a 
three-dimensional space, unless underlined. If a tensor is underlined 
like, e.g., q or N, it is defined in a two-dimensional space and follows this 
spaces algebra. A linear mapping of a second-order by a fourth-order 
tensor reads as A = C[B]. The scalar product reads as A - B. The 
tensor power, i.e., the k-th dyadic product of, e.g., a first order tensor 
a is denoted by a®* yielding, e.g., a?? = a & a & a. An orthonormal 
basis is denoted by {e;} with e; : e; = ó;; and the Kronecker delta ó;;. If 
a matrix of tensor coefficients is used in mixed notation, the coefficient 
matrix is directly followed by the tensor basis where the first index of the 
basis corresponds to the rows of the coefficients matrix, the second one 
to the columns. Summation convention applies, unless otherwise stated. 
Representations in index notation always refer to an orthonormal basis. 
The Rayleigh product is used to represent an active rotation of a physical 
quantity and for a first order tensor is defined by Q x n = n; Q e;. Sets, 
i.e., collections of quantities, are denoted by calligraphic symbols, e.g., 
F and are constructed by curly braces. Inside the curly braces, elements 
of the set are given explicitly, or by a generator expression following the 
pattern {quantity | condition fulfilled inside set}. Although this work 
and related code is based on Harris et al. (2020); Meurer et al. (2017), 
numbering and indices follow the continuum mechanics convention 


starting at one. 
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Acronyms 

CODF Crystal orientation distribution function 
CT Computer tomography 

FODF Fiber orientation distribution function 
FOT Fiber orientation tensor 

IBOF Invariant-based optimal fitting 

iso Isotropic 

ORF Orthotropic fitted 

ortho Orthotropic 

perm Permutation 

plan Planar 

RV Reference volume 

RVE Representative volume element 

SMC Sheet molding compound 

sym Symmetric 

transv Transversely isotropic 

unidirect Unidirectional 


Latin letters 


a,b,c, A, B, C Scalar quantities 

a,b,c Tensors of first order 

A,B,C Tensors of second order 

A,B,C Tensors of fourth order 

Ar), Bary, Ct Tensors of k-th order 

Cf Fiber volume content 

e Mass specific internal energy 
w* Complementary energy density 
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1.6 Latin letters 


= % 
2 


Z*E top 


(CMrTlinearStiffness 
CMrTlinearCompliance 
(CMTOAB 

F 


I 
JS 


Direction-dependent Young’s modulus 

Shear modulus 

Compression modulus 

Direction-dependent generalized compression 
modulus 

Singular surface 

Reference volume 

Normal vector, i.e., direction 

Heat flux 

Stress vector 

Displacement 

Position in current configuration 

Position in reference configuration 

Distortion rate tensor 

Structure tensor 

Deformation gradient 

Displacement gradient 

Second-order fiber orientation tensor of Kanatani 
first kind 

Rotation, i.e., element of SO(3) 

Strain localization tensor of the i-th phase in the 
single inclusion problem 

Elastic stiffness 

Reference stiffness 

Effective stiffnesss of two-step Hashin Shtrikman 
homogenization scheme 

Average of phase stiffnesses 

Inverse of the average of phase compliances 
Effective stiffnesss following Benveniste (1987) 
Structure tensor 

Identity on fourth-order tensors 

Identity on symmetric fourth-order tensors 
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Greek letters 


> S53 


ES 


qa ar & v 


Operators 


Fourth-order fiber orientation tensor of Kanatani 
first kind 

Isotropic projectors 

Hill's polarization tensor 

Fiber orientation tensor of Kanatani third kind 
and k-th order 


Weight of the i-th transversely isotropic structure 
tensor 

Entropy 

Temperature 

i-th eigenvalue of N 

Density 

Fiber orientation distribution function 
Movement 

Infinitesimal strain tensor 

Cauchy stress 


Determinant of a quantity 
Deviatoric part of a quantity 
Symmetric part of a quantity 
Material time derivative 
Symmetry group of a quantity 
Jump of a field quantity 
Quantity in 2D framework 


1.6 Operators 


k-th dyadic product of a quantity 
Transpose of a second-order tensor 


Right transpose of a fourth-order tensor 
Rayleigh product operator 
Orthonormal coordinate system 
Orthonormal eigensystem of N 
Quantity at a singular surface 
Fiber quantity 

Matrix quantity 

Effective quantity 

Volume average over RVE 
Generalized average 
Advani-Tucker average 


Advani-Tucker average explicitly formulated in 


fiber orientation tensor N 


Advani-Tucker average in formulation of Gold- 


berg et al. (2017) 


Advani-Tucker average in numerical 2D maxi- 


mum entropy formulation 
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Chapter 2 


Continuum Mechanics 


2.1 Motivation 


Abstracted and restricted to a classical mechanical point of view, motion 
and deformation of physical bodies in an Euclidean space at different 
size-scales can be used to describe the environment we are living in. 
Continuum mechanics allows for a mathematically well-defined de- 
scription of bodies, their motion and their deformation on a chosen 
size scale. The flexibility in choosing an appropriate size scale for engi- 
neering applications is a striking feature of continuum mechanics. For 
engineering tasks on macroscopic technical components, a consideration 
on the scale of molecules or atoms is not purposeful. For this reason, 
numerical implementations based on models, motivated by continuum 
mechanics, are ubiquitous in engineering. The best known application is 
structural analysis using the finite element method, to analyze stresses 
and deformations in technical components. The present work is focused 
on the analysis of microstructures, i.e., the structure of heterogeneous 
components, specifically fiber-reinforced polymers, and the effect of the 
microstructure on macroscopic material properties. Due to this focus, 
only a very small part of the continuum mechanics fundamentals is 
relevant for this work. In the following chapter, important basic terms 
are introduced and a classification in the broad context of continuum 
mechanics is given. This chapter in large parts follows Bóhlke and 
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Frohnapfel (2019); Bóhlke (2014) and is influenced by Bertram and Glüge 
(2013); Bertram (2012). 


2.2 Kinematics and deformation 


We choose a closed region in a three-dimensional Euclidean space at 
a randomly chosen reference time and call this region material body. 
Points of a material body in the initial placement, called reference place- 
ment, can be identified by a position vector X. The current placement of 
a material point is denoted by x. The movement of a material point is 
described by a function 


x (X,t) =x (X,t), (2.1) 


taken to be invertible and twice continuously differentiable. Therefore, 
x (X,t) at a given time t, defines a field specifying the current placement 
of a body composed of material points. The first partial derivative of the 
movement x (X, t) with respect to the reference placement 


(2.2) 


is called deformation gradient and states the basis for several possible 
deformation measures capable of describing large deformations. How- 
ever, in this work only small deformations are considered, which can be 
described by a linearized theory, developed hereafter. The displacement 
field is defined as 

u(X,t) = x (X,t) - X. (2.3) 


Taking the first partial derivative of the displacement field u (X, t) with 
respect to the reference placement X, the displacement gradient is 
defined as 

. Ou(X,t) 


H (X,1) = — (2.4) 


18 


2.3 Time derivatives 


Based on the displacement gradient, a unique deformation measure for 
small deformations exists and is given by the infinitesimal strain tensor 
€ (x, t), which can be calculated for a known field of displacement u (x, t) 


by 
«o i (EE (hem 


= l (H (x,t) + H7 (x,t). (2.5) 


The transposition AT of a second-order tensor A is defined by 
v:-(ATw)=w-(Av) Vv,w (2.6) 


with tensors of first order v and w. Theories, based on the infinitesimal 
strain tensor, which is linear in the deformation gradient, are called 
geometrically linear. Deformations are small, if ||H|| = VH-H « 1 
holds and imply x « X. 


2.3 Time derivatives 


The material time derivative of a field quantity Y(X, t) is given by 


X,t) = —— : 
ix) = I, Q7) 
whereas the material time derivative of the same quantity parameterized 
as a function of the current placement y(x, t) is given by 


(x,t) = ents) m u Dt) 


(2.8) 
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following the chain rule with the velocity of a material point 


-— Ox(X, t) 


m (2.9) 


2.4 Stress 


Mechanical interaction between bodies can be described by the concept 
of forces. Similarly, mechanical interaction between material points 
located inside or at the boundary of a body can be described by force 
densities. The introduction of the principle of cuts (Bertram and Glüge, 
2013) allows the definition of stress tensors. Let a body be mentally 
divided into two new bodies along a cut surface. The mechanical 
interaction that existed before the hypothetical separation process can 
be represented by a stress vector t on each of the two, newly created cut 
surfaces. The stress vectors on both surfaces are of equal magnitude and 
oriented in opposite directions. According to Cauchy's theorem, a stress 
vector on the cut surface with normal n is defined by 


t—on (2.10) 


and is associated with a symmetric second-order tensor, the Cauchy 
stress ø. The Cauchy stress can be interpreted as a tensor-valued field 
o (X) which assigns each point X on a material body a local stress tensor. 


2. Balance equations 


The object of interest is an abstract material body which consists of 
material points. As this body might move in space, the volume occupied 
by the body Y(t) might be a function of time. Points at which fields 
are smooth are called regular points. Smooth means continuous and 
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continuously differentiable. This smoothness of the fields is not given 
perpendicular to so-called singular surfaces S(t). Singular surfaces 
(Bertram, 2012; Cermelli et al., 2005; Prahs, 2020) describe places where 
jumps of field quantities occur and separate areas of regular points. An 
example of singular surfaces are the interfaces of different phases in 
inhomogeneous materials, e.g., the interface between fibers and matrix 
in fiber reinforced composites. 


2.5.1 Divergence theorem 


For a given vector-valued field quantity b and a domain V(t) with 
boundary V(t) and singular surfaces S(t), the surface integral over b 
expressed by a volume integral is 


/ b-ndA— 2 [b]: ndA = / div(b) dV (2.11) 
ov(t) S(t) V(t) 


with the jump of b at the singular surface [b] = bt — b^ and the normal 
n at S(t) pointing from the side where b has the value b~ towards 
the side where b is equal to bt. The divergence operator div (-) = 
tr (grad (-)) is based on the trace operator tr (-) and the gradient denoted 
by grad (-). This relation is essential for transforming integral balances 
into local balances. 


2.5.2 Transport theorem 


Let Y (t) be the volume of a certain domain, for example a body. This 
volume may change and therefore is a function of time. The time 
derivative of an integral over a domain depends both on the change 
of the integrated field quantity v and the change of the volume Y (t) 
occupied by the domain. For a volume V associated with a material, this 
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relation is described by 


m V= oe dV + uv -ndA 
dt Jve) va) Ot av(t) 
= ) [vv dA (2.12) 
S(t) 


with v being the velocity field of the body, [y] the jump of v at S(t) and 
v$ the velocity of the singular surface normal to S(t). In the absense of 
singular surfaces, combining Equation (2.12) with Equation (2.11) yields 


d 7 ap... 
ái f wo f o DF + div (yv) dV (2.13) 


being valid only at regular points. 


2.5.3 General balance equation 


We consider balancing a time-dependent field quantity v (x,t) over a 

material body. The volume occupied by the material body Y (t) is a 

function of time t. The field quantity is said to be additive in nature and 

has an associated 

e Production density py (x,t) representing the local production of y per 
unit volume 

* Supply density s; (x,t) representing the local supply of v per unit 
volume 

e Non-convective flux qy (x,t) representing the local flux of v per unit 
area. 

Fluxes caused by other physical processes than convection are called 

non-convective fluxes. Heat flux and Cauchy stress are examples of 
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non-convective fluxes. The balance of 7 over V reads as 


al vdV = / po + Sy dV + | qy -ndA. (2.14) 
dt Jy) v(t) av(t) 


The local form of this balance restricted to regular points is obtained 
using Equation (2.13) leading to 
OW 


Bt + div (dv) = py + sy + div (qy) (2.15) 


with wv being the convective flux. If the volume in which Equation (2.14) 
is evaluated is completely shrunk to the singular surfaces, a jump 
condition for the field quantity at S(t) is obtained and reads 


[V(v:n —v$)] - [ay]: n= 0. (2.16) 


In the case of a material singular surface, such as those present at phase 


boundaries in heterogeneous materials, v? — v holds and the jump 


condition in Equation (2.16) reduces to a requirement on qy- 


2.5.4 Specific balance equations at regular points 


Balance equations for mass, linear momentum, angular momentum, en- 
ergy, and entropy can be obtained from Equation (2.14) or Equation (2.15) 
using the quantities specified in Table 2.1. The remaining considerations 
within this chapter are restricted to regular points. Jump conditions at 
singular surfaces can be found in Bóhlke (2014); Bertram (2012). 


Mass 


The balance of mass is given by 


p+ pdiv (v) = 0 (2.17) 
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Balance Quantity: 4% | Produc.: py Supply: sy Flux: qy 
Mass p 0 0 0 
Lin. Momentum pv 0 pb c 
Ang. Momentum || (x — xo) x pv O0 (x— xo) x pb|(x —x0) xa 
Energy ped &v-v 0| pw-cpb.v| -q+o'v 
Entropy pn PPn zw = ad 


Table 2.1: Building blocks to derive specific balance equations from the general balance 
equation. The ansatz for supply and flux of entropy are of constitute nature. The 
abbreviation produc. stands for production. 


with the mass density p. 


Linear momentum 


The axiom of the balance of linear momentum fits experimental observa- 
tions for non-relativistic cases and leads to 


pv = pb + div (o). (2.18) 
with a mass force density b. 
Angular momentum 


The balance of angular momentum, is an axiom, which for a Boltzmann 
continuum, i.e., a continuum without point-wise rotational degrees of 
freedom, results in a system of algebraic equations denoted by 


g —3!. (2.19) 
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Total energy 


Total energy is the sum of kinetic and internal energy and its balance is 


pe + m (v. v) = pw — div (q) + pb - v + div (oTv) (2.20) 
with the heat source density w and the heat flux q. The total energy is a 
conservation quantity because the production term vanishes. The local 
balances used in this section apply to material points. A material point 
can be considered as a thermodynamic system. The energy due to the 
motion of this system is its kinetic energy. The mass specific internal 
energy e contains the energy which is independent of the system’s 
motion. The internal energy states a thermodynamic potential. 


Kinetic energy 


The balance of kinetic energy is obtained by scalar multiplication of 
Equation (2.18) with the velocity field v leading to 


5p(v-v) = pb- v div (atv) —o-D (2.21) 
with the distortion rate tensor D = sym(grad (v)). This procedure 


implies that the balance of kinetic energy is another representation of 
information already included in the balance of linear momentum. 
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Internal energy 


The balance of internal energy, given in Equation (2.22), is obtained by 
subtracting the balance of kinetic energy from the balance of total energy. 


perso (v. v) =pw — div (q) + pb: v + div (oTv) tot. (2.20) 

j^ (v.v) = pb:v-rdiv(cTv) —e-D kin. (2.21) 

pe =pw — div (q) +oa-D (2.22) 
Entropy 


The balance of entropy is given by 


. pw X. (d 
pi = = div ( 2) Hep (2.23) 


with entropy 7, temperature 0 and the entropy production p,, which has 
to be non-negative to satisfy the second law of thermodynamics. 


2.6 Material modeling and closure of balance 
equations 


The following considerations are limited to small deformations and 
in consequence, deformation is described by the infinitesimal strain 
tensor e. Further consequences are that the reference placement is 
approximately equal to the current placement (X = x), the density is 
approximately constant (p ~ po) and the symmetric part of the velocity 
gradient can be approximated by the rate of e, i.e., D ~ & holds. 


Applying the restrictions of small deformations and regular points, the 
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resulting balance equations 


p ™ po (2.24) 
g —g! (2.25) 
pu = pb + div (c) (2.26) 
pé = pw — div(q) +0 -è (2.27) 
a Oe avid 
p — 7g div ( 3) + PPn (2.28) 


state the starting point to derive a closed system of field equations 
for initial boundary value problems. The quantities p, b and w are 
known for a given initial boundary value problem. Displacement u, 
temperature 0 and a vector of internal variables a are selected as base 
fields B = {u, 0, a}. The base fields state the solution space of the initial 
boundary value problems. 

From Equations (2.26) and (2.27) follows that the remaining set of quan- 
tities is K = (0,6,q, n}. These quantities are called constitutive quan- 
tities and are unknown and not yet linked to the base fields. From 
an engineering point of view, the entropy is usually not of practical 
interest. However, including entropy as an auxiliary quantity in the set 
of constitutive quantities permits to incorporate restrictions stated by the 
second law of thermodynamics. In order to generate a closed system of 
equations, a set of material functions M has to be defined, expressing the 
constitutive quantities as functions of the base fields. As the constitutive 
quantities K may also depend on gradients and time derivatives of the 
base fields, a set of input quantities £ is defined containing the base 
fields as well as their gradients and derivatives. 

Material modeling is the process of defining the material functions M 
such that K = M(E) closes the system of Equations (2.24) to (2.28) while 
* Fulfilling restrictions stated by material theory 


* Not violating the second law of thermodynamics and 
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* Reproducing experimental observations. 

Equations (2.24) and (2.25) state algebraic equations and can be fulfilled 
identically. Equations (2.26) and (2.27) state the primary partial differen- 
tial equations and Equation (2.28) states the basis for restrictions caused 
by the second law of thermodynamics. 


2. Linear elasticity 


Within this work, only the simplest material model, linear elasticity for 
small deformations is considered. Linear elasticity implies, that there 
is no interaction between temperature and the displacement field. In 
addition, no state variables exist, as the current stress is solely deter- 
mined by the current deformation. In consequence, the only non-trivial 
field equation is the balance of linear momentum in Equation (2.26) and 
material modeling is reduced to finding the relation between o and e. 


Linear elasticity postulates this relation to be a linear mapping 
c = C [e] (2.29) 


defining the material stiffness C. As the stiffness maps between symmet- 
ric tensors of second-order, without loss of generality, the stiffness can 
be assumed to have both minor symmetries. A fourth-order tensor C is 
called minor symmetric, if 


A-C[B] = AT. C[BJ| A. C[B'] VA, B (2.30) 


holds with A and B being tensors of second order and their transposition 
is defined in Equation (2.6). Hyperelasticity, i.e., postulating the existence 
of a strain energy density function, implies the main symmetry of the 
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stiffness C. A fourth-order tensor C is major symmetric, if 
A-C[B]| 2B. C[A] VA, B (2.31) 


holds. Based on the assumptions of linear elasticity and small deforma- 
tions, the interaction between averaged microstructure descriptors and 
the linear elastic anisotropic effective material behavior is investigated 
within this work. 
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Chapter 3 


Tensoralgebra 


3.1 Motivation 


Euclidean tensors, i.e., tensor components 4;; associated with a fixed 
orthonormal coordinate system {e;}, defining A = A;je; Q ej, are 
established descriptors for physical quantities in continuum mechan- 
ics. Second- and fourth-order tensors have a special significance in 
continuum mechancis, since stress and strain are described by second- 
order tensors and linear constitutive laws can be represented as linear 
mappings between tensors of second order. Tensors of fourth order 
map second-order tensors linearly. A major contribution of this work is 
the identification of suitable parameterizations for special fourth-order 
tensors, namely fiber orientation tensors. The most important building 
blocks are linear invariant decompositions, irreducible tensors and 
material symmetry. Following Adams et al. (1992); Onat and Leckie 
(1988); Jerphagnon et al. (1978), material modeling can be reduced to 
the identification of irreducible tensors and definition of corresponding 
evolution equations. 


3.2 Linear invariant decomposition 


Linear invariant decomposition of symmetric second-order tensors is a 
common topic of continuum mechanics lessons. Applied to the stress 
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tensor, a linear invariant decomposition into the sum of a spherical and 
a deviatoric part reads as 


o= tr (o) I+ dev (o) (3.1) 


= —pI + dev (ø) (3.2) 


and can be motivated by the physical interpretation of the pressure p and 
remaining deviatoric stresses dev (c) = ø — itr (c) I with the identity 
on second-order tensors I defined by I [v] = v for all first-order tensors 
v. Both parts are orthogonal and transform differently with respect to 
rotations, i.e., actions of elements of the special orthogonal group SO(3). 
The spherical part is isotropic by definition, i.e., 


Qx(-pPD=-pI YQ € SO(3) (3.3) 


holds, whereas the deviatoric stress transforms with respect to Q € 
SO(3) by 


Qx dev (c) = Q dev (o) Q7}. (3.4) 


The scalar pressure and the deviatoric stress do not interact with each 
other under the action of elements of SO(3) and are called irreducible 
tensors, being of order zero and two. Linear invariant decompositions 
into irreducible tensors exist for tensors of arbitrary order and symmetry. 
However, in this work fourth-order tensors with at least minor and 
major symmetry are of special interest. 

Different linear invariant decompositions are given in the literature, 
see Rychlewski (2000); Forte and Vianello (1996); Olive et al. (2017). A 
detailed discussion on specific decompositions defined by the selec- 
tion of two pairs of symmetrization operators is given in Rychlewski 
(2000) and utilized in Section 4.6.1. However, the so-called harmonic 
decomposition is the most common decomposition (Forte and Vianello, 
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1996; Rychlewski, 2000; Bóhlke and Brüggemann, 2001; Fernández and 
Bóhlke, 2019; Lobos et al., 2017; Backus, 1970; Baerheim, 1993; Boehler 
et al., 1994; Mochizuki, 1988; Rychlewski, 2001; Lubarda and Chen, 2008) 
and decomposes a minor and major symmetric fourth-order tensor C 


into 
C (hi, ho, H;, H», H) ES 
hi Pi + h3 Po + (I8 Hi +H, 8I) 4 J [H5] + H (3.5) 
with 
P2 
hı = Pı s C, ha = 2" C, H = dev (C) (3.6) 
|[IPa | 
and 
1 
H: = 7 [5 dev (C) — 4 dev (V)] (3.7) 
1 
H, = 7 [3 dev (C) — 2 dev (V)] (3.8) 
based on the dilatational modulus C and the Voigt tensor V which are 
given by 
(C); = (C); (3.9) 
My = (C) iik; (3.10) 


see Cowin (1989). The isotropic projectors Pı and P» are given by 


1 
P; = 3l QI, P, =I — P; (3.11) 


with P= (I-II), I = e; 8 €j Q ei Q 6j. (3.12) 


1 
2 


The right index transpose operator for fourth-order tensors (-)™ is 
defined by (A ® B)'* = A & B! for arbitrary second-order tensors A, 
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B and the transpose of a second-order tensor (-)!. The deviatoric part 
of a fourth-order tensor is defined in Spencer (1970); Kanatani (1984); 
Bóhlke (2001) and can be written as 


dev (C) = C— sym(C @I)+ Ž C Isym(I6G I). (3.13) 


with the intermediate quantity C = C [I] and the operator sym(-) taking 
the completely symmetric part of a tensor (Spencer, 1970). The isotropic 
sixth-order tensor J in Equation (3.5) is defined by 


1 
J= 5 (ei 8 €; te; 8 ei) 
e (ej Q ek + ek Q ej) 
@ (ex @ ei + ei 8 ex), (3.14) 


see Forte and Vianello (1996); Böhlke and Brüggemann (2001). 
Consequently, any minor and major symmetric fourth-order tensor C can 
be decomposed into a unique set of irreducible tensors hı, ha, H4, H2, H 
of orders zero, zero, two, two and four. Following Forte and Vianello 
(1996, Equation (9)), the harmonic parts obtained from Equation (3.5) 
transform with respect to a rotation Q € SO(3) with 


Q*C = (hi, ha, Qx H,, Qx H5, Q«H). (3.15) 


33 Irreducible tensors 


Irreducible tensors are completely (index) symmetric and traceless 
(Spencer, 1970; Forte and Vianello, 1996; Jerphagnon et al., 1978; Adams 
et al., 1992; Cowin, 1989; Rychlewski, 2000). Due to the complete 
symmetry and the condition of the vanishing trace, the number of 
independent components of the tensors obtained by the decomposition 
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in Equation (3.5) are 


dim (hy) = 1, dim (A3) = 1, dim (H1) = 5, (3.16) 
dim (H5) = 5, dim (IH) — 9 (3.17) 
and reflect the dimension of the associated subspaces. In the absence of 


additional (material) symmetry a generic irreducible tensor of second 
order reads as 


ĝi îs ĝa 
F (di; d», d3, da, ds) = do da e; © ej. (3.18) 
sym —(d1 + d) 


and an irreducible tensor of fourth order is given by 


F (fi, ...f9) = (8.19) 
—(pit+ pe) pi Da —V2(pa+ ps) V2 Pe V2 fs 
-(Bı + Bs) Pa V2 pa —V2 (fs + bz) V2 Po 
—(p2 + Bs) | V2 Bs V2 pr — V2 (ps + po) 
completely symmetric 
Be & Be 


in Kelvin-Mandel notation and short hand notation for completely 
symmetric tensors defined in Appendix B.1 and derived in Section 4.6. 
Representations of irreducible tensors in Equations (3.18) and (3.19) 
combined with the harmonic decomposition in Equation (3.5) can be 
used to parameterize fourth-order tensors. 


3.4 Material symmetry 


The concept of material symmetries is an important tool for understand- 
ing consequences of the structured composition of matter (Jerphagnon 
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et al., 1978). Even in the absence of such a structured composition, artifi- 
cial assumption of material symmetry offers the possibility to reduce the 
dimensionality of a problem and to gain in this way simplified insights 
into underlying principles. 
A quantity C is symmetric with respect to a rotation, if the rotation does 
not affect the quantity, i.e., 


Q+C=C (3.20) 


holds. If Equation (3.20) is true for a specific Q, this rotation is said to be 
part of the symmetry group of C, denoted by S(C). Consequently, the 
symmetry group S(C) of a quantity C is defined by 


Q+C=C vQes(C) (3.21) 


Forte and Vianello (1996) showed that eight material symmetries, i.e., 
eight distinct closed groups of elements of SO(3), exist. These material 
symmetries are also referred to as material symmetry classes and their 
hierarchy is depicted in Fig. 3.1, following Bóna et al. (2007); Olive et al. 
(2022); Abramian et al. (2020); Weber et al. (2019). Zheng and Boehler 
(1994) discuss connections to the concept of point-groups. Visualizations 
of the material symmetry classes based on symmetry planes are given in 
Table 3.1 following Francois et al. (1998). 

Material symmetry of a quantity simplifies its parameterization by 
reducing the number of independent components. Following Forte 
and Vianello (1996, Equation (10)), the symmetry group of a tensorial 
quantity C is the intersection of the symmetry groups of its harmonic 
parts, i.e., 


S (C) = S (H1) N S (H2) NS (H). (3.22) 


Therefore, a physical quantity C might be parameterized by the combina- 
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isotropic 
transv. isotropic cubic 
tetragonal 
trigonal 


orthotropic 


monoclinic 


triclinic 


OO 


Figure 3.1: Hierarchy of material symmetries following Abramian et al. (2020); Francois 
et al. (1998). Arrows indicate inclusion, i.e., an arrow points from a symmetry classes with 
strong symmetry to a symmetry class which is included. 


tion of a linear invariant decomposition and a set of irreducible tensors 
potentially constrained by material symmetry. Essential building blocks 
for such a parameterization are irreducible second- and fourth-order 
tensors constrained by material symmetry. 


3.5 Irreducible tensors constrained by material 
symmetry 


This section is based on personal communication with Prof. Dr.-Ing. 
habil. Thomas Bóhlke, (Bóhlke, 2021). Derived from the generic parame- 
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terizations of irreducible tensors in Equations (3.18) and (3.19), material 
symmetries lead to the tensors listed below. 


3.5.1 Second-order tensors 


Triclinic, monoclinic and orthotropic irreducible second-order tensors 
pid, pmono, Fortho coincide and can not be distinguished, i.e., F = 
prid — pmono — porto holds. In other words, orthotropy is the weakest 
possible material symmetry an irreducible second-order tensor can have. 
This observation is directly related to the observation that any symmetric 
second-order tensor can be diagonalized. If an ordering convention for 
the eigenvectors based on eigenvalues applies, exactly four right-handed 
and four left-handed coordinate systems exist (Hasan et al., 2001; Bro 
et al., 2008). These are called eigensystems and lead to representations 
of Equation (3.18) with vanishing off-diagonal entries in the component 


matrix reading as 


F (q1, 92, Q (q3, qa, 95)) = (3.23) 
qi 0 0 
q 0 Qe;& Qe; 
sym —(q + q2) 


with an orthogonal tensor 


Q (qs. qa, 5) = Vi Q ei (3.24) 


defined by three scalars and mapping the arbitrary but fixed basis {e;} 
onto an eigensystem {v;}. This observation highlights the relevance of 
Equation (3.22), since the symmetry group of one single harmonic part 
does not define the symmetry group of the decomposed quantity. How- 
ever, eigensystems of a decomposed quantity, such as C in Equation (3.5) 
are also eigensystems of the harmonic parts. Hereafter, the dependence 
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of the eigensystem on usually three scalar parameters is no longer stated 
explicitly and the representations are given directly in an eigensystem 
te. 

A transversely isotropic irreducible tensor of second order in an eigen- 
system which first axis v; is aligned with the transversely isotropic axis 
is given by 


1 0 0 
iransv (q) =q —1/2 0 9; @V;. (3.25) 
sym 1/2 


The definition of a valid eigensystem in the transversely isotropic case 
only requires two scalars, since any rotation about the transversely 
isotropic axis leads to an equally admissible eigensystem. The represen- 
tation in Equation (3.25) coincides with those of trigonal and tetragonal 
irreducible tensors of second order, i.e., F'ansv = ptetra — ptrigo holds. 
Irreducible tensors of second order with isotropic or cubic material 
symmetry vanish, i.e., Fbic = piso = 0. 


3.5.2 Fourth-order tensors 


Departing from the generic, i.e., triclinic, case in Equation (3.19) a 
monoclinic irreducible tensor of fourth-order is given by 


pones (pi; P2, P3, P8, Po) = (3.26) 
—(pı +p2) pr pa 0 0 v2ps 
—(pi-c ps) P3 0 0 v2ps 
—(pa ps) | 0 0 -v2 (ps + po) 
completely symmetric 
B? 8 BY 
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in a coordinate system {v;} with v3 aligned with the normal of the 
monoclinic symmetry plane, see Table 3.1. The constraints of material 
symmetry for weaker symmetries lead to the irreducible tensors 


ortho (1, pa, p3) = (3.27) 
—(pitpe) m pa 0.0 0 
—(pi- ps) P3 0.0 0 
-(pe+ps) |0 0 0 ? e BY, 
completely symmetric 
FM8°*Y (p3, pg) = (3.28) 


8p3 —4p3 —4p3 0 0 0 
3p3 — pa 0 0 V2po9 
3p; |0 0 —V2pg 


B? @ Bi, 
completely symmetric 
p'ea (p1, pg) = (3.29) 
-2pı n n 000 
—(pitps) ps 000 
—(pi+p3)|0 0 0 ? e BY, 
completely symmetric 
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8 -4 -4/;0 0 0 
3 il 0 0 0 
3 |o 0 0 — 
transv x (pa) = ps B; & BY, (3.30) 


completely | symmetric 


and 
—2 1 1 0.0 0 
—2 1 0.0 0 
-2/0 0 0 


F™ (p1) = pi Y@Bi (31) 


completely | symmetric 


represented in coordinate systems aligned along characteristic axes of 
the material. Isotropic irreducible tensors vanish by definition. The 
superscripts of the irreducible tensors [Fono z, jetrigo xy, pietra x and IFtransv x 
indicate, that for these material symmetry groups multiple minimal 
representations by tensor components exist. These representation are 
equivalent, but differ by the mapping of the axes of the coordinate 
system {v;} and the charactersitic axes of the material, see Table 3.1. 
The representations in Equations (3.26) to (3.31) fit the alignment of the 
visualizations in Table 3.1. The aligment of the dominant axes of the 
trigonal, tetragonal and transversely isotropic representations with vı 
and the alignment of the symmetry plane normal of the monoclinic case 
with v3 proves to be practical later in this work. Representations of fiber 
orientation tensors of planar fiber-reinforced composites are directly 
connected to the selected representations if the out-of-plane normal of 
the plane of planarity is aligned with the local v3 axis and interpreted 
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as monoclinic symmetry plane. The fiber direction of unidirectional- 
reinforced composites is typically associated with either the v3 or vı axis. 
In this work, the latter convention is followed. 
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Material Symmetry Symmetry Planes | Generators 
Triclinic I 
Monoclinic e: 
Orthotropic MES 
Z| 
v3 E 
Trigonal a = esi “le 
p ^v | 
EH t» 
3 E 
Tetragonal amp, v2 or eus 
eae Vi 
Transversely isotropic a QS, 
i 37/2 27/3 
Cubic BL ow EN 
Isotropic Qi 


Table 3.1: Visualization of symmetry planes of material symmetry classes. The specified 
non-unique generators are rotations, i.e., elements of SO(3), which generate the complete 
symmetry group, if combined exhaustively (Weber et al., 2019). Two rotations Q1 and Q» 
are combined to a third one Qs by linear contraction, i.e., Q3 = Q1Q». The rotation QZ, 
rotates about the axis e3 by the angle m and the rotation Q7 rotates about the arbitrary 
axis v by the arbitrary but non-vanishing angle y. 
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Variety of Fiber Orientation 
Tensors ' 


4.1 Introduction 


The microstructure of fiber reinforced composites significantly affects the 
effective mechanical properties of the compound. Two local microstruc- 
ture descriptors commonly used to predict the effective properties are 
the fiber volume fraction and the fiber orientation distribution. The local 
fiber orientation distribution is completely defined by a fiber orientation 
distribution function (FODF). However, in practical applications the 
FODF is usually approximated by a limited number of fiber orientation 
tensors, as tensor representations fit into the continuum mechanics 
framework. Multiple kinds of orientation tensors exist (Kanatani, 1984). 
Practical applications related to fiber orientations include both identifica- 
tion of the fiber orientation information by experimental methods (Pinter 
et al., 2018) or process simulation (Meyer et al., 2020a) and usage of the 
orientation information in full field (Górthofer et al., 2020), mean field 
(Miiller and Bohlke, 2016; Kehrer et al., 2020) or damage simulations 
(Schemmann et al., 2018b). Both groups of applications, identification 


1 This chapter reproduces (Bauer and Bóhlke, 2022c), i.e., Bauer, J. K., Bóhlke, T., 2022. 
Variety of fiber orientation tensors. Mathematics and Mechanics of Solids 27 (7), 
1185-1211, 10.1177/10812865211057602. Reproduced with permission. ©2021 The 
Authors. Published by SAGE Publications Ltd under CC BY-NC 4.0 
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and usage of orientation information, benefit from a well known variety 
of fiber orientation tensors. If the variety of fiber orientation tensors is 
well defined, the identified orientation tensors can be assessed and the 
space of input parameters of methods predicting effective mechanical 
properties is set. The variety of second-order fiber orientation tensors 
is well known (Cintra Jr and Tucker III, 1995; Chung and Kwon, 2002; 
Cowin, 1985) and used, e.g., in Kóbler et al. (2018); Górthofer et al. 
(2020); Kóbler et al. (2021) to generate clear pictures of the dependence 
of mechanical properties on the second-order directional measures. 
However, concise descriptions and parameterizations of fourth-order 
fiber orientation tensors are rare. The variety of transversely isotropic 
fourth-order orientation tensors is identified by Nomura et al. (1970) 
and used in Müller and Bóhlke (2016); Jack and Smith (2008). Müller 
and Bóhlke (2016) showed that the identification of effective mechanical 
properties based on second-order orientation tensors yields insufficient 
accuracy (Breuer et al., 2019). For a given second-order orientation tensor, 
closure approximations (Advani and Tucker III, 1990; Han and Im, 1999; 
Cintra Jr and Tucker IIL, 1995; Chung and Kwon, 2002; Montgomery- 
Smith et al., 2011a) identify a corresponding fourth-order orientation 
tensor based on assumptions. Closure approximations are used by, e.g., 
Jack and Smith (2007; 2008); Buck et al. (2015); Müller and Bóhlke (2016); 
Goldberg et al. (2017) to model mechanical properties. Breuer et al. (2019) 
assess the reconstruction of FODFs based on closure approximations and 
minimum entropy method used by Bóhlke (2005); Müller and Bóhlke 
(2016). This paper addresses the research question on the variety of 
fiber orientation tensors of order two and four. The paper has review 
character and is structured as follows: The definition of FODF and fiber 
orientation tensors is followed by a review on the variety of second-order 
orientation tensors, leading to the orientation triangle. The classical 
parameterization of the orientation triangle is framed and parameteri- 
zations which highlight the isotropic state and material symmetry are 
motivated. Implications of index symmetry of fourth-order tensors 
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in Mandel (Mandel, 1965) notation are discussed, before fourth-order 
orientation tensors are treated as special Hooke tensors and harmonic de- 
composition is applied. This translates the reasoning of Kanatani (1984) 
into continuum mechanics notation, leading to simplified notation with 
focus on the isotropic state. A parameterization of a generic fourth-order 
orientation tensor is given based on a parameterization of the second- 
order orientation tensor in combination with a triclinic structure tensor. 
Admissible orientation states are identified by demanding positive semi- 
definiteness. Results of Nomura et al. (1970); Müller and Bóhlke (2016) 
on admissible transversely isotropic orientation tensors are reproduced 
and framed by coefficient-wise constraints. Minimal sets of discrete fiber 
orientations leading to special transversely isotropic orientation states 
are presented and discussed. The variety of orthotropic fourth-order 
fiber orientation tensors is visualized. Transversely isotropic and planar 
orthotropic fourth-order orientation tensors are highlighted as limiting 
orthotropic cases and admissible parameter ranges are specified. A 
compact parameterization of planar orthotropic fourth-order orientation 
tensors is given. Minimal sets of discrete planar orthotropic fiber orien- 
tations visualize the character of fourth-order orientation information. 
The variety of planar fourth-order fiber orientation tensors including a 
parameterization and admissible ranges is given. This paper closes with 
notes on a small set of closure approximations which are contrasted by 
the variety of fourth-order orientation tensors. 


4.2 Fiber orientation 
Given a reference volume (RV) of arbitrary size and shape, which might 


be interpreted as the part of a body B centered around position x. The 
orientation of axisymmetric fibers inside this reference volume can be 
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described by the fiber orientation distribution function (FODF) 
~:S?>R, with S? = (n € R? |||n|| = 1] (4.1) 


mapping any direction n onto a scalar value (n). S? is the two- 
dimensional surface of a unit sphere parameterized by, e.g., a unit vector 
n. The function v (n) is non-negative, i.e. 


w(n)>0, Ynes? (4.2) 


holds and normalization of v (n) implies 


/ y(n) dn=1. (4.3) 

S2 

As fibers have a direction but no attitude, w (n) is symmetric, i.e. 
v(-n)-v(n), Ynes? (4.4) 


holds (see Górthofer et al. (2020); Advani and Tucker III (1987)). The 
fiber orientation distribution function solely describes the orientation of 
the fibers. Being a one-point statistic information, the FODF contains 
no additional information on, e.g., the spatial arrangement of the fibers 
inside the reference volume. In a heterogeneous material with spatially 
varying microstructure, the FODF usually is influenced by the size of 
the reference volume (see, e.g., size parameter in Górthofer et al. (2019, 
Figure 4) ). If the fibers have identical volumes, e.g., because of identical 
length and constant cross section, the volume fraction cr of those fibers 
pointing through F CS ? on all fibers, is given by 


CF E av= f win) dn (4.5) 


Ve Jvr 
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with the volume of all fibers V; and Vz being the volume of fibers 
pointing through F (see Schemmann et al. (2015); Müller and Böhlke 
(2016)). Equation (4.5) demonstrates the transformation of a volume 
average, e.g., over Vz, into an average over corresponding parts of 
the unit sphere. Volume averages on a reference volume Vay of a field 
quantity q (x) which is axisymmetric about a spatially varying principle 
axis n (x) can be transformed into an average over S? 


1 


Xf aco av - f val) an, (46 
RV JVay S? 


if q can be parameterized solely in n. Similar reasoning for volume 
averages over the special orthogonal group SO(3) in context of crys- 
tal orientation distribution functions (CODF) are extensively used in 
literature, see, e.g., Böhlke (2005); Fernández and Böhlke (2019). ODF 
averages are limited to axisymmetric quantities since FODF is defined 
on directions. In contrast, CODF can be interpreted as rotations of 
coordinate systems and therefore induce no restriction on the quantity 
which is to be averaged. The FODF reflects the material symmetry of the 
microstructure. As a consequence, 


V(S*n)-v(n) YS € Sky C SO(3) (4.7) 


holds. Equation (4.7) implies that w (n) is symmetric with respect to all 
rotations S which are contained in the symmetry group of the microstruc- 
ture Spy being a subset of SO(3). Following Forte and Vianello (1996), 
eight different material symmetries exist and tensor representations of 
Hooke tensors following these symmetries are given in, e.g. Bohlke 
(2001). An inclusion scheme, depicting the hierarchy of the symmetry 
classes combined with illustrations of the symmetry planes, can be found 
in Francois et al. (1998, Fig. 4). Implications of material symmetries on 
second-order tensors are summarized in Appendix A.1. 
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4.3 Orientation tensors of first kind 


In the standard framework of continuum mechanics, physical quantities 
are expressed as invariant tensors. This enables the use of known 
transformation rules and simplifies the storage in computer memory. 
Kanatani (1984) approximates experimentally obtained directional data 
by tensor series and Advani and Tucker III (1987) defines orientation 
tensors of Kanatani first kind (Kanatani, 1984) by 


Nay = y (n) n®* dn (4.8) 
Ss? 


with n®* being the k-th moment of n. For example, the resulting second- 
and fourth-order tensors are 


N = Nv) = i weaned (4.9) 
N = Nyy =| v(n)n@n@n@ndn. (4.10) 
S2 


Equation (4.8) can be interpreted as a weighted summation of moment 
tensors. The moment tensor represents the tensorial character of a spe- 
cific direction and weights are specified by the distribution information 
of the FODF. It follows from Equation (4.8) that N and N are completely 
symmetric with respect to index permutations. Odd orientation tensors 
vanish due to the symmetry of v (n) and higher order tensors contain 
all tensors of lower order as 


Nq-2, =N [I (4.11) 


holds for 2 < k with the identity on second-order tensors I. As a 
consequence of the normalization of both w (n) and n, the limiting case 
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of Equation (4.11) yields 
N-I-tr(N)— 1. (4.12) 


Orientation tensors of Kanatani first kind (Kanatani, 1984) are commonly 
used to represent experimentally obtained directional data, e.g., from 
computer tomography scans or results of flow simulations (Górthofer 
et al., 2019). 


4.4 Orientation tensors of third kind 


The FODF can be expressed as a tensorial Fourier series 


le 
= =>) (x) n®* (4.13) 


T 
k=0 


which is called spherical harmonic expansion (Kanatani, 1984, page 154) 
and introduces orientation tensors of Kanatani third kind (Kanatani, 
1984) 


A 2k +1 (2k 
(k) = gr ( x dev (Nu): (4.14) 


Constructing Dy k) in Equation (4.14) contains two steps, first taking 
the deviatoric part of N(,) and second, scaling with factors for the 
series expansion given by 241 (?") . Skipping the second step, defines 
orientation tensors of third aa as done by, e.g., Müller and Bóhlke 
(2016), as 


(k) — dev (Na) . (4.15) 
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In combination with Equations (4.11) and (4.12) this leads for k € [2,4] 
to 


D, = dev (N) = N— „ltr (N) (4.16) 


=N- N“ (4.17) 


(4) = dev(N)=N- sym ((N [1]) & I) 


+ sym (I @ I) tr (N [I]) (4.18) 
6 3 
=N- zeym (N & I) 4 gg3ym (Ic I) (4.19) 
with the isotropic orientation tensor of second order 


Nise — 1 (4.20) 


The operator dev(-) extracts the deviatoric part and the operator sym(-) 
extracts the totally symmetric part with respect to index symmetry. 
Definitions of the operators sym(-) and dev(-) for higher order tensors 
are given in Spencer (1970). In other communities, alternative repre- 
sentations of directional data are common. The connection to spherical 
harmonics in the context of quantitative texture analysis is discussed in 
Fernandez and Bohlke (2019); Breuer et al. (2019); Kanatani (1984); Jack 
and Smith (2008). 


4.5 Variety of second-order orientation tensors 


Equations (4.2) and (4.8) imply that orientation tensors of Kanatani first 
kind (Kanatani, 1984) and second order N are symmetric and positive 
semi-definite. As a consequence, N can be diagonalized, i.e., pairs 
of eigenvalues à; with A; > 0 and orthonormal eigenvectors v; for 
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i € [1, 2,3] exist, such that 


3 
N= Nf? e, & e; = 5 Av OVi = Ag 0 v;& v; (4.21) 
= sym Ag 


holds and there exists a rotation defined by an orthogonal tensor 
Q =v;®de (4.22) 


mapping the arbitrary but fixed basis {e;} onto the basis {v;}. In the 
following, the orthonormal basis {v;} spanned by the eigenvectors is 
called orientation coordinate system. The arbitrary ordering convention 
of the eigenvalues 


Ag < A2 <Aı (4.23) 


is common in the literature (see Goldberg et al. (2017); Mentges et al. 
(2021)). The constraint in Equation (4.12) is equivalent to 


MA 3-21 (4.24) 


and reduces the number of independent components of N from six 
to five. The coordinate system {v;} or equivalently the mapping Q is 
defined by three scalars, leaving two scalars specifying structural infor- 
mation on fiber orientations inside the orientation coordinate system. As 
a consequence, the variety of N can be expressed as a two dimensional 
space known as the orientation triangle, e.g., parameterized in pairs 
(Ai, Aa), in combination with a mapping Q which defines the orientation 
coordinate system. Information of N is limited to the definition of a 
coordinate system and two half axes of an ellipsoid aligned with theses 
coordinate axes, see Cowin (1985). The third half axis of the ellipsoid is 
defined by Equation (4.24). The weakest material symmetry which can 
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vlrwir 


B 
0 T T T CAEI | | 
g 7 T 2 1 1 2 1 Aı 
3 2 3 3 3 3 
= Ao = planar = transv. iso. V3 
=  d2=)s3 = transv. iso. vı === orthotropic 
—A-1-A41—-232-0 Q isotropic @ unidirectional 
mm (Aa € M) A (Aa € A2) ^ (0 € Aa) Q planar isotropic 


(a) (b) 


Figure 4.1: (a) Constraints defining the orientation triangle (b) Material symmetries in the 
orientation triangle (Cintra Jr and Tucker III, 1995) 


be described by N is orthotropy (see Cowin (1985) or Appendix A.1). 

Representations of orientation triangles are given, e.g., in Cintra Jr and 
Tucker III (1995); Chung and Kwon (2002); Goldberg et al. (2017); Kóbler 
et al. (2018) as well as in Figure 4.1a. An alternative visualization 
is given by the orientation invariant map in Chung and Kwon (2002, 
Figure 1b). The triangle in Figure 4.1a is only one of possible orientation 
triangles and is called standard orientation triangle. The boundaries of 
the triangle are labeled in Figure 4.1a and follow directly from Equations 
(4.23) and (4.24) and 0 € As. Following Cintra Jr and Tucker III (1995), 
vertices and edges of the triangle are discussed, based on Figure 4.1b. 
The vertices A, B and C are given by the extremal orientation states, 
which are isotropic, unidirectional and planar isotropic. Starting from a 
planar isotropic state (Ay = Aa = 1/2, A3 = 0) at vertex C, moving along 
the orange edge towards vertex A, A; increases while A; and A; decrease 
uniformly until the isotropic state (A; = Aa = A3) is reached. Continuing 
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from vertex A along the edge towards vertex B, the largest eigenvalue 
A1 increases while Az and A3 decrease uniformly until the unidirectional 
state (A; = 1, Ag = As = 0) is reached. All states along the edges CA and 
AB have at least two identical eigenvalues and therefore are transversely 
isotropic with principle axes being v3 and v, respectively. The material 
symmetries of points of the orientation triangle are visualized in 4.1b 
and discussions can be found in Cowin (1985); Cintra Jr and Tucker III 
(1995). The states along the edge BC are planar as A3 vanishes. Starting 
from vertex B along the edge towards vertex C, ^, decreases and Aa 
increases until both are equal. All points inside the triangle including 
the edge between vertex B and C are orthotropic. In summary, two edges 
are transversely isotropic, one edge is planar and one point is isotropic. 


The triangles in Figures 4.1a and 4.1b are projections of one sixth of a 
orientation plane in the three-dimensional space spanned by {v;} given 
in Figure 4.2a. For each point in the orientation triangle, five correspond- 
ing points with identical structural properties and different ordering 
conventions of eigenvalues exist. To illustrate the redundancy which is 
inherent in the orientation plane, a randomly chosen orientation state 
with eigenvalues being any permutation of (1/2, 1/3, 1/6) is marked in 
orange multiple times in Figures 4.2a and 4.2b. For each eigenvector 
in {v;}, transversely isotropic orientation states with principle axis v; 
are located along a straight line being the interSection of the orientation 
plane and a plane with Am = An for i 4 m Æ n. Planar states are found 
at the interSection of the orientation plane and planes of one vanishing 
eigenvalue. The central character of the isotropic state is reflected 
by Figure 4.2a. The three-dimensional representation in Figure 4.2a 
extends existing visualizations, e.g., in Goldberg et al. (2017, Figure 2), by 
adding material symmetries and motivates alternative parameterizations 
introduced in the next section. 
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T T T T T 
3 2’ X o 9 1 2o 
3 3 3 3 
(a) (b) 
IN| € 1 mem plane oo triangle = transv. iso. v4 
== transv. iso. V3 === transv. iso. V2 
Q planar iso. vi 


Q isotropic 
@ unidirectional va 


@ planariso.v3 © 


@ unidirectional vı 
© planar iso. va @ unidirectional v3 
perm(1/2, 1/3, 1/6) 


(c) 


Figure 4.2: (a) Orientation plane, triangle and selected orientation states in the space 
spanned by the eigenvalues of N (b) Plane, triangle, selected states and states constrained 


by the norm of N in parameter space (o, a3) (c) Shared legend of Figures 4.2a and 4.2b 


4.5.1 Parameterizations of the orientation triangle 


The classic parameterization of the orientation triangle in eigenvalues A; 
and A» is 


N (Ai, Aa) = ivi ® V4 + rove Q v2 + (1 — Ay — Aa) Va Q Va 


(4.25) 
with 


<A <1 


Wi] = 


and =(1 = i) < Aa < min (Ai, 1— Aı) : (4.26) 
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Introducing two transversely isotropic deviatoric structure tensors 


1 0 0 
pransvl _ -1/2 0 | vi& vj (4.27) 
sym —1/2 
-1/2 0 0 
pne = -1/2 0| vi& vj (4.28) 
sym 1 


enables an alternative parameterization of the orientation triangle by 


N (a, a3) 5 Nise 4. ake! s jb (4.29) 
with 5 i 
Q 
0<a< and Ta ea (4.30) 


The parameterization in Equation (4.29) highlights the central role of the 
isotropic state, which is reached for a; = a3 = 0. Positive values of aı 
lead to a deviation from the isotropic state towards the unidirectional 
state in direction vı. Deviation towards the planar isotropic state with 
principle axis v3 , i.e., away from the unidirectional state in direction v3 is 
described by a3. Figure 4.2b shows the orientation triangle as part of the 
orientation plane in the parameter space (a1, a3). The orientation plane 
itself is part of those pairs of (01,03) which lead to orientations with 
norm less than one. If the separation into two transversely directions 
is not required, F'?^*! and F'?"5? can be combined to a generally 
orthotropic structure tensor of second order, leading to 


N (à, ê) = Ni? + portho (â, ê). (4.31) 
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with 


£g» 
j=) 


possis sje — (â+ 6) 
sym 


Vi & Vj (4.32) 


œ oO o 


A parameterization reflecting the arbitrariness of the ordering conven- 
tion on A; and deploying barycentric coordinates is given as a function 
of three non-negative weights 0 < ù; with i € [1,2,3] and ûs = va Wi 
by 


N (4, i9, 3) = N (^ o need s (4.33) 
wS 


requiring Ws > 0. See Goldberg et al. (2017, Figure 2) for a visualization 
of this parameterization. Each parameterization given in equations 
(4.25), (4.29), (4.31) and (4.33) in combination with a rotation, following 
Equation (4.22), can be used to represent all possible second-order 
orientation tensors. A table with parameters of special points in all 
parameterizations is given in Appendix A.2. 

Parameterizations of all transversely isotropic N, not being restricted to 
the standard orientation triangle, are given by 


NESS (A1) = N (A4, A2 = (1 — A1) /2) (4.34) 
Nee (ay) =N (a1, Q3 = 0) . (4.35) 


A parameterization of transversely isotropic N inside the orientation 
triangle requires a piece-wise definition due to the kink at the isotropic 
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state leading to 


Ntransv (A2) = (4.36) 


(1 — 2X2) v1 Q v1 + À2V2 Q V2 + AQV3 ® V3 if0 € As «i 
1 
2 


A2V1 G V1 + À2V2 Q V2 + (1 — 203) va @ v3 ificAXz 


Planar second-order orientation tensors inside the orientation triangle 
can be parameterized by, e.g. 


1 
planar (a1) =N c a3 = ar _ 3) (4.37) 


with 0 € o4 < 2/3. 


4.6 Variety of fourth-order orientation tensors 


Equation (4.10) implies that N is completely (or totally) symmetric, i.e., 
any permutation of the indices 

NS = NS un (4.38) 
holds for a representation with tensor coefficients N — NS je; ej ® 
ep © ei. Complete (index) symmetry reduces the number of independent 
components of a generic fourth-order tensor in three dimensions from 
81 = 31 to 15, as there are 15 = ("^7 !) unordered ways of choosing a 
combination of k = 4 elements from a set of n = 3 elements. Mandel 
notation, introduced in Mandel (1965) and also known as normalized 
Voigt notation (Mehrabadi and Cowin, 1990; Bóhlke, 2001) enables 
compact two-dimensional representations of fourth-order tensors with at 
least minor symmetry. A fourth-order tensor A = Aij; €i © ej & ex 8 ei 
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is minor symmetric if it has both minor symmetries, i.e. Ajjxi = Ajikt = 


A;jır holds. Introducing base tensors in an arbitrary Cartesian basis {e;} 


by 
Bı =e, 8e1, 
By = ez Q es, 
B3 = ea 8 es, 


Bı = 


B; = 


Bg = 


(eo & es - 


(ez & ea 4 


SEELE 


- e3 @ ea), 


Fea Q ei), 


(eo C ei = 


-e1 ea), 


(4.39) 


a minor symmetric tensor A is represented by a six by six matrix of 


coefficients Ag; 


A = Ajjri €i Q €j Q ex Ber = Ac; Be @ Be = 


A1111 
A2211 
A3311 
V2A2311 
V2Ai311 
V2A121 


Be 8 Be 


with € and ¢ summing from 1 to 6. A Hooke tensor 


A1122 
A2222 


A3322 
V2A2322 


V2A1322 
V2A1222 


A1133 
A2233 


A3333 
V2A2333 


V2 A1333 
V 241233 


V2A1123 
V2.42223 
V243353 
242323 
241323 
241223 


V2A11s 
V242213 
V 23313 
242313 
241313 
241213 


V2A1112 
V242212 
V2A3312 
242312 
241312 
241212 


(4.40) 


is minor and major 


symmetric, i.e., in addition to the minor symmetry condition it holds 
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that Bijkı = Bruj. The coefficient matrix of a Hooke tensor 


V2Bi4 
V2Bis 
V2Bi6 


Be 8 Be 


V2B»4 
V2B»s 
V2B»s 


By V2Bi4 V2Bis 
Bog V2Bo,  V2Bss 
B33 V2Bs4  V2B35 


V2Ba 2Bas 
V2Bss  2Bas 
V2Bsg 2B4s 


2Das 
2 Bss 
2 Dsg 


V2Bi¢ 
V2 Bag 
V2B36 
2D46 
2Bs6 
2 B66 


(4.41) 


is symmetric and contains 21 independent parameters. Complete index 


symmetry of N implies 


Niv ONE NG v ANI) UNIO VINIS 
NY NG VNG VING VANS 
NGANG) Vans VIND 

on) 2NG 2NG 

major symmetric on = an 
2NO 


B; 8 Be. 


(4.42) 


In Equation (4.42) indices of redundant tensor coefficients are colored. 


The redundancy implies that six coefficients in the upper left quadrant 


and nine coefficients in the upper right quadrant of the coefficients in 


Mandel notation define a completely symmetric tensor. This motivates a 


short hand notation „completely symmetric", which to the best of the 


authors' knowledge has not been used in the literature so far, see, e.g., 
Equation (4.55). As N contains N, the constraint on the trace of N in 


Equation (4.12) reduces the number of independent components of N by 


one to 14. In the literature erroneous implications of Equation (4.12) on 


the number of independent components of N are found, see, e.g., Dray 
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et al. (2007). Expressing N by contraction of N, i.e., 


4 4 4 4 4 4 4 4 4 
Nie uic NÉS T Ne ay NE ae Nee NE T Noe Tv Na | 


4) 


N 
4 4 4 4 4 4 4 4 4 
| NT T NG a D ADT er Noe By Nye Wie zu o au i ] 
A 4 4 4 4 
| N ER N. "WS NS NS m Ni EE NY Ni Me NS x NS 


reveals the implication of Equation (4.12) on N to be N. d +N. (9 +N. (9 + 
2NG +2NG) e 3NS = 1. 


4.6.1 Harmonic decomposition 


Linear invariant decompositions, including the classic harmonic decom- 
position, are frequently used on Hooke tensors. The fourth-order orien- 
tation tensor N is a Hooke tensor and has additional index symmetry. 
Therefore, linear invariant decompositions can be applied to N to study 
its structure. Literature on invariant decompositions of fourth-order 
tensors is extensive, see, e.g. Backus (1970); Mochizuki (1988); Baerheim 
(1993); Boehler et al. (1994); Forte and Vianello (1996); Olive et al. (2017) 
or those focusing on Hooke tensors and addressing an engineering 
audience Rychlewski (2000; 2001); Böhlke (2001). 

Following Rychlewski (2000), a Hooke tensor H with 21 independent 
components can be split into five parts (K, G, H1, H5, dev (H)) leading 
to 


H = 3KIP + 2GP» + sym, (I ®H,+H, Q T) (4.44) 
+ sym, (I & H5 + H5 8 I) + dev (H) 


with isotropic projectors P; and Po, two symmetrization operators 
sym, (-), sym, (-) and the deviatoric operator dev(-). The symmetrization 
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operators have to follow a special structure which is described in 
Rychlewski (2000). Two common choices of the symmetrization 
operators are used in the following. The numbers of independent 
parameters of the five parts (K, G, Hi, H2, dev (H)) are (1, 1, 5, 5, 9) 
with the scalar compression modulus K and shear modulus G. The 
isotropic parts are obtained by projections, i.e. 


1 


- I Bu=!1P-H. (4.45) 
IP2| 


3K =P, -H, 2G : 


The parts H; and Hz depend on the choice of the operators sym, (-) and 
sym, (-) but for any choice are only functions of the dilatational modulus 
C and the Voigt tensor V which are given in index notation by 


Ci; = Hijkk and Vij = Hikjk (4.46) 


with C = C;i;ei 9e;, V = Vije; & e; and H = Hijy ei $e; 8 ex 8 er. 
The classic harmonic decomposition, e.g., used in Bóhlke (2001), is 
obtained for sym, (-) being the identity and sym, (-) symmetrizing the 
only remaining index asymmetry of Hooke tensors, leading to 


Hæhssic — L (5dev (C) — 4 dev (V)) (4.47) 


classic __ 
H5 = 


(3dev (C) — 2 dev (V)). (4.48) 


NR Ne 


Applying the classic harmonic decomposition to N reveals that the 
dilatational modulus and the Voigt modulus of N coincide and both 
are N. Following Cowin (1989), a Hooke tensor whose dilatational 
modulus and Voigt tensor coincide, is said to fulfill the Cauchy relations. 
Any completely symmetric tensor and thus also any orientation tensor 
fulfills the Cauchy relations. Referencing Musgrave (1970), Cowin (1989) 
states about the Cauchy relations ,, [...] the relations do not hold for most 
elastic materials, but only for materials which can be described as having 
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u 


central-force laws operating between points of a simple lattice, [...] 
This observation fits the model of fibers pointing towards the origin of a 
unit sphere. A stiffness which is linear in the orientation tensors does not 
have to fulfill the Cauchy relations as its deviatoric parts (Hı and H3 in 
Equation (4.44)) might be multiples of N and therefore do not coincide. 
Motivated by their orthogonality, Rychlewski (2000) proposes an alterna- 
tive pair of symmetrization operators sym, (-) = sym(-) and sym,(-) = 
I° — sym, (-) focusing on the index symmetry with the identity on sym- 
metric fourth-order tensors I°. This choice leads to the corresponding 
parts 


Erde _ : (dev (C) + 2 dev (V)) (4.49) 
Hirdex — dev (C) — dev (V). (4.50) 


As a consequence, N can be decomposed into 


" X 3 3 
N=3KP,+2GP.+ sym (x & „dev (N) + dev (N) ® 1) 


+ dev (N) (4.51) 


= 3K P, + 2G P5 +2 sym (Faew (N) ® 1) + dev (N) (4.52) 


— 3K P, - 28 P5 + : sym (dev (N) & I) + dev (N). (4.53) 


Due to the isotropy of the isotropic projectors, the linearity of the scalar 
product and the normalization of 4 (n), the projection of N on either of 
the isotropic projectors is 


P; -N= P; | v (n) n® dn = P; . n94 with i € [1,2]. (4.54) 
S2 


Consequently, the isotropic part of N can be determined analyzing the 
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fourth-order moment n?* which is given in Mandel notation by 


ni hing ning Vantnanz Van?nz V2nns 
ng n2n2 V2n3ng V2ninjns V2n4n3 
n2t ni V2nan3 V2nin3 V2nin2n3 
completely symmetric 
Be B; (4.55) 
leading to 
11 1000 
11 1000 
R 111 1 100 0 
Ska P; n8t = B: & B, - (n84 
i 3/0 0 0 0 o o| 2B) 
000 0 0 0 
00 0 0 0 0 
1 
= 3 (ni + Zning + Znin; + ng + 2ngn3 + n3) 
1 
=L 4.56 
; (4.56) 
and 
2 -1 -1 0 0 0 
-1 2 -1 0 0 0 
^ P» 1 |-1 -1 2 000 
2G = ine ec B; & B, - (n®? 
(Pall? Blo 0 0 3 0 of BeBe) 
0 0 0 03 0 
0 0 0 00 3 
2 
= —. 4.57 
15 (4.57) 


65 


4 Variety of Fiber Orientation Tensors 


Thus, the isotropic part of any fourth-order fiber orientation tensor is 
fixed due to the normalization of both v (n) and n and is given by 


311000 
131000 
l . i 11113000 
N59 = 3KP, + 2GP, = — B:®B 
: = dB Ig. 0020 0o 57s 
000020 
000002 
7 
= gym (Lo 1) (4.58) 


leading to a compact representation of Equation (4.53) by 
N=N®+ "sym (dev (N) @ I) + dev (N). (4.59) 


This representation emphasizes isotropy as the one element of directional 
measures. Equation (4.59) reveals that in contrast to the well known har- 
monic decomposition of Hooke tensors, the harmonic decomposition of 
fourth-order orientation tensors only contains one irreducible subspace 
of second order and the isotropic subspace degenerates from two scalars 
to a constant value. Coincidence of Equations (4.59) and (4.19) is given 
by 

sym (I ® dev (N)) = sym (I 8 N) — ssym (IGI). (4.60) 


Interpretation of N as a Hooke tensor and applying the harmonic decom- 
position translates the reasoning of Kanatani (1984) into the language 
of continuum mechanics. N can be parameterized by two deviators, 
one being of second and the other of fourth order. Both represent the 
deviation from the isotropic state. Five of the 14 independent parameters 
of N define the second-order deviator and nine define the fourth-order 
deviator. Three independent parameters of the second-order deviator 
define an orientation coordinate system in which both deviators have 
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simplified representations (see the structure tensor in Equation (4.31) 
and the following section). 


4.6.2 Parameterizations and admissible parameter ranges 


The space of fourth-order orientation tensors N, which fulfill the alge- 
braic constraints imposed by Equation (4.8) and which can be derived 
from an FODF, remains to be identified. The variety of N is known from 
Section 4.5. Therefore, Equation (4.59) reveals that the analysis of the 
variety of N is reduced to the identification of a fourth-order deviator 
dev(N). This fourth-order deviator reflects the material symmetry of 
the FODF and consequently also of N and N. As a consequence, special 
representations of the fourth-order deviator in the orientation coordinate 
system exist. Knowledge on material symmetries of Hooke tensors is 
applied to identify admissible ranges of the fourth-order deviator. The 
most general fourth-order orientation tensor has the structure 


: 6 ; 
N (N, dı, do) = N? + = sym (dev (N) @ I) + Fl (di... dg) (4.61) 


with a triclinic deviatoric structure tensor 


FFtrid (di, ..., dg) = (4.62) 
—(di + d2) di da —vV2 (d4 + d5) V2de V2 dg 
—(dı + d3) d3 V2 d4 — V2 (de + dz) V2 do 
— (da + da) | V2 ds V2 d; — V2 (dg + do) 
completely symmetric 
¿8 BY 


which is a function of nine scalar parameters. The Mandel basis Br 8 BY 
is spanned in the orientation coordinate system, i.e., for example the 
first definition in Equation (4.39) becomes B, = vı & vı. Without loss of 
generality, any orthotropic parameterization of N, e.g., Equation (4.35), 
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Material symmetry reflected by 

Symmetry class of || N N 

isotropic isotropic isotropic 

cubic cubic isotropic 

transv. iso. transv. iso. || transv. iso. 
tetragonal tetragonal | transv. iso. 
trigonal trigonal transv. iso. 
orthotropic orthotropic || orthotropic 
monoclinic monoclinic || orthotropic 
triclinic triclinic orthotropic 


Table 4.1: Implications of symmetry classes of the FODF w on the second- and fourth-order 
fiber orientation tensors 


can be combined with Equation (4.62) yielding a set of parameters, 
e.g., (ai, a3, di, ..., dg). If the FODF has a material symmetry, this 
symmetry also applies to the orientation tensors and reduces the number 
of independent parameters, see, e.g., the implications of transverse 
isotropy on the parameterization of N in Equation (4.35). For a given 
material symmetry class, the number of independent parameters is fixed. 
All eight material symmetries imply different constraints on N. However, 
e.g., triclinic, monoclinic and orthotropic symmetry of 4 lead to identical 
constraints on N. In other words, triclinic, monoclinic and orthotropic 
symmetric second-order tensors coincide and are usually referred to as 
orthotropic (see Appendix A.1). Implications of the symmetry of the 
FODF on the symmetry of N and N are listed in Table 4.1. The hierarchy 
of material symmetries including visualizations of symmetry planes is 
given in Francois et al. (1998, figure. 4). As a consequence, index and 
material symmetry constrain the space of N. The question on the variety 
of fiber orientation tensors is reduced to the determination of admissible 
parameter combinations. The flexible parameterization of fourth-order 
fiber orientation tensors in Equation (4.61) is new, highlights the isotropic 
state and can be easily adapted to stronger material symmetries. The re- 
maining algebraic constraint on N is positive semi-definiteness. Positive 
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semi-definiteness implies that 
0<N-(A@A) (4.63) 


holds for any second-order tensor A. Without loss of generality, A can be 
restricted to symmetric second-order tensors as N is symmetric and skew, 
and symmetric parts are orthogonal. As a consequence, Equation (4.63) 
can be written in Mandel notation stating a root finding problem of the 
characteristic polynomial of N of degree six. N is positive semi-definite 
if all of its eigenvalues A; are non-negative, i.e. 


0 € A; fori € [1, ..., 6) (4.64) 


holds. Explicit representations of the eigenvalues of N are given for a 
limited number of material symmetries in Mehrabadi and Cowin (1990). 
Eigenvalues of fourth-order tensors can not be identified from coefficient 
matrices in Voigt notation. Mandel notation has to be used (Mehrabadi 
and Cowin, 1990; Bóhlke, 2001; Mandel, 1965). 


4.6.3 Transversely isotropic case 


A transversely isotropic harmonic, i.e., completely symmetric and com- 
pletely traceless, tensor of fourth order with principle axis in direction 
ei; is a multiple of the structure tensor 


8 —4 -4/0 0 0 
3 1 0.0 0 
0 
rtransvi = 3 0 0 B. E B, (4.65) 


completely | symmetric 


B Fed (4, —4,1,0,0,0,0,0, 0), 
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see, e.g., Fernández and Bóhlke (2019, (75)), Lobos Fernández (2018, 
A.3) or Müller and Bóhlke (2016, (11)) and can be expressed in terms of 
the triclinic structure tensor defined in Equation (4.62). In consequence, 
transversely isotropic fourth-order orientation tensors are parameterized 
in the orientation coordinate system by 


New (a, p) = NP? +a i sym [pios e I) + p Füansvi — (4.66) 
6a/7+8p a/l4—4p a/l4—4p |0 0 0 
—3a/7+3p -a/7+p |0 0 0 
Nio 4 -3o/T--3p|0 0 0 YoBy. 
completely symmetric 


The choice of the principle axis only effects the mapping Q, defined in 
Equation (4.22), and parameterizations around v2 or v3 are equivalent, 
see, e.g., Müller and Böhlke (2016). Demanding positive eigenvalues 
leads to the set of admissible transversely isotropic fourth-order orienta- 
tion tensors 


2 
N ransv = {Ne (a, p) | —--~<a< 3” (4.67) 


13 T ee 
~ "99^ ^ 56% 60 


given, e.g., in Nomura et al. (1970); Müller and Böhlke (2016). The set 
Ns combined with the parameterization in Equation (4.66) represents 
the variety of transversely isotropic fourth-order orientation tensors 
(Müller and Bóhlke, 2016). 

An alternative derivation of A/**'*Y is given by Nomura et al. (1970, 
Equation (54)) deploying the Schwarz inequality on directional cosines 
to identify a relationship between N and N. This relationship implies 
constraints on the coefficients of N which are reformulated in the notation 
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of this work as 
(NQ vi & vi) "«NÜBreBrzNÜv,ev, (4.68) 
for k € [1,2,3] without summation convention. To be precise, the set 
Nees, = [wem (4.69) 
(NQ vy ® vi) < NY BY @ BY «NO v.a ve} 


is identical to the set "SY. However, for weaker symmetries, positive 
definiteness is a stronger constraint than the conditions in Equation 
(4.68). Comparison of second- and fourth-order moment tensors (see, 
e.g., Equation (4.55)) motivates the reasoning behind Equation (4.68). 
In order to frame the implications of positive semi-definiteness, Figure 
4.3 contains additional subsets, lines and special points discussed in 
the following. Interpretation of N as weighted summation of moment 
tensors leads to coefficient-wise constraints on N (see Appendix A.3). 
The set AE?"*", defined in Equation (A.15), fulfills these coefficient-wise 
constraints. The set of transversely isotropic orientation states bounded 
by the norm of a singular orientation is denoted by 


AD = {Nias | 0 < IIN] Frobenius $ 1} z (4.70) 


Figure 4.3 visualizes the relation Nav c Mias cC Nav, Spe- 
cial orientation states, such as the planar isotropic and unidirectional 
state are marked in Figure 4.3 for reference. The isotropic state N*9, 
defined in Equation (4.58), is located at a = p = 0 and is part of a 
straight vertical line with a = 0 which contains orientation states with 
isotropic second-order tensor. Therefore, the set of transversely isotropic 
orientation states with isotropic second-order orientation tensor is a 


one dimensional subspace of N "A"sv with limiting orientation states 
P 8 
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p 
ee a pee Aftransv 
9» fi] transv 
linear 
1 transv 
60 E Nason 
3/280 4 — ND=0 
= NO 
NO =0 
"" — NO =0 
— (4) — 
i N@=1 
— (4) — 
= N» = 1/4 
=== [BOF closure 
=== ORF closure 
TA === linear closure 
= @ isotropic 
@ unidirec. vı 
3 ; 
-5 1 _ X Q plan. iso. vi 
— e ^ 


N 
m 
© 


1a  @B 
6 


Figure 4.3: Transversely isotropic fourth-order orientation states following Nomura et al. 
(1970) 


visualized in Figure 4.3 and given by 


A= Atransv (a = 0, p = 1/60) 


EE 0 0 000 | 
1/4 1/2|]0 0 0 
1/4 |0 00 


= Ye BY (4.71) 


completely symmetric 
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and 


= Ayfransv (a = 0,p = —1/90) 


1/9 1/9 1/9 |0 0 0 
1/6 1/18 0 0 0 
1/6 10 0 0 x " 
= £ @ Be. (4.72) 
completely symmetric 


The states A and B can not be distinguished by means of second-order 
orientation tensors as 


A [I] = B [I] = N“ (4.73) 


holds, but the states differ significantly. Figure 4.3 shows that for the spe- 
cial case of transversely isotropy, only two second-order fiber orientation 
tensors, namely those defined by a = —1/3 and a = 2/3, are connected 
to a corresponding fourth-order fiber orientation tensor by a one-to-one 
relation. The subsequent sections show that for weaker symmetries a 
one-to-one correspondence solely remains for the unidirectional case. 
Closure approximations, shortly discussed in Section 4.7, construct a 
one-to-one mapping between a given second-order fiber orientation 
tensor and an unknown fourth-order fiber orientation tensor. In Figure 
4.3, three closures which are discussed in Section 4.7 are visualized 
and indicate that the variety of transversely isotropic fourth-order fiber 
orientation tensors is not fully reflected by the selected closures. 


An intuitive interpretation of the orientation states described by A and 


in Equations (4.71) and (4.72) is obtained by minimal sets of discrete 


fibers which are described by these fiber orientation tensors. The single 
fibers can be visualized and represent minimal realisations of states in 
N™rsv_ To the best of the authors’ knowledge, visualization of fiber 


orientation tensors by minimal discrete fiber arrangements is new. For 
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both the state A and the state B, an FODF can be build by a sum of delta 
functions with homogeneous weights 


1 1 
y (n) = N2 5 6 n ni) +ô (n + nj) (4.74) 
based on a set of discrete fiber directions (n; with i € [1, ..., N]}. Mini- 


mal sets of discrete fiber directions leading to orientation tensors À and 


are 
2r(i— 1 
Da = 4 |cos (6) | with 6; = in ) tori € ll...) (4.75) 
sin (¢;) 
1 
+34 10 and 
0 
cos (7/4) boe of 
Dg = cos (¢;) | with ¢; = emo D fori € [1, ..., 6] (4.76) 
sin (¢;) 


given in the orientation coordinate system {v;}. Visualizations of the 
sets of directions are given in Figures 4.4b and 4.4c. Comparison of 
the vector sets and inspection of Figures 4.4b, 4.4c and 4.3 allow for an 
interpretation of the parameter p in Equation (4.66). The lower the value 
of p, the more pronounced are directions towards diagonals between the 
principle axis of the orientation coordinate system. 


4.6.4 Orthotropic case 
An orthotropic harmonic tensor of fourth order is given by 


poto (di, da, d3) = Be (di, d», d3, 0, 0, 0, 0, 0, 0) (4.77) 
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V3 


v2 Vi 


(a) (b) (c) 


Figure 4.4: Visualization of minimal sets of discrete directions representing extremal 
transversely isotropic orientation states which are N-isotropic: (a), (b) Two views on Da 
leading to A, (b) Dg leading to B 


with Fe! defined in Equation (4.62). As the weakest material symmetry 
of second-order tensors is orthotropy, a parameterization of orthotropic 
fourth-order orientation tensors is given by 


Ne (N, di, da, da) = 


i 6 
N59 + z sym (dev (N) & I) + F&° (dı, da, ds). (4.78) 
The set of admissible orthotropic fourth-order orientation tensors is 


Northo = m (N, di, də, ds) | (4.79) 


0 < A; V A; in eigenvalues of M. 


see Section 4.6.2. 
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Parameterization of N following Equation (4.29) leads to 


Northo (an, a3, dı, do, d3) — Nie. (4.80) 
3 (201 - a3) - dı - d2 dj (a1 - 203) + di i4 (a1 +43) + da 000 
-3 (a1 +3) -dı -d3 d; (-201+a03)+d3 |000 
3 (-a1 + 2a3) - d2 - d3|0 0 0 


completely symmetric 


e£ G Be. 


Three eigenvalues of N?'t^? in the six dimensional space of the Mandel 
representation are directly given by 


1 1 1 

jM = J5 (—204 + a3) + da + ig (4.81) 
1 1 1 

gfe =i (ay + a3) + d2 + Ip’ (4.82) 
1 1 1 

EN = aa 4. 
5s 14 (ay 203) + di + 15° ( 83) 


due to the diagonal form of the lower right quadrant of N" in Equation 
(4.80). Following Mehrabadi and Cowin (1990, Equation (5.35)), the re- 
maining eigenvalues (A4, A5, A6) are the roots of the characteristic poly- 
nomial 


det(M — QI) =0 (4.84) 


being cubic in Q, with the identity matrix of dimension three I and the 


matrix 


M= (4.85) 


3 (2a1 - a3) - dı -d2 + $ Zum - 2a3)+ di +% — d (oi o3) - do + 
-3 (a1 +a3) -dı -d3 + $ d; (-2a1 + a3) + da + $ 
symmetric 3 (-a1 + 2a3) - da - da +4 


w] 
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For transversely isotropic orientation tensors, a compact and explicit 
representation of the admissible parameters, is given by Equation (4.67). 
For the orthotropic case, explicit formulas specifying M% by combina- 
tions of the parameters of N°"!"°, can be identified analytically, e.g., using 
complex or trigonometric expressions. As an alternative, non-negativity 
of eigenvalues can be demanded by alternating signs of monomials of 
the cubic characteristic following Vieta’s formula preventing the need to 
explicitly calculate the eigenvalues. However, the resulting expressions 
are lengthy and therefore are not given here. For applications N "!"° can 
be calculated numerically. 

However, for a given specific second-order orientation tensor N, com- 
pact representations of sets of admissible tuples (di ‚da , d3), i.e., tuples 
leading to No (N, di, do, d3) € N°™°, exist. For reference, an ex- 
plicit parameterization of admissible tuples for the special case of an 
isotropic second-order orientation tensor N'? is given in Appendix A.4. 
Admissible tuples for specific N along three paths are visualized in 
Figures 4.6, 4.7 and 4.8. The paths are defined in Figure 4.5a. Figure 
4.6 reveals that the admissible region of (di, da, d3) for an isotropic 
second-order orientation tensor is symmetric and comparably large. 
Departing from the isotropic second-order orientation tensor along the 
blue path along v, (see Figure 4.5a) towards the unidirectional state, 
the admissible region of (di, d2, d3) changes its shape and position in 
the space spanned by dı, d» and d3 and finally degenerates to a single 
point representing the unidirectional state. As all points on the blue 
path are transversely isotropic with respect to the v;-axis, the admissible 
ranges of the parameters d; and d» are identical along this path (see 
Ftnsvl in Equation (4.27)). Figure 4.7 contains views on the admissible 
tuples for specific second-order orientation tensors along the orange 
path, in 4.5a. Second-order orientation tensors along the orange path 
between the isotropic and the planar isotropic state are transversely 
isotropic with principle axis v3. Therefore, admissible ranges of the 
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dy 
a3 19/140 ~ 
m 
ziad 34/315 4 
—1/6 4 41/560 4 
1/4 4 
use] 3/280 -| 
T T T T T 04 
0 1/6 1/3 1/2 2/3 Ay 
-€ path along vı -€ path along v3 
-€ mixed path —— planar states 
Q isotropic @ unidirectional vı 
@ planar iso. v3 mem A Planar ortho un 
4 1.rowfigure9 @ 2. row figure 9 
0 2/9 1/3 1/2 2/3 1 
(a) 
(b) 


Figure 4.5: (a) Paths on the orientation triangle defining specific second-order orientation 
tensors used in Figures 4.6, 4.7 and 4.8 (b) Admissible ranges of dı for planar orthotropic 
orientation states and specific orientation states used in Figure 4.9. The legend is shared 
with Figure 4.5a. 


parameters dz and d3 change homogeneously and shrink along the path 
towards the planar isotropic second-order orientation tensor. For the 
planar isotropic second-order orientation tensor, the parameters d; and 
dą are fixed and equal —3/70. Therefore, at this second-order state the 
admissible region degenerates to a line and the remaining degree of 
freedom of an orthotropic fourth-order orientation tensor with planar 
isotropic Nplanar iso is d} with admissible range —4/35 < d, < 19/140. 


E a=0,a3=0 

E 4=1/6,03=0 
E @1=1/3,a3=0 
E a1=1/2,03=0 
E 0,219/30,25-0 


0.1- 


d200- 


-01+ 0.1} 


La aaa aaa aa O 1 L eo re " 
-0.1 0.0 0.1 -0.1 0.0 0.1 -0.1 0.0 0.1 
di di d2 


4 Ada aaa daa 


Figure 4.6: Three views on admissible ranges of (di, d2, d3) for specific second-order 
orientation tensors along the blue path (along v1) in Figure 4.5a 
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E =0,03=0 
f E a,20,a52-1/12 
al E a4=0,a5=-1/6 
E a1=0,03=-1/4 


1 a,=0,a3=-19/60 


935.0- 


0.14 


| PA RE BC RE BE BE BR ENS 
-0.1 0.0 0.1 
d2 


Figure 4.7: Three views on admissible ranges of (di, da, d3) for specific second-order 
orientation tensors along the orange path (along v3) in Figure 4.5a 


E 0,-0,a5-0 

E a4=1/6,04=-1/12 

© aı=1/4,a3=-1/8 
a4=19/60,03=-19/120 


oit 


-04r -04r 


L ae er ia n i Lag a aaa aaa ia 
-0.1 0.0 0.1 -0.1 0.0 0.1 -0.1 0.0 0.1 
di di d2 


Figure 4.8: Three views on admissible ranges of (di, d2, d3) for specific second-order 
orientation tensors along the black path on orthotropic states in Figure 4.5a 


The planar isotropic orientation state is located at dı = 3/280. Second- 
order orientation tensors of planar orientation states are located at the 
purple line in Figure 4.5a. This line connects the planar isotropic and 
the unidirectional states. Planar orthotropic fourth-order orientation 
tensors are degenerated in the sense that the parameters d» and d are 
directly related to N, ensuring that the orientation state is planar. As 
a consequence, the remaining degree of freedom of planar orthotropic 
fourth-order orientation tensors is solely determined by the remaining 
parameter dı. In consequence, a parameterization of planar orthotropic 
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fourth-order orientation tensors is given by 
Nplanar ortho (a1, di) = Northo (v = planar (o1) ! di, (4.86) 


d2 


_ 15-6 , 15m -6 
- — 3M °° 14 


with NP!" (64) following Equation (4.37). The set of admissible planar 
orthotropic orientation tensors 


N planar ortho __ {eae ortho (o1 ! di) | (4.87) 


2 4 19 9 
DS <a dp ea 
Hx en gt) 


is visualized in Figure 4.5b. For reference, the planar isotropic and unidi- 
rectional states are highlighted in Figure 4.5b. To the best of the authors' 
knowledge, the set of admissible orthotropic fiber orientation tensors of 
fourth order has not been presented in the literature before. Second-order 
orientation tensors which describe planar orientation states in a three- 
dimensional space have four variable parameters. Three of these four 
parameters define the mapping Q (see Equation (4.22)) which defines 
the orientation coordinate system and the fourth parameter defines the 
orientation state. One possible choice for this fourth parameter is o, (see 
Equation (4.37)). Orthotropic fourth-order orientation tensors, which 
represent planar orientation states, add one additional independent 
parameter, e.g., dı. Figure 4.5b visualizes the admissible ranges of dı 
as a function of a1, assuming orthotropy. The admissible range of this 
fourth-order orientation tensor parameter d; degenerates to a single 
value in case of the unidirectional state. This degeneration reveals a 
valuable view on the information which is contained in N but not in N. 
The second-order orientation tensor N contains directional measures in 
the directions defined by the axes of the orientation coordinate system. 
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(b) (d) (£) (8) 


Figure 4.9: Visualizations of minimal discrete sets of directions in the orientation 
coordinate system {v;} for given planar orthotropic orientation tensors which are 
marked in Figure 4.9 and parameterized by a; and dı (see Equation (4.86)). (a) 
a, =0, dı=19/140 with D@ = {[1/V2,1/V2 — 8,0], [1/V2, —(1/V2 — £), 0]} 
with «<1 (b) o1 =0, dı = —4/35 with D®) = {[1,0,0], [0,1,0]} (c) o1 = 2/9, 
dı = 34/315 with DO = ((/2//3,1/V3,0], [V2/V3, -1/V3,0]) (d) a1 = 2/9, 
dı = —4/35 with D@ = {2 x [1,0,0], [0,1,0]} (e) aı=1/3, dı = 41/560 
with D(9 = ((/3/2,1/2,0],[V3/2, —1/2,0]) (f) aı=1/3, dı=-4/35 with 
DM = (3 x [1,0,0], [0, 1, 0]} (g) a1 = 2/3, dı = —4/35 with D(9 = {[1, 0, 0]} 


An orthotropic fourth-order orientation tensor contains additional direc- 
tional information representing the spread around the directions defined 
by the axes of the orientation coordinate system. For the unidirectional 
state, i.e., if all fibers are aligned in one direction, there is no spread. In 
the planar orthotropic case, there is only one spread. This mental model 
is a simplification of directional measures by spherical harmonics but 
suits engineering applications and can be used to interpret Figures 4.6 to 
4.8. In the following, examples of planar orthotropic orientation states 
are discussed. Figures 4.9a to 4.9g visualize sets of discrete directions 
for planar orthotropic orientation states marked in figure 4.5b. The 
upper (lower) row of sub Figures in Figure 4.9 represents discrete sets of 
directions which lead to planar and orthotropic fourth-order orientation 
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tensors with maximum (minimum) admissible values d; (see Figure 
4.5b for reference). A minimal set of discrete directions for each sub 
Figure in Figure 4.9 is given in the caption. The discrete sets of directions 
labeled D® for i € |a, b, c, d, e, f, g] in the caption of Figure 4.9 can 
be combined with Equation (4.74) to identify the corresponding FODF 
and orientation tensors. For reference, limiting cases of orthotropic 
fourth order orientation tensors which contract to the planar isotropic 
second-order orientation tensor are given by 


NPlanar ortho (ay =0,d, = 19/140) = (4.88) 


1/4 1/4 0/0 0 0 
1/4 0/0 0 0 


0/0 0 0 v v 
£8 B? 
completely | symmetric 
and 
NPlanar ortho (o4 =(,dı = —4/35) = (4.89) 

1/2 0 o]0 0 0 

1/2 0/0 0 0 
0/0 0 0 v v 


completely | symmetric 


4.6.5 Planar case 


Without loss of generality, inspection of planar states is restricted to 
orientation states which are located inside the plane spanned by v; and 
v2. Following this convention, some coefficients of moment tensors 
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represented in the orientation coordinate system (see, e.g., equation 
(4.55) or (A.12)) vanish, because out of plane components of the unit 
vector, which builds the moment tensors, equal zero. Combining the 
knowledge of vanishing coefficients of NP!*™** with Equation (4.61) leads 
to the parameterization 


wplanar (o4. dy, dg) = (4.90) 
djo1—dit 25 dict x 0/0 0 V2ds 
-3ai-di- 25 0/0 0 -V2ds 
0000 YoBY 
completely symmetric 


which fits into the general triclinic framework given by equations (4.61) 
and (4.62) with 


Nplanar (o4. di, dg) = N (tans (01) , d, (4.91) 
=f = = 
də = 501 a ds = 15a, a 
140 140 


0, 0, 0, 0, dg, do = Z3 


Similar to the planar orthotropic case, the parameters dz and d3 are 
directly related to N, ensuring that out of plane coefficients in the 
upper left quadrant in Equation (4.90) vanish. Removing the orthotropic 
constraint adds degrees of freedom. However, the planar plane intrin- 
sically acts as a symmetry plane and planar orientation states have at 
least monoclinic material symmetry. As a consequence of this material 
symmetry as well as due to the planarity itself, the parameters d4, 
ds, dg and d; in Equation (4.62), which is utilized in Equation (4.91), 
equal zero. Due to the planarity, dy = —dg holds in Equation (4.91) as 
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the coefficient v2 (dg + dg) in Equation (4.62) has to vanish. This can 
be seen by comparison of Equation (4.62) and the moment tensor in 
Equation (A.12) where the planarity implies 0 = 7. Demanding positive 
eigenvalues of NPlara in Equation (4.90) leads to the set of admissible 
planar fourth-order orientation tensors 


2 
Aims — {NPI (oy, di, de) |O< o1 < 5, (4.92) 
4 19 9 4 
a x qa TU uec. _ , d1) < ds € d 
z $4 Sio En f (a1, di) € da < f(a, )] 
with 
f (oa, di) = (4.93) 
1 


Em 304 — 1260 a? + 420 d, — 11025 a? dı — 19600 d? . 


To the best of the author's knowledge, this set is new in the literature 
and of great practical importance for shell-like components reinforced 
with long fibers. The set of admissible planar fourth-order orientation 
tensors is visualized in Figures 4.10a and 4.10b. The connection to the 
admissible ranges in the planar orthotropic case, which is illustrated in 
Figure 4.5b, is given by a projection of A/P^"?' onto the plane dg = 0. In 
Figure 4.10c, the plotting range of the parameter o, is extended from 
one unidirectional state along a planar boundary of the orientation plane 
in Figure 4.2b towards a second unidirectional orientation state. This 
illustrates the continuity of the admissible ranges and the redundancy 
of the orientation plane. 


4.7 Anote on closure approximations 


For a given second-order orientation tensor N, closure approximations 
identify an unknown fourth-order orientation tensor N"nkew following 
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-0.1 5. "m 
0.0 i ET 
dp 05 oo 02 Pd 

h 


(a) (b) (c) 


Figure 4.10: (a), (b) Views on N Planer (c) Transition of admissible parameters of planar 
states between two unidirectional orientation states visualized by varying the parameter 
a1 in Equation (4.92) in the range —2/3 < ay < 2/3. 


a specific closure approach such that 
Nunknown p] = N (4.94) 


holds. Based on the quantity which is directly affected by a closure’s 
assumptions, closure approximations might be classified into three 
groups. The first group contains closures based on algebraic assump- 
tions. Closures based on assumptions on the material symmetry of 
orientation tensors form the second group. The third group contains 
closures which state assumptions on the FODF. The linear closure is 
given by 

Niinear (N) = Nis? + sym (dev (N) @1), (4.95) 


belongs to the first group and is defined in Han and Im (1999, Equation 
(13)) or Advani and Tucker III (1990). The representation in (4.95) benefits 
from the notation introduced in Equations (4.58) and (4.59). The linear 
closure states that the unknown fourth-order orientation tensor is linear 
in the second-order orientation tensor. This implies vanishing fourth- 
order harmonic part, i.e., dev (N'nk^own) = O. Figure 4.3 contains the 
linear closure as dashed light green line and clearly shows the limitations 
of this closure as the line is not completely within the admissible param- 
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eter space. The second group contains, e.g., orthotropic fitted closures 
which combine analytical insights on the material symmetry of fiber 
orientation tensors with fitting procedures reproducing selected flow 
simulation results. Examples are given in Cintra Jr and Tucker III (1995); 
Chung and Kwon (2001; 2002) and both the orthotropic fitted closure 
(ORF) (Cintra Jr and Tucker III, 1995) and the invariant-based optimal 
fitting closure (IBOF) (Chung and Kwon, 2002) are visualized in Figure 
4.3 using nilsmeyerkit and JulianKarlBauer (2021). The third group 
contains closures, e.g., Montgomery-Smith et al. (2011a); Górthofer 
et al. (2020), which enable exact reconstructions of the unknown FODF, 
based on the assumption that the ODF belongs to a special class of 
FODFs. Therefore, assumptions are made on the FODF and tensors 
of any order can be constructed based on the reconstructed FODF. 
Although this short note on closures does not claim to be exhaustive, 
reference is made to additional closures, e.g., VerWeyst et al. (1999); 
Montgomery-Smith et al. (2011b); Jack et al. (2010). Reproduction of 
the FODF by maximization of the information entropy, e.g., Müller 
and Bóhlke (2016) (see also Bóhlke (2005)) states an alternative and its 
wording explicitly addresses the problem of gaining information by 
assumptions. Assessment of closure approximations is challenging since 
the definition of an application-independent metric for the quality of the 
closure is difficult. Nevertheless, assessments of closure approximations 
are given in Breuer et al. (2019); Jack and Smith (2007; 2008). 


4.8 Summary and conclusions 


Application of linear invariant decomposition with focus on index sym- 
metry (Rychlewski, 2000) transfers the results of Kanatani (1984) into the 
continuum mechanics framework and leads to a compact representation 
of fourth-order fiber orientation tensors in Equations (4.58) and (4.59). A 
harmonic, i.e. completely symmetric and completely traceless, triclinic 
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fourth-order structure tensor is introduced in Equation (4.62). This struc- 
ture tensor has a simplified structure in the coordinate system spanned 
by the eigenvectors of the second-order fiber orientation tensor and leads 
to a parameterization of generic fourth-order fiber orientation tensors in 
Equation (4.61). Material symmetries of the orientation distribution func- 
tion reduce the number of independent parameters. In the triclinic case, 
the independent parameters are a second-order fiber orientation tensor 
plus nine scalars. Admissible ranges of the independent parameters are 
discussed in detail for orthotropy, transversely isotropy and planar cases. 
Insufficiency of coefficient-wise constraints are demonstrated for the 
transversely isotropic case. The variety of fourth-order fiber orientation 
tensors is given by the set of positive-definite tensors which can be 
expressed by Equation (4.61). Inspection of planar orthotropic states 
illustrates the character of fourth-order orientation information. Notes 
on closure approximations demonstrate their limitations on reflecting 
the variety of fourth-order fiber orientation tensors. 

As asummary, the parameterizations of generic second-order orientation 
tensors N 


N (a, a3) = Nis0 4 o4 Fransvt ED ag Dc see (4.29) 
generic fourth-order orientation tensors N 


N (N, di, ..., dg) = see (4.61) 


Niso EN ° sym (dev (N) & I + Fic (di, it's dg) ; 


orthotropic fourth-order orientation tensors NO"! 


Noho (N, di, da, d3) = see (4.78) 


Nise 4. sym (dev (N) & I) + ortho (di, da, da) 
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and planar fourth-order orientation tensors Nplanar 


NPlanar (o, di, dg) = (om (a1), di, see (4.91) 
2] = = 
ds = 504 2 ds = 15a4 d 
140 140 


0, 0, 0, 0, dg, do — -4) 


are repeated, including the structure tensors 


1 0 0 
piani _ —1/2 0 Vi 8 Vj, see (4.27) 
sym —1/2 
—1/2 0 0 
pans _ -1/2 0| v; vj, see (4.28) 
sym 1 
Fr (d4, ..., dg) = see (4.62) 
—(dı + d2) di d2 — V2 (d4 + d5) V2 de V2 da 
—(di + da) ds V2 da — V2 (de + d7) V2 dg 
—(da + da) | V2 ds V2 dz — V2 (da + do) 
completely symmetric 


tent 
and 
Fo (dy, da, d3) = IF"? (di, da, da, 0, 0, 0, 0, 0, 0). see (4.77) 


Admissible parameter ranges follow from the requirement of positive 
semi-definiteness (see Equation (4.64)) and are explicitly given for special 
cases of material symmetry in Equations (4.67), (4.87), (4.92) and (A.16). 
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Beside the review contribution, novel technical and scientific contribu- 


tions in this work are listed hereafter. 


An invariant framework for parameterizations of fourth-order fiber 
orientation tensors based on deviators from the isotropic orientation 
state is developed. The most general case is obtained by an orthotropic 
deviator of second order and an triclinic deviator of fourth order. 
Within this framework, deviators of stronger material symmetries are 
obtained as special cases with constrained parameters. 

Explicit parameterizations and admissible parameter ranges present 
in the literature are extended by orthotropic as well as planar fiber 
orientation tensors in a three-dimensional framework. 

Visualizations of admissible parameter sets are given for the first time 
for the orthotropic states in Figures 4.5 to 4.8 and for the planar states 
in Figures 4.10a. 

Positive semi-definiteness is identified to be a stronger constraint 
on admissible orientation tensors than linear constraints on tensor 
coefficients which themselves are identified to be stronger constraints 
than a bonded norm. 

A short hand notation for completely symmetric fourth-order tensors 
is introduced. 

Representations of discrete fiber sets (Figures 4.4 and 4.9) are identified 
to be an valuable visualization of orientation tensors. 


These results enable engineers to parameterize the space of admissi- 


ble fiber orientation tensors obtaining a valuable tool for engineering 


applications like model inspection or data validation. 
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Chapter 5 


Fiber Orientation Distributions 
Based on Planar Fiber Orientation 
Tensors of Fourth Order ' 


5.1 Introduction 


Suitable microstructure descriptors are essential for the prediction of 
effective properties of fiber reinforced composites. Fiber orientation 
distribution functions (FODF) represent exact microstructure descriptors 
of the orientation of axisymmetric fibers within a specified volume of 
a fiber reinforced composite (Kanatani, 1984; Advani and Tucker III, 
1987). However, in practice the exact distribution of fibers is commonly 
approximated by fiber orientation tensors (Bóhlke et al., 2019; Górthofer 
et al., 2019) which represent averaged directional measures and can 
be directly determined by non-destructive analysis methods, such as 
computer tomography (Pinter et al., 2018; Schóttl et al., 2020). Further- 
more, fiber orientation tensors fit into the tensor-based framework of 
continuum mechanics which is frequently used in material modeling 


1 This chapter reproduces (Bauer and Böhlke, 20222), i.e., Bauer, J. K., Böhlke, T., 2022. 
Fiber orientation distributions based on planar fiber orientation tensors of fourth 
order. Mathematics and Mechanics of Solids, online first, 10.1177/10812865221093958. 
Reproduced with permission. €2022 The Authors. Published by SAGE Publications Ltd 
under CC BY-NC 4.0 
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(Kehrer et al., 2020; Brylka, 2017; Hessman et al., 2021). The dependence 
of material models on the set of admissible second-order fiber orientation 
tensors is studied numerically in, e.g., Goldberg et al. (2017); Gorthofer 
et al. (2020). However, the dependence of material models or other 
quantities of interest such as reconstructed FODF itself on fourth-order 
fiber orientation tensors is rarely studied analytically. In addition, for 
some applications, e.g., damage modeling (Schemmann et al., 2018b) or 
averaging schemes (Hessman et al., 2021; Brylka, 2017), the identification 
of a FODF based on given fiber orientation tensors is beneficial. Due 
to the averaged character of fiber orientation tensors, no one-to-one 
correspondence to an FODF exists. Nevertheless, for a given fiber 
orientation tensor, any associated FODF is of interest and eases the 
indirect visualization of fiber orientation tensors. Many fiber-reinforced 
composites are plate-like and if the mean fiber length is larger than the 
plate thickness the resulting fiber distribution is approximately planar. 
This holds, e.g., for sheet molding compound (Górthofer et al., 2019). 
For such planar fiber distributions, Bauer and Bohlke (2022c) identifies a 
set of all admissible fourth-order fiber orientation tensors. Based on this 
set and reconstruction methods following Kanatani (1984); Miiller and 
Bohlke (2016), the variety of reconstructed FODF based on planar fiber 
orientation tensors of fourth order is studied in the current work. 

This paper is structured as follows: Definitions of FODF and fiber 
orientation tensors are followed by a reformulation of planar fiber 
orientation tensors following Bauer and Bóhlke (2022c). Within the 
admissible set of planar fourth-order fiber orientation tensors given by 
Bauer and Bohlke (2022c), a subset of distinct planar fourth-order fiber 
orientation tensors is identified. FODF approximations by truncated 
Fourier series with planar leading fiber orientation tensors in a three- 
dimensional framework are identified as non-planar and motivate a 
two-dimensional framework. The reconstruction of FODF based on 
the maximum entropy method following Müller and Bóhlke (2016) is 
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recast into this two-dimensional framework. Discrete slices and points 
of the set of admissible and distinct planar fourth-order fiber orientation 
tensors are used to visualize reconstructed FODF based on fourth-order 
fiber orientation tensors. A note on FODF reconstruction solely based 
on second-order fiber orientation tensors including the exact closure 
(Montgomery-Smith et al., 2011a) closes this paper. 


5.2 Directional measures as microstructure 
descriptors 


Taking the average of a tensorial quantity over orientations requires 
a directional measure which quantifies the orientations. Established 
directional measures of axisymmetric fibers are the fiber orientation 
distribution function (FODF) and fiber orientation tensors of several 
kinds and orders. Both FODF and fiber orientation tensors quantify 
orientations inside a reference volume which might be interpreted as 
a section of specified size at a specific position x inside a structural 
component. Consequently, the directional measurement depends on 
the size of this reference volume which usually represents a scaling 
parameter in measurement algorithms, see e.g., Górthofer et al. (2019, 
Figure 4) or Schöttl et al. (2020). This section introduces basic quantities 
and in large parts follows Bauer and Bóhlke (2022c, sections 2.1 and 2.2). 


5.2.1 Fiber orientation distribution function 


The fiber orientation distribution function w at a given position inside a 
component 


v:S? >R, with S? = (n € R? ||n|| 2 1] (5.1) 
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maps any direction n onto a scalar value y (n). S? is the two-dimensional 
surface of a unit sphere in three dimensions parameterized by, e.g., a 
unit vector n. The function v (n) is non-negative, i.e., 


w(n)>0, Ynes? (5.2) 


holds and normalization of v (n) implies 


f (n) dn — 1. (5.3) 

S2 

As fibers have a direction but no attitude, v/ (n) is symmetric, i.e., 
$(-n)-v(n), Ynes? (5.4) 

holds (see Advani and Tucker III (1987); Górthofer et al. (2020)). 

5.2.2 Fiber orientation tensors 


Fiber orientation tensors of Kanatani first kind (Kanatani, 1984) are 
defined by 


Nu, = f(a) n?* dn (5.5) 


with n®* being the k-th moment of n. In this work, only the second- and 
fourth-order fiber orientation tensors 


N = Ng = v (n)n @ndn, (5.6) 
S2 
N-Na- f v(n)n@n@n@ndn (5.7) 
S2 
are used. More details on the properties of N,,, can be found, e.g., in 


Bauer and Bohlke (2022c). Higher order fiber orientation tensors of 
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Kanatani first kind contain all tensors of lower order as 
Nx—2) = Na [I] (5.8) 


holds for 2 < k with the identity on second-order tensors I = 6;; ei & ej. 
In addition, for the second-order orientation tensor 


N-I=tr(N) =1 (5.9) 


holds due to normalization, see Bauer and Böhlke (2022c). Orientation 
tensors of Kanatani first kind are commonly used to represent direc- 
tional data obtained by computer tomography scans or results of flow 
simulations (Görthofer et al., 2019). Basic properties of the second-order 
orientation tensor N are briefly summarized following Bauer and Böhlke 
(2022c). N is symmetric and positive semi-definite. In consequence, 
N can be diagonalized, i.e., pairs of eigenvalues X; with A; > 0 and 
orthonormal eigenvectors v; for i € [1, 2, 3] exist, such that 


3 Aı 0 0 
N = Ni) ei & ej = 5 AiVá C) VQ = Aa 0 Vi Q Vj (5.10) 
i=l sym Aa 


holds and there exists a rotation defined by an orthogonal tensor 


mapping the arbitrary but fixed basis {e;} onto the basis {v;}. The 
orthonormal basis {v;} spanned by the eigenvectors is, in the following, 
called orientation coordinate system. Based on the ordering convention 


A3 € A2 < Ai, (5.12) 
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an established classification of structurally differing N exists (see Cin- 
tra Jr and Tucker III (1995); Chung and Kwon (2002); Goldberg et al. 
(2017); Kóbler et al. (2018)). A visualization corresponding to this 
classification exists and is called orientation triangle. Any second- 
order orientation tensor can be represented by a pair (Ai, A5), which is 
connected to a point inside the orientation triangle and by a mapping 
Q which defines the orientation coordinate system. The FODF can be 
expressed in terms of fiber orientation tensors as a three-dimensional 
tensorial Fourier series 


V(n)— = > an (x) dev (Na) x (5.13) 
k= 


which is called spherical harmonic expansion (Kanatani, 1984, page 
154). The operator dev(-) extracts the deviatoric part and its definition 
for higher order tensors in a three-dimensional framework is given 


1 
2 
1 
ae 
"eg 
6 
0 
d8 ` 0.2 
@ isotropic @ unidirectional (b) Visualization of the set of admissible 
© planar isotropic === planar planar fourth-order fiber orientation 
tensors VP!" in the parameter space 
(a) Orientation triangle visualizin A1, di, dg} using the parameterization 
8 8 8 the p 
second-order orientation tensors of Nplanar in Equation (5.15) 


Figure 5.1 
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in Spencer (1970). Equation (5.13) implies, that in general an infinite 
number of orientation tensors of increasing order is required to express 
w (n). However, for many practical applications only second and fourth- 
order orientation tensors are available at a specific location inside a 
component. Recent developments in computer tomography methods for 
fiber reinforced composites lead towards the identification of the FODF 
instead of fiber orientation tensors, see, e.g., Schóttl et al. (2021). 


5.3 Admissible and distinct planar fiber orien- 
tation tensors 


A parameterization of planar fiber orientation tensors of second order 
NPplanar is given by 


Planar (41) = Ar v4 Q vi + (1— A1) v2 @ v2 (5.14) 
à 0 0 
= 1—A; 0 | Vi ®v; 
sym 0 


in the orientation coordinate system {v;} defined in Equation (5.11) and 
with 1/2 < A, < 1. Bauer and Böhlke (2022c, Equation (87)) identify 
planar fiber orientation tensors of fourth order NPlar@ (a1, dı, ds) which 
can be reparameterized in A, with o; = $ (A1 — 4) leading to 


lanar 
NP??? (Ai, di, dg) = (5.15) 
à-d- dtm 0|0 0 vV2ds 
(1—31)-di—& 0|0 0 -V2ds 
35 
0|0 0 O0 
č 8 BY 
completely | symmetric 


97 


5 Fiber Orientation Distributions Based on Planar Fiber Orientation Tensors 


in Kelvin-Mandel (Thomson, 1856; Mandel, 1965; Mehrabadi and Cowin, 
1990) notation, which is explained in detail in Appendix B.1. The 
Kelvin-Mandel basis B? & BY is spanned in the orientation coordinate 
system {v;}, i.e., for example the fourth basis vector in Equation (B.1) 
becomes BY = v2 (v2 8 v3 -- va & vo). In Equation (5.15), a short 
hand notation for completely symmetric fourth-order tensors following 
Bauer and Bóhlke (2022c) is used, which is explained in Appendix B.1. 
Demanding positive eigenvalues of NP!*"" in Equation (5.15) leads to 


the set of admissible planar fourth-order orientation tensors 


Afplanar _ Lon (as di de) | Sm Ax (5.16) 
— f (ài, di) € dg € f (^, a) | 
with 


f(a, dı) = (5.17) 


3 y —16—280 dı —1225 d?+140 41 —140 A1?+1225 dı 41 —1225 di Aı?, 


see Bauer and Böhlke (2022c, Equation (89)). The set \/P!*"*" implicitly 
defines a body in the parameter space (A1, di, dg}. This body is visual- 
ized in Figure 5.1b and contains all admissible planar fourth-order fiber 
orientation tensors. Any point inside this body represents a fourth-order 
orientation tensor. 

As Bauer and Böhlke (2022c) derive N Plarar solely from algebraic prop- 
erties of NPlanar, no planar fourth-order fiber orientation tensor outside 
AfPlanar. exists, Consequently, the parameterization Planar (A1, di, ds) 
combined with the parameter space N Planar allows for a complete study 
of the influence of fiber orientation on the mechanical response of 
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material models, which are based on fourth-order fiber orientation 
tensors of planar fiber architectures. However, not necessarily all points 
inside N Plana represent orientation states which differ structurally and 


therefore lead to different mechanical behavior. Two tensors A and B are 


called structurally identical if they differ solely by a rotation, i.e., 


LU 


Q € SO(3) with Q x A = (5.18) 


holds with the special orthogonal group in three dimensions SO(3). The 
parameterization of NP^??' in Equation (5.15) and the corresponding set 
of admissible parameter combinations A/P'?^?r in Equation (5.16) is based 
on the orientation coordinate system which is introduced in Equations 
(5.10) to (5.12). However, for orientation states with planar isotropic 
second-order orientation tensor, i.e., ^; = 1/2, the eigenvalue problem 
identifying v; and vz is ill-posed. As the eigenvalues of N (A; = 1/2) 
are 1/2, 1/2 and 0, any pair of two orthonormal vectors being pairwise 
perpendicular to v3 is a valid choice for the eigenvectors. The ambiguity 
of the orientation coordinate system introduces a redundancy in the 
set of admissible fourth-order orientation tensors given in Equation 
(5.16), i.e., multiple parameter combinations lead to structurally identical 
Nplanar, The shape of the body A/P*?"?' in Figure 5.1b motivates the 
reparameterization of NPlanar (\,,dı,dg) as a function of A, f, B with 


dy = di + > [3 — 22] = (5.19) 
d, — £sin (8) (5.20) 
dg = f cos (3) (5.21) 


The parameters £ and are visualized in figure 5.2a and lead to a 


99 


5 Fiber Orientation Distributions Based on Planar Fiber Orientation Tensors 


reformulation of Equation (5.15) with 


NPlanar (s. f, 8) = (5.22) 


- 24% 0|0 0 V2rcos (8) 
-rsin (Â) + 4-3 +1 0.0 0 — 2f cos (8) 
0 


0 0 0 
completely symmetric 
BY © BY 
and Equation (5.16) by 
Nplanar _ {Nee On. f, 8) | ! <) <1, (5.23) 
oses 2 [s - A], 
ÜDege an} 


For the special case A; = 1/2, the reparameterized fourth-order fiber 
orientation tensor reads as 
Nplanar Qu = 1/2, f, 8) = (5.24) 
-fsin (8) +2 f sin (B + 0/0 0 vf cos (ô 
-rsin(3)+3 0/0 0 —Vv2cos (8) 
01:00 0 


completely symmetric 


z 8 Bc 


and comparison with representations of rotations in Kelvin-Mandel 
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notation in Cowin and Mehrabadi (1995, section 3) and Mehrabadi and 
Cowin (1990, section 3) indicate a rotational redundancy. Active rotation 
of NPlaner( A, = 1/2, 7, 3) by an angle of —5/4 around the axis v3 leads 
to 


Q" (y = - JA) « NPanar (1 = 1/2,58) (5.25) 
N Ài (x =1/2,7,8 = 0) 
3/8 1/8 0/0 0 Var 
3/8 0|0 0 -v2 
0 0 0 Bio BY 


completely symmetric 


with a rotation around the axis v3 parameterized by the angle y with 


cos(y) -sin(y) 0 
Q (y)= | sin(y) cos (7) 0 [vi v. (5.26) 
0 0 1 


As the angle B in Equation (5.25) is arbitrary and the right hand side 
of Equation (5.25) is independent of ĝ, any NPlner(A, = 1/2, f, 9) can 
be expressed by a reference tensor Nplanar Qu = 1/2,f, B = 0) combined 
with a rotation by 


Nplanar (Aı = 1/2, f, 8) - (5.27) 
Q“ (y= 3/4) x NPlanar (= =1/2,7 B=0). 


In consequence, coincidence or symmetry of the eigenvalues A; and Aa 
in the special case of \; = 1/2, degenerates the space of structurally 
differing planar fiber orientation tensors of fourth order from a circle to 
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(b) As an example, an admissible parameter 

triplet (A1 , d,, dg) is exemplarily selected to be 
(2/3, 1/[16 V2], 1/[16 /2]). The polar plot (1) 

(a) The parameter space NPlanar shows the projection of NPlanar (A! , d^ , dh) onto the 
spanned by (A1, di, dg) or vi-v2-plane, i.e., NPlanar (A. q^ , q^) NOT (y, 7/2) 
equivalently (A1, f B ) following with the fourth-order moment OR (p,0) ofa 
Equations (5.19) to (5.21). The line unit vector n(y, 0) in polar coordinates, e.g., 

(Ay = 1/2, f,0) with 0 € # < 1/8 defined in Equation (5.33). The polar plot (2) shows 
is colored in red. Nplanar(\ | q^. — d^) -n®4(p, 7/2). 


—i 
8 


Figure 5.2 


a line. Fiber orientation tensors along this line act as reference tensors 
following Equation (5.27). The line is, e.g., parameterized by f or 
equivalently by dg and one arbitrary choice is colored red in Figure 
5.2a. 

If Nplanar serves as a directional measure, the full parameter range of ds, 
specified in Equation (5.16), is necessary to account for all admissible 
orientation states. If, in contrast, the influence of NPlanar onto mechanical 
properties or reconstructed FODF is to be studied, a minimal parameter 
set of structurally differing planar fiber orientation tensors is of interest. 
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This set is given by 


. 1 
Afplanar = [wee (Ai, di, dg) | 5 «AXI, 
4 4 
-—<d<-24+h-A1, 
0 < ds < f (An, di) } 


u {Near (Ay = 1/2, di, ds) | dı =0, 0 < ds < 1/8} (528) 


with Nplanar and f (Ai, dı) defined in Equations (5.15) and (5.17) and 
f(Ar = 1/2, dı = 0) = 1/8. Equation (5.28) accounts for the special 
case A; = 1/2 and removes a second redundancy which is visualized 
in Figure 5.2b and derived by the following observation. Two fiber 
orientation tensors of fourth order NPlara (\,, dı, dg) which differ solely 
by the sign of the parameter dg are structurally identical and only differ 


by the rotation 
10 0 
(A mus 0 —1 0 Vi ® Vj (5.29) 
00 -—1 


which rotates any physical quantity by 180° around the axis vı. Due to 
its symmetries, for any planar fourth-order fiber orientation tensor 


gone x „ Nplanar Dus dı, dg) u Nplanar (Ai, di, —da) (5.30) 


holds and motivates restriction to positive values of dg in the minimal set 
specified in Equation (5.28). This set can be used to study the influence 
of planar fourth-order orientation tensors on derived quantities. 
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5.4 Reconstructed fiber orientation distribu- 
tion functions 


In the previous section distinct and admissible planar fiber orientation 
tensors of fourth order are identified. The question of interest is, which 
fiber orientation distributions are associated to these tensors? It is 
evident from Equation (5.13) that for a given leading fiber orientation 
tensor there is non one-to-one correspondence to an FODF. However, 
identification of any FODF which is connected to the given leading fiber 
orientation tensor is of interest. An approximation of a FODF by lead- 
ing fiber orientation tensors up to fourth order in a three-dimensional 
framework is given by 

) (n, N) = li des (N) n8? + ae (N):n9*|, (531) 

Am 2 8 

following Equation (5.13), where N can be expressed by N based on 
Equation (5.8). Müller and Bóhlke (2016) discuss that the approximation 
ý (n, N) is not necessarily non-negative, i.e., the condition stated on 7 (n) 
in Equation (5.2) does not hold for the approximation 7) (n, N), due to the 
truncation after the fourth-order term. For the special case of planar fiber 
orientation tensors NP!*"*" the approximation in Equation (5.31) leads to 


QV" Co, 9, X,, di, dg) = Ô (n (ip, 0) ‚Nplanar (A. di, ds)) (5.32) 
1 
UIT [201604 sin? (0) sin (4y) 
+ 3 [72 — 6720d;] sin? (0) cos (4) 
— 3 [1120A, — 560] [3 cos (20) + 1] sin? (0) cos (2) 


— 420 cos (20) + 945 cos (40) + 435 
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(a) (b) (c) 


Figure 5.3: Approximation of the FODF by leading fiber orientation tensors up to fourth 
order following Equation (5.32), i.e., QN (p, 6, NPlanar (A, di, ds)) for three selected 


points in MPlanar, Positive values are plotted in blue and negative values are plotted in red. 
The parameters in the format (A1, di, dg) are (a) (1/2, 3/280, 0) (b) (2/3, —1/315, 0) (c) 
G, —4/ 35, 0) 


with a parameterization of the unit vector n in two spherical angles 
n (9,0) = sin(0) cos(p)vı + sin(0) sin(p)va + cos(0)va. (5.33) 


The approximation pee is not planar Figures 53a to 5.3c show 
three-dimensional spherical plots of QN (0, p, A1, di, dg) for three 
selected combinations of the parameters A1, dı and dg. Representation 
of non-smooth functions by Fourier series require an infinite number 
of non-vanishing Fourier coefficients. For any given planar FODF, the 
transition from non-vanishing values in the plane spanned by v, and vz, 
to vanishing values outside this plane is not smooth. In consequence, this 
transition is inadequately represented by the truncated spherical Fourier 
series N”, As a result, planarity of a given fiber orientation tensor in 
a three-dimensional framework does not imply that derived quantities, 
such as reconstructed FODF, are planar as well. Following Bauer and 
Bóhlke (2022c, Equation (59)), any fourth-order fiber orientation tensor 
in three dimensions can be expressed by two deviators dev (N) and 
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dev (N) fulfilling 
N = N + sym (dev (N [I]) & I) + dev (N). (5.34) 


Vanishing deviators, i.e., dev (N) = 0 and dev (N) = 0 describe the 
isotropic state N*^ = -sym (I & I) which corresponds to the FODF 
Y (n, Niso) = — and which is not planar. For Nplanar the deviators from 
the isotropic state which are also utilized in Equation (5.31) read as 


Aı-1/3 0 0 

dev (NP) - 2/0- X, 0 vi@v;, (535) 
sym —1/3 

dev (e) E (5.36) 
à1/7-d -4 di 1/35 — Aı/7 0 0 V2ds 
(1—21)/7—di — 1/35 1/35—(1—31)/7 | 0 0 -V2dg 
3/35 0 0 0 
completely | symmetric 


BY @ BY 


and demonstrate the drawbacks of expressing planar fiber orientation 
tensors in a three-dimensional framework since the planarity has to be 
enforced by deviation from the isotropic state. Note that the fourth-order 
deviator depends on the parameter A, which defines the second-order 
fiber orientation tensor. 


5.4.1 Transition from 3D into 2D 
Motivated by the previous section and following Kanatani (1984); Bli- 


nowski et al. (1996); Vianello (1997); Desmorat and Desmorat (2015); 
Vannucci (2018), planar quantities are expressed in a two-dimensional 
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framework with representations 


2 
a= p» a,e; = a,e1 + (5€2 +0 e3, (5.37) 
i=1 
2 2 
A=),) Aue: ® 8; (5.38) 
i=1 j=1 
2 2 3 3 
= S Aijei @ ej +) 0eiSes + M 0esGe;, 
isi j=l i=] j=l 
3 3 3 3 
A= A,B; e B, = 5 5 AgcBe 8 Be (5.39) 
€=1¢=1 £c[1,2, 4] CE[1, 2, 4] 


with generic tensors a, A, A. This notation directly connects objects 
in R? and R? and requires the reader to select the appropriate two or 
three dimensional view onto the planar physical quantity of interest, 
which is part of the three-dimensional reality. For the two-dimensional 
Kelvin-Mandel bases B, = Bi, B, = B» and B, = B, holds. To be 
explicit, for indices of tensor components in the 2D-framework i, j € 
[1, 2] and £, ¢ € [1, 2, 3] holds in contrast to i, j € [1, 2, 3] and £, C € 
[1, 2, 3, 4, 5, 6] in the 3D-framework. 

If planar orientation tensors are derived from a two-dimensional frame- 
work, naturally no out-of-plane tensor components exist. Basic planar 
isotropic tensors in the 2D framework are given by 


1 0 
I= 0;; e; ® ej = | san, 4 e; 8 £j, (5.40) 
1 0 0 
I? = õe B, 6B, = 1 0 | Be@By, (5.41) 
sym 1 
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and 


P; = 


IGI, P, =I” — P, (5.42) 


N| = 


following Blinowski et al. (1996, Equation (2.3)), Aßmus et al. (2017, 
Equation (38)). Deviator operators are defined by 
dev (A) = A— 5 (A - DL, (5.43) 
dev (A) = sym (A) — sym (sym (A) [I] @ I) (5.44) 


ak iym (I6 I) (I- A |I]) 


following Kanatani (1984, (7.13)) and Vianello (1997, section 3). The 
central fiber orientation tensor in the two-dimensional framework is 
identified as 


. D p 1 1 
Nie _ | =- a1) P, + = n°*) P ,=3P +3 (545 
a 


i 3 1 
=> 3 J0 B; 8 Be 
compl. | sym. 


= Ner, = 1/2, dı = 0, ds = 0) 


and is called planar isotropic fiber orientation tensor of fourth order. It 
should be noted that this is not the only fourth-order fiber orientation 
tensor which contracts to the planar isotropic second-order fiber orienta- 
tion tensor. Deploying harmonic decomposition in the 2D-framework 
following Blinowski et al. (1996, Equation (2.25)) or Desmorat and 
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Desmorat (2015, Equation (18)) with 


A= (50 ale + (2h -a)B tsymidev(AN)oD 649 


Pj 
+ dev (A) 


and knowledge on irreducible tensors, any fourth-order fiber orientation 
tensor is parameterized by 


N (Au, pi, p2) = NPS? + sym (E (1) @ D) + E (pı, p2) (5.47) 
with 
1 0 
F (Ai) = (Ai - 1/2) | V; 8 Vj, (5.48) 
sym —l 
—] 1 0 
F (p1, p2) = pı =1 |0 Bz @ BY (5.49) 
compl. | sym. 
00 [V2 
+ p2 0 -Y2 | BY 8 BY. 


compl. | sym. 


Coincidence with the parameterization in the 3D-framework is given 
by the prefactor in Equation (5.48) and the specific choice of the factors 
pı = dı — 3/280 and pə = ds leading to 


N (Ai, di, dg) = (5.50) 
ài — dı — 4/35 dı + 4/35 V2dg 
(1— 1) — di — 4/35 | -V2ds | BY @ BY. 
compl. | sym. 
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The shift by 3/280 is required due to the different expansion points in 2D 
and 3D. As the point of expansion in Equation (5.47), i.e., NPP9 is planar 
itself, the deviators are compact and formatted into the 3D-framework 


read as 
Aı-1/2 0 0 
dev (N) — 1/2— àı 0 | vj vj, (5.51) 
sym 0 
dev (N) — (5.52) 
280 dı (80 d) 0/00 V2 ds 
xg di 0/0 0 -vV2ds 
0 0 0 7 " 
completely symmetric 


The deviators specified in Equations (5.35), (5.36), (5.51) and (5.52) 
directly enter approximations of the FODF in terms of truncated Fourier 


series. 


5.4.2 Truncated fiber orientation distribution function in 
a 2D-framework 


In analogy to the three-dimensional case in Equation (5.13), FODF can 
be reconstructed in a two-dimensional framework. A FODF is given in 
terms of fiber orientation tensors by a tensorial Fourier series (Kanatani, 
1984, page 158) 


oo 


v (n) = = Y dev (Nyy) n° (5.53) 


k=0 
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with a unit tensor of first order n and fiber orientation tensors of k-th 
order in two dimensions N,;,. An approximation of the FODF y (n) by 
leading fiber orientation tensors up to fourth order is given by 


^ 1 
V (n,N) = = [1 + 4 dev (N) -n® + 16 dev (N) - n?^] , (5.54) 
and combined with N (A1, p1, p2) and n(p) = cos(p) v, +sin(y) v4 leads 
to 
" 1 
È (P; Ar, pi Ba) = z= [1 + (4 Ar — 2) cos (29) (5.55) 


— 16 pı cos (4y) + 16 pa sin (4p)]. 
Introducing the coordinate transformation 
p; — fsin (8) , Pa =fcos (3) (5.56) 


yields the following formulation 


Q (e. ED 8) = = E d (^ = 5) cos (2y) (5.57) 
— 16 sin (3 = 4e)] 


The parameterization of [D in Equation (5.57) separates second and 
fourth-order contributions into one trigonometric summand each. The 
cosinus-contribution of the second-order fiber orientation tensor is al- 
ways aligned with the coordinate axes of the orientation coordinate 
system and its amplitude scales with A;, being zero for the planar 
isotropic case, i.e., Aı = 1/2 and the frequency is two in the interval 
p € [0, 27). The fourth-order summand is a sinusoidal function of 
double the frequency of the second-order term, with phase shift B 
and an amplitude which scales with f. The frequency implies that 
this fourth-order term has four maxima and four minima within the 
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parameter range of o € [0, 27). In consequence, in the planar isotropic 
case, Equation (5.57) degenerates to 


) (v. = 2 f, 5) = = E — 16f sin (à = 4e)] (5.58) 


which directly indicates, that the parameter B implies a rotation of D by 
an angle of B /4 around the axis v3 of the 3D-framework pointing out 
of the 2D-plane. If ^; # 1/2, the interference of the second- and fourth- 
order summands, which have different frequencies, leads to the variety 
of y. The structure of Equation (5.57) shows, that ) may be represented 
by circular harmonics. Further details on harmonic representations of 
planar tensors of fourth order can be found, e.g., in Vannucci (2018, 
chapter 4) and Forte and Vianello (2014). 


5.4.5 Maximum entropy reconstruction 


Identification of a representative FODF based on a given leading fiber 
orientation tensor is of importance for several applications, e.g., numeri- 
cal calculation of orientation averages of direction dependent mechanical 
properties (Hessman et al., 2021; Schemmann et al., 2018b; Brylka, 2017). 
Truncated Fourier series are used in the literature to identify an FODF 
based on fiber orientation tensors (Schóttl et al., 2020; Jack and Smith, 
2004; Eik et al., 2016) although the identified functions do not meet the 
non-negativity requirement of a FODF. Naive averaging with partly 
negative FODF leads to non-physical results. Müller and Bóhlke (2016) 
give a solution to the reconstruction problem and identify FODF based 
on leading fiber orientation tensors by maximizing the information- 
theoretic entropy fulfilling normalization and non-negativity constraints. 
For limited available information, the entropy principle yields the most 
likely FODF fulfilling specified constraints. In this context "most likely" 
corresponds to "maximizing entropy". The procedure of Müller and 
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Bóhlke (2016) is briefly repeated for the special case of fourth-order fiber 
orientation tensors and tailored to planar fiber orientation tensors. The 
maximum entropy approximation p(n) maximizes the information- 


theoretic entropy E(w) 
Ei) =— [bay inca) as (5.59) 
and fulfills the constraints 
G=Geos = = w(n) dS — 1, (5.60) 
G = Geos = = i/(n) dev (n8?) dS — dev (N), (5.61) 
G=Ger = w ý(n) dev (n®) dS — dev (N) (5.62) 


where the constraint G ensures the normalization and G and G ensure 
that the solution 7"*(n) meets the given orientation tensors. The domain 
of the integral depends on the dimensionality of the given orientation 
tensor and is either the surface of the unit sphere 5? for planar 3D tensors 
or the circumference of the unit circle S for 2D tensors. The problem 


max E() with Gz4. =0<a> foro c[0,2, 4] (5.63) 
Y 


corresponds to the Lagrange functional 


LD) = E(w) = 5 Lea> + Gea> (5.64) 


a€ [0, 2, 4] 


with Lagrange multipliers L = Leos, L = Leos and L = L<4>. The 


solution is given by the root of the first variation of £, i.e., Ó£(v) = 0 
and is formulated in the Lagrange multipliers by 


Ņ™ (n, L,L,L) = exp(-1— L — L:dev(n9?) — L- dev(n®*)) (5.65) 
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^ A 


19/140 19/1404 
gd 6/6 d, ? 6/6 
5/6 5/6 
—4/35 4/6 Ai —4/35 4/6 Ar 
3/6 3/6 
° dg 1/8 j dg 1/8 


(a) Paths and points in the plane A; = 3/6 (b) Points in the plane \; = 4/6 used in 
used in Figure 5.8 figure 5.9 


Figure 5.4: Definition of representative points in the parameter space VP!" of planar 
fourth-order fiber orientation tensors. The outer points in radial directions do not lie on 
the boundary of the admissible region. Parameters of each point are listed in Table B.1. 


which fulfills the non-negativity condition identically. The Lagrange 
multipliers are obtained by solving the system of equations stated by 
the constraints in Equations (5.60) to (5.62) with 7(n) being repaced 
by Ņ™ (n, L, L, L) from Equation (5.65). Due to the projection of the 
Lagrange multipliers onto deviators of the moment tensors in Equation 
(5.65), the Lagrange multipliers are without of loss generality irreducible, 
i.e., completely symmetric and traceless (Rychlewski, 2000; Spencer, 
1970). Table 5.1 lists the number of independent components of irre- 
ducible tensors which depend on the tensor order and the dimenision- 
ality of the vector space. The number of independent components of 


Tensor order || 2 
2 2 
4 2 


Table 5.1: Number of independent component of irreducible tensors depending on order 
and dimensionality of the vector space 


D | 3D 
5 
9 
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the Lagrange multipliers L, Land L sum up to 1 +5 + 9 = 15 and 
1+2-+2= 5 in three and two dimensions respectively and correspond 
to the fourteen independent components of N and four independent 
components of N plus one normalization constraint each. In tensor 
notation, Equations (5.60) to (5.62) read as 1+9+81 = 91 scalar equations, 
whereas in symmetric Kelvin-Mandel notation the identical equations 
yield 1 + 6 + 36 scalar equations. These naive views on the set of equa- 
tions illustrate that a minimal and redundancy-free parameterization 
deploying index symmetry of the Lagrange multipliers is required. Such 
a parameterization is necessary to reduce the dimensionality of the 
system of equations to at most 15 or 5 depending on the dimensionality 
of used tensor framework. Three out of five independent components 
of N in 3D and one out of two independent components of N in 2D 
define the orientation coordinate system. In the orientation coordinate 
system, the second-order orientation tensor has diagonal structure. Off- 
diagonal components of the second order Lagrange multiplier ensure 
that the orientation coordinate systems of the reconstructed FODF and 
the original fiber orientation tensors coincide. Based on the orientation 
coordinate system, a generic representation of the Lagrange multipliers 
of order two and four in three dimensions is given by 


| bi bs b 
L(bi, m bs) = ba b3 Vi & Vj, (5.66) 
sym — (bı + bz) 
L(cı, ee, C9) = (5.67) 
—(e1 te2) c1 c2 —V2(c4 +c5) V2cg V2 cg 
—(c1 +¢3) ea V2 c4 — v2 (cg +c7) V2cg 
—(c2 + e3) 2 c5 V2 c7 — V2 (cg + c9) 
completely | symmetric 


By @ By. 


Starting from the triclinic case in Equations (5.66) and (5.67), material 
symmetries may be deployed to further reduce the number of indepen- 
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dent components of the Lagrange multipliers, as the fiber orientation 
tensors and the reconstructed FODF share their symmetry group, i.e., 
S" = S with 


$(Q «n) = ó(n) vQes*, (5.68) 
Q+N=N vOcs" (5.69) 


which follows directly from the definition of N in Equation (5.7). Similar 
to the FODF reconstruction by truncated Fourier series, the ansatz ("* 
in Equation (5.65) contains deviators from a central state, which in 3D is 
the isotropic state. In order to reconstruct planar FODF based on planar 
fiber orientation tensors, it is beneficial to deploy a two-dimensional 
framework with deviators deviating from the planar isotropic state. In 
the orientation coordinate system the Lagrange multipliers of order two 
and four read as 


L (fi, fe) = | fi E Vj @V;, (5.70) 
sym -fı 
-91 gı V2 92 
L(g, 92) = -gı | -V292 | BZ @ BŽ. (5.71) 
compl. | sym. 


In consequence, combination of Equations (5.70) and (5.71) leads to 
$^ (o. L, fi, fa, 91, 92) = exp( —1- L- ficos (29) (5.72) 
— fasin (2) + gi cos (4g) — g2 sin (4). 


For a given N(A;, di, ds) following Equation (5.50) and a specific se- 
lection of independent tensor components, the resulting five scalar 
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equations are 
n VO L, fi, fa, 91,92) dS — 1 — 0, (5.73) 
sus; L, fi, fo, 91, 92) (30 o) ds = (^ = 5) =0, (5.74) 


(n, L, f ‚f2 91, 92) sim 20) dS — 0 = 0, (5.75) 


a go a a 
& 


€ 


de L, fı, f2, 91, 92) (poao) ds — (« i= 3) =9 ‚ (5.76) 
( 


8 
(n, L, f , f2, 91, 92) v2 sin 0) ds — V2 dg = 0. (5.77) 


e 


The equations are solved numerically for a given tripled (Ai, dı, dg). 
Numerical integration on the unit circle is done with nschloe (2021) 
based on Krylov and Stroud (2006). 


5.4.4 Visualization of reconstructed planar fiber orienta- 
tion distribution functions 


Fiber orientation tensors represent averaged properties of an underlying 
FODE. In consequence, a complete reconstruction of the underlying 
FODF is not possible. It is of interest to visualize possible shapes of 
reconstructed FODFs based on admissible mean values. The admissi- 
ble mean values are given by the set of distinct and admissible fiber 
orientation tensors of fourth order A/ Planar which: can be combined with 
the developed reconstruction methods ) and $ °. The parameter space 
AFP is discretized and for selected points in VP", i.e., for selected 
N, reconstructed FODFs are visualized. The discretization is based on 
several slices through N Pla" which are given in Figures 5.4 and 5.5. 
The arrangement of the visualizations of the reconstructed FODFs in 
figures 5.7, 5.8, 5.9 and 5.10 mimic the position of the points on the slices 
and reuse the color coding of the overview plots. In consequence, the 
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19/1404 


dy 


EIE 
—113/2520 | 
0 
—4/35* ns 
A F mra 1/8 
3/6 4/6 5/6 6/6 
Ai ^ ds 


(b) Ten points in the plane of planar orthotropic fiber 
(a) Points in the plane orientation tensors of fourth order being a slice of 
Aı = 5/6 used in figure 5.10 NM Planat, The points are used in figure 5.7 


Figure 5.5: Definition of additional representative points, see also Figure 5.4 


influence of the position inside VP*"* on the shape of the reconstructed 
FODF is obtained. The developed visualizations directly show the 
variety of reconstructed FODFs and as a side effect act as an intuitive 
view on fiber orientation tensors of fourth order. As the truncated FODF 
reconstruction does not fulfill the non-negativity constraint, negative 
values of ý are highlighted in red, as defined in the shared legend in 
Figure 5.6. Parameters of the utilized orientation tensors are listed in 
Table B.1. The scaling of all polar plots of FODFs in Figures 5.7, 5.8, 5.9 
and 5.10 is homogeneous with the lower radial limit being zero and the 
upper radial limit 0.9. 

A slice containing orthotropic (Bauer and Bohlke, 2022c) fiber orientation 
tensors is shown in Figure 5.5b and the corresponding FODFs are given 
in Figure 5.7. The transition from the planar isotropic state in Figure 


— Positive values of ) (p) — Negative values of ) (y) 


Figure 5.6: Shared legend for figures 5.7, 5.9, 5.9 and 5.10. ) (p) and p(y) are defined in 
Equations (5.54) and (5.72) respectively. 
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5.7(4) towards the unidirectional state in Figure 5.7(10) demonstrates 
the degeneration of the variety of fourth-order fiber orientation tensors 
towards the unidirectional state, which is defined completely by ^, = 1. 
Figures 5.7(1) to 5.7(3) correspond to orientation states represented 
by orange points in Figure 5.5b. These points keep a distance in the 
dı-direction from the boundary of the admissible body, as orientation 
states at the boundary itself, correspond to Dirac distributions into two 
directions each (see Bauer and Bóhlke (2022c, Figure 9)). In analogy, 
Figures 5.7(7) to 5.7(9) correspond to points highlighted in magenta in 
Figure 5.5b and are also located off the boundary of the admissible region. 
Figure 9 in Bauer and Bóhlke (2022c) contains FODF approximations 
based on a limited number of discrete fibers aligned along the axes of 
the orientation coordinate system which correspond to points on the 


Figure 5.7: Reconstructed FODF at specific points which are defined in Figure 5.5b. A 
legend is given in Figure 5.6 
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Figure 5.8: Reconstructed FODF at specific points along paths defined in Figure 5.4a. The 
radial spacing of the plots along the green and blue paths are not to true scale to allow for 
larger plots A legend is given in Figure 5.6 
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(02) ‘ee » 


Figure 5.9: Reconstructed FODF along several paths. The radial spacing of the plots along 
the green and blue paths are not to true scale to allow for larger plots. A legend is given 
in Figure 5.6 and parameters of the points in MPlarar, j e,, Nplanar, are listed in Table B.1. 
Paths are defined in Figure 5.4b. 
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Figure 5.10: Reconstructed FODF along several paths. The radial spacing of the plots along 
the green and blue paths are not to true scale to allow for larger plots. A legend is given 
in Figure 5.6 and parameters of the points in M planar i e, Nplanar, are listed in Table B.1. 
Paths are defined in Figure 5.5a. 
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boundary of the admissible region specified by the line di = —4/35 
and dg = 0 in Figure 5.5b. Planar fiber orientation tensors at the 
boundary of the admissible region directly correspond to orientation 
states of a limited number of Dirac distributions. This observation 
explains why experimentally identified fiber orientation tensors obtained 
from comparably large fiber arrangements may not fill the complete 
admissible region but fill only a (equally shaped) subregion within 
the admissible region. In a comparably large fiber arrangement, Dirac 
distributions are unlikely. The entropy reconstruction scheme contains 
a basis function involving a direction depended exponential function. 
Although the slope of an exponential function is rather steep, extreme 
values of the Lagrange multipliers are necessary to approximate Dirac 
distributions. Therefore, the entropy reconstruction scheme, which 
includes numerical integration and a numerical solver, is not completely 
stable when operating on orientation states at the boundary of the 
admissible region. As the unidirectional case in Figure 5.7(10) lies 
on the radial boundary of the admissible orientation states, the cor- 
responding entropy-based reconstruction FODF degenerates to a single 
Dirac distribution. The approximation $ in this figure is localized, i.e., 
vanishes for all angles except those pointing along axes of the orientation 
coordinate system. 

The slices defined in Figures 5.4, 5.5a and the corresponding FODF 
approximations in figures 5.8, 5.9 and 5.10 visualize the variety of 
reconstructed FODF among fourth-order fiber orientation tensors which 
contract to identical second-order fiber orientation tensors. Any FODF 
in Figure 5.8 leads to the planar isotropic second-order fiber orientation 
tensor N(Xı = 1/2). In analogy, any FODF in Figure 5.9 has an identical 
second-order fiber orientation tensor being N(Xı = 2/3). From this 
viewpoint, the figures in any of the figures 5.8, 5.9 or 5.10 visualize 
the variety of planar fourth-order fiber orientation tensors for a fixed 
second-order orientation tensor. Figure 5.8 shows the redundancy inside 
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M tt ^i 


4/8 5/8 6/8 7/8 9999/10000 
— Positive values of (p) — Negative values of 1}(y) 
Gestell je ( p) = =s je closure ( » 


Figure 5.11: Influence of Aı on reconstructed FODF for selected reconstruction methods 
based on second-order fiber orientation tensors. For Aj = 1/2 all approximations coincide. 


NPlanar identified in Equation (5.27) as FODFs with identical radial 
distance to the centered Figure 5.8(x) only differ by a rotation. In contrast, 
interference of the non-vanishing second-order contribution with the 
fourth-order contribution in figures 5.9 and 5.10 leads to a variety of 
FODF approximations. The developed views on JVP!*"*" based on slices 
and selected points can be reused to analyse and visualize any quantity 
which depends on planar fourth-order fiber orientation tensors. 


5.4.5 Reconstruction solely based on second-order fiber 
orientation tensors 


If fourth-order fiber orientation information is not available, the param- 
eter space of reconstructed FODFs degenerates to one parameter, e.g., 
1/2 € 4 € 1, which defines the second-order fiber orientation tensor. 
Missing fourth-order fiber orientation information implies vanishing 
fourth-order deviators and therefore leads to FODF reconstructions with 


v (v, Ai, Dpi-— 0, P2 = 0) in Equation (5.53) and $" (v, L, fi, fos 91 = 
0,g2 = 0) in Equation (5.72). Figure 5.11 shows the influence of A, 
on reconstructed FODFs for the methods used in the previous section 


and in addition for the FODF approximation obtained by the exact 
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closure approximation (Montgomery-Smith et al., 2011a). The exact 
closure approximation is frequently used, e.g., in Goldberg et al. (2017); 
Schneider (2017); Hessman et al. (2021); Górthofer et al. (2020) and 
postulates a one-to-one correspondence between a given second-order 
fiber orientation tensor and an a priori unknown FODF 7% “lesere of 
special shape. In the general three-dimensional case this FODF can be 
identified by solving elliptic integrals (Montgomery-Smith et al., 2011a). 
For the planar case, Górthofer et al. (2020) gives explicit formulas which, 
combined with Equation (5.14), lead to 


oe son Ad) = > (1 = Aı)Aı A (5.78) 
T A? + (1 — 2A1) cos (p) 

Visualization of Equation (5.78) for five values of ^, is given in Figure 
5.11. If missing fourth-order orientation information is compensated by 
deploying closure approximations, the variety of resulting FODFs shown 
in Figures 5, 5.9 and 5.10 is not reachable. Instead, the parameter space 
of structurally differing resulting FODFs degenerates to the parameter 
space of the second-order orientation tensor within the orientation 
coordinate system. If closure approximations are used in the context 
of planar fiber orientations, the complete orientation information is 
contained within two scalar parameters. One scalar parameters specifies 
the orientation coordinate system and the other parameter specifies the 
structure of the second-order fiber orientation tensor. In contrast, if the 
discrete orientation of each single fiber is known, e.g., based on direct 
fiber simulation (Meyer et al., 2020b), the variety of resulting FODF is 
generic. However, averaged information of such a discrete set of fibers 
expressed by fiber orientation tensors is still limited. 
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5.5 Summary and conclusion 


Planar fourth-order fiber orientation tensors describe the fiber orienta- 


tion in many sheet-like fiber reinforced composites. The variety of these 


tensors is known (Bauer and Bohlke, 2022c). The variety of reconstructed 


fiber orientation distribution functions (FODFs) based on planar fourth- 


order fiber orientation tensors is studied in this work leading to the 


following insights: 


1. 
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Based on the set of admissible planar fiber orientation tensors in 
Bauer and Bóhlke (2022c), a minimal set of admissible and distinct 
planar fiber orientation tensors of fourth order VP!" is derived 
and given in Equation (5.28). This set is the basis for studies on the 
influence of planar fourth-order fiber orientation tensors on any 
physical quantity or material function which is formulated in this 
directional measure. 


. The variety of reconstructed FODF is visualized by an arrangement 


of polar plots which mimics the shape of the admissible and 
distinct parameter space. This arrangement or view is generic 
and may be applied to study the dependence of other quantities 
on planar fourth-order fiber orientation tensors. 


. Reconstructed FODF based on truncated series expansion within a 


three-dimensional framework are identified to be not planar as the 
central state in three dimensions is isotropic. 


. A two-dimensional formulation of planar fiber orientation tensors 


is introduced and linked to parameterizations of planar fiber ori- 
entation tensors in three dimensions. The central, i.e., “isotropic”, 
state in two dimensions is planar isotropic. 


. Within the two dimensional framework it is shown that interfer- 


ence of second and fourth-order contributions leads to the variety 
of reconstructed FODF based on truncated Fourier series. 


. Visualizations of truncated FODF reconstructions in Figures 5.7, 5.9 


5.5 Summary and conclusion 


and 5.10 highlight their limitations and motivate more advanced 
reconstruction methods. The maximum entropy reconstruction 
of Müller and Bóhlke (2016) is explicitly formulated for the gen- 
eral three-dimensional case and recast for the planar case in a 
two-dimensional framework which leads to a low-dimensional 
optimization problem. Resulting FODF approximations are nor- 
malized and non negative. 

. For given reconstruction methods, the structural variety of re- 
constructed FODF based on planar fourth-order fiber orientation 
tensors is limited. Throughout this work, separation of rigid-body 
rotations, and thus orientation in space, from structural properties 
of represented FODFs is accomplished by representations in the 
orientation coordinate system. 

. Visualization of reconstructed FODFs solely based on second-order 
fiber orientation tensors, including a reconstruction method based 
on the exact closure (Montgomery-Smith et al., 2011a; Górthofer 
et al., 2020), close this paper and motivate higher order directional 


measures. 
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Chapter 6 


On the Dependence of Orientation 
Averaging Mean Field 
Homogenization on Planar 
Fourth-Order Fiber Orientation 
Tensors ' 


6.1 Introduction 


The effective mechanical properties of fiber reinforced composites highly 
depend on the microstructure. The microstructure descriptors fiber 
volume content and fiber orientation tensors are commonly used in 
homogenization techniques for two-phase composites (Buck et al., 2015; 
Brylka, 2017; Hessman et al., 2021; Schemmann et al., 2018b; Kehrer 
et al., 2020). Experimental or numerical identification and subsequent 
use of second-order fiber orientation tensors is well established in pro- 
cess chains (Górthofer et al., 2019) and commercial software. The use 
of higher-order fiber orientation tensors or the direct use of a fiber 


1 This chapter reproduces (Bauer and Bóhlke, 2022b), i.e., Bauer, J. K., Bóhlke, T., 
2022. On the dependence of orientation-averaging mean field homogenization on 
planar fourth-order fiber orientation tensors. Mechanics of Materials 170, 104307, 
10.1016/j.mechmat.2022.104307. Reproduced with permission. ©2022 The Authors. 
Published by Elsevier Ltd under CC BY 4.0 
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orientation distribution function, is progressing (Meyer et al., 2020b). 
Nevertheless, closure approximations, which identify higher-order 
orientation information based on assumptions instead of information, 
are still used in material modeling. For the special case of linear elasticity, 
mean field homogenization is popular due to its simplicity (Kehrer 
et al., 2020). A unified formulation of several mean field approxi- 
mations is given in Hessman et al. (2021). The influence of model 
parameters, the fiber volume content or the aspect ratio of either the 
fiber’s shape or the two-point correlation functions, on the predicted 
effective mechanical properties obtained by mean field homogenization, 
is studied extensively, e.g., in Kehrer et al. (2020); Buck et al. (2015); 
Brylka (2017); Hessman et al. (2021); Müller (2016); Trauth et al. (2021); 
Kehrer (2019). However, the influence of fiber orientation distributions 
on the effective properties has not been studied systematically. In this 
work, the influence of fiber orientation distributions on the effective 
properties of long-fiber-reinforced composites with large aspect ratios, 
such as sheet molding compound (SMC), is studied. Fiber distributions 
of SMC are identified to be in most cases approximately planar and 
therefore correspond to fiber orientation tensors which are planar as 
well. A set of all admissible and structurally differing planar fiber 
orientation tensors of fourth order is identified by Bauer and Bohlke 
(2022c;a) within a three-dimensional framework. This set is the basis for 
a complete study on the influence of the fiber orientation distribution 
on the effective mechanical properties for selected orientation-averaging 
mean field homogenizations within this work. This complete study is 
new and of high relevance as it visualizes all possible effective stiffness 
tensors which can be obtained for planar distributions and selected 
homogenizations. The methodology may be applied to any other ho- 
mogenization or any effective quantity, which is a function of a planar 
fourth-order fiber orientation tensor. As a side effect, the limitations 
of the use of second-order orientation tensors only are visualized. The 
developed visualization method may be adopted to higher-dimensional 
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parameter spaces of fourth-order fiber orientation tensors with lower 
symmetry. 

This paper starts with definitions of fiber orientation tensors before 
orientation states of SMC specimen obtained by computer tomography 
(Pinter et al., 2018; Schottl et al., 2020) are identified to be approximately 
planar. In consequence, the microstructure of SMC can be parameterized 
in the three-dimensional space of admissible planar fourth-order fiber 
orientation tensors. A parameterization of this space is obtained in 
Bauer and Bóhlke (2022c) and combined with a minimal parameter 
set of structurally distinct planar fourth-order fiber orientation tensors 
given in Bauer and Bóhlke (2022a). As the Advani-Tucker orientation 
average (Advani and Tucker III, 1987) is an essential building block of 
various mean field schemes, but the implication of its linearity in both 
the orientation tensor and the averaged quantity is seldom addressed 
explicitly, a reformulation based on the harmonic decomposition (lin- 
ear invariant decomposition) is given. A new numerical variant of 
the Advani-Tucker orientation average based on reconstructed FODF 
following the maximum entropy method, developed in Müller and 
Böhlke (2016); Bauer and Böhlke (20222), is given and shown to have 
a good convergence for non-localized orientation states. In a next 
step, four mean field homogenization approximations (Kehrer et al., 
2020; Walpole, 1966a;b; Benveniste, 1987) are reviewed, partially re- 
formulated and investigated. Implementations are made available at 
Bauer (2022). The first approximations are orientation-averaging Mori- 
Tanaka following Benveniste (1987) and a two step Hashin-Shtrikman 
homogenization scheme (Kehrer et al., 2020) both in formulations based 
on the orientation average following Advani and Tucker III (1987). The 
remaining two approximations are direct Advani-Tucker orientation 
averages of either a unidirectional stiffness or compliance obtained 
by the Mori-Tanaka approximation (Mori and Tanaka, 1973). Views 
on the effective stiffnesses obtained by the selected schemes are given 
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based on Young’s modulus and generalized bulk modulus following 
Bohlke and Briiggemann (2001). Observations on structural differences 
of the effective stiffnesses obtained by schemes averaging stiffness- or 
compliance-like quantities close this work. 


6.2 Sheet molding compound and planar mi- 
crostructures 


Sheet molding compound (SMC) is a material class with a thermoset 
polymer matrix enforced by long glass fibers and is of special interest 
in this work. Due to the manufacturing process (Bohlke et al., 2019), 
components made from SMC are shell-shaped, i.e., at each point of the 
component, the thickness is significantly smaller than the remaining two 
dimensions of the component. As the fiber length exceeds the component 
thickness, alignment of fibers perpendicular to the local plane of the 
shell is limited to local fiber bending. In consequence, local directional 
measures which describe the orientation of fibers in SMC components 
are approximately planar. In order to elaborate the consequences of this 
planarity, directional measures are introduced, closely following a more 
comprehensive discussion in Bauer and Bohlke (2022c, sections 2.1 and 
2.2). 

A fiber orientation density function (FODF) is an established directional 
measure for the orientation of axisymmetric fibers contained inside a 
reference volume. Such a reference volume can be interpreted as a 
section of a structural component of specified size. The size of the 
reference volume influences the directional measurement and represents 
a scaling parameter in measurement algorithms, e.g., in Górthofer et al. 
(2019, Figure 4) or Schóttl et al. (2020). The FODF 


v:S? >R, with S? = {n € R? ||nl|| = 1] (6.1) 
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at a given position inside a component maps any direction n, being 
part of the two-dimensional surface of a unit sphere in three dimen- 
sions S?, onto a scalar value 7 (n). The function v/ (n) is non-negative, 
normalized and symmetric (see Advani and Tucker III (1987); Bauer 
and Böhlke (20226); Górthofer et al. (2020); Bauer and Böhlke (2022a)). 
The FODF contains the complete directional information, which can be 
contained in a one-point correlation function, but is usually unknown. In 
contrast, orientation tensors only contain a limited amount of averaged 
information but can be obtained by computer tomography scans or flow 
simulations (Górthofer et al., 2019) and fit into the tensor framework of 
continuum mechanics. In this work, fiber orientation tensors of Kanatani 
first kind (Kanatani, 1984) of second N and fourth order N are used and 
defined by 


N- v (n) n $9 ndn, (6.2) 
S2 


N= w(n)n@n@n@ndn. (6.3) 
S2 


The fiber orientation tensors represent averaged directional measures 
and are coefficients of a three-dimensional tensorial Fourier series repre- 
sentation of the FODF 


1 «€ 2k-4-1/2k 
=0 


called spherical harmonic expansion (Kanatani, 1984, page 154), with the 
operator dev(-) extracting the deviatoric part, see Spencer (1970). Higher 
order fiber orientation tensors contain all tensors of lower order which 
implies N = N [I] for the tensors orders two and four, with the identity 
on second-order tensors I = 6;; e; Q ej. 

In order to define planarity of the directional measures, the basic prop- 
erties of the second-order orientation tensor N are briefly summarized 
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following Bauer and Bohlke (2022c). As N is symmetric and positive 
semi-definite, it can be diagonalized, i.e., pairs of eigenvalues A; with 
A; > 0 and orthonormal eigenvectors v; for i € [1,2, 3] exist, such that 


N= NE Je, @e; = Y^ vi Q vi (6.5) 
à 0 0 
= Aa 0 Vi e Vj 
sym Ag 


holds with the ordering convention A3 < Aa < A; and there exists a 
rotation defined by an orthogonal tensor 


Q = Vi & e; (6.6) 


mapping the arbitrary but fixed basis {e;} onto the basis {v;}. A visual 
classification of structurally differing orientation tensors is given by the 
so called orientation triangle (Cintra Jr and Tucker III, 1995; Chung and 
Kwon, 2002; Goldberg et al., 2017; Kóbler et al., 2018). As A3 = 1— Aı — Aa 
holds due to the normalization of 7, any second-order orientation tensor 
can be represented by a pair (A1, A»), which is connected to a point 
inside the orientation triangle and by a mapping Q which defines the 
orthonormal basis {v;} spanned by the eigenvectors. This basis is called 
orientation coordinate system (Bauer and Bóhlke, 2022c) and can be used 
to define the term planarity for FODF and fiber orientation tensors. A 
FODF is planar, if 


v(n(o, 9)) =0 Vp, VOS T (6.7) 


holds with a unit vector parameterized in two spherical angles in the 
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1 
2 
1 A 
3 1/8! 
"d 
6 d, 0 
LAE T 
1 4 5 1 — 6/6 
6 
p = -1/81 Fi 
1 
Aı 0 1/8 3/6 
@ isotropic @ unidirectional ds 


© planar isotropic e SMC CT scan data (b) Visualization of the set of admissible 


and distinct planar fourth-order fiber 
orientation tensors VP!" in the 
parameter space D» di, ds } using the 


(a) Orientation triangle visualizing second- 
order orientation tensors NS obtained 
by computer tomography (CT) scans at 


several positions of a SMC specimen. For parametrization of NPlanar in Equation 


reference, the vertices of the orientation (6.11). A derivation of this space is given in 
triangle are labeled in the legend. Bauer and Bóhlke (2022a) 
Figure 6.1 


orientation coordinate system 
n(y, 0) = sin(0) cos(p)vı + sin(0) sin(p)va + cos(0)va. (6.8) 


A fiber orientation tensor of arbitrary order is planar if its successive 
contraction yields a second-order fiber orientation tensor with vanishing 
third eigenvalue. 

The assumption of the planarity of SMC microstructures is assessed by 
inspection of a biaxial tensile specimen made of SMC. Second-order fiber 
orientation tensors following equation (6.5) at several locations inside 
this specimen are obtained by computer tomography following (Schóttl 
et al., 2020) and are referenced as NGjc . The obtained fiber orientation 
tensors are visualized as points inside an orientation triangle in Figure 
6.1a. The orientation states along the upper right border of the triangle 
in Figure 6.1a represent perfectly planar orientations, as the third eigen- 
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value A3 = 1 — A1 — A» vanishes. The measured fiber orientation tensors 
align well onto the planar border of the orientation triangle in Figure 6.1a 
and therefore confirm the planarity assumption. However, the spatial 
resolution of CT images as well as the choice of algorithm parameters for 
the processing of the CT images, influence the planarity of the resulting 
fiber orientation tensors. Not all fiber orientation tensors obtained for 
SMC in the literature are perfectly planar. As an example, Kehrer et al. 
(2020, Figure 5) contains non-planar fiber orientation states, obtained by 
CT scans, indicating a high spatial resolution of the corresponding CT 
images combined with small extraction cylinders (see Górthofer et al. 
(2019, Figure 4)) in the post-processing of the CT images. Visualizations 
of SMC microstructure data obtained by computer tomography (CT) 
are, e.g., given in Trauth et al. (2021, Figure 2). Artificially generated 
SMC microstructures can be found in Górthofer et al. (2020, Figure 3b 
and 5). The algorithm, which Górthofer et al. (2020) used to generate 
artificial microstructures, hints towards the planarity of the resulting 
directional measures. Throughout this work, fiber orientation tensors 
of SMC are assumed to be perfectly planar and hence, the dependence 
of the mechanical response of shell-like SMC structures on a varying 
microstructure can be investigated based on planar fiber orientation 
tensors. In this work, the SMC-specific bundle structure which is, e.g., 
considered in Schóttl et al. (2021); Górthofer et al. (2020), is not taken 
into account. 
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6.3 Planar fourth-order fiber orientation ten- 
sors 


A parameterization of planar fiber orientation tensors of second-order 
Nplanar is given by 


NPlanar (A4) = Ai v4 Q vi + (1 — A1) V2 Q v2 (6.9) 
Ài 0 0 
= 1-Xı 0 | ;®v,; 
sym 0 


in the orientation coordinate system {v;} which is defined in Equation 
(6.6). Positive definiteness and normalization of the trace combined 
with the ordering convention of the eigenvalues demand the parameter 
range 1/2 < A, € 1. Bauer and Böhlke (2022c) discuss that planar 
fiber orientation tensors of fourth order only depend on three inde- 
pendent parameters and derive admissible ranges of these parameters 
demanding positive definiteness of the orientation tensors. Among the 
admissible parameters, Bauer and Bóhlke (2022a) identify a subset of 
parameter combinations which represent all admissible and structurally 
distinct planar fourth-order fiber orientation tensors. Following Bauer 
and Böhlke (20222), two tensors A and B are structurally distinct if 


AQ € SO(3) with Q« A = (6.10) 
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with the special orthogonal group SO(3), which contains all proper 
rotations. These results are summarized by the parameterization 


Nplanar (A5. 8) e (6.11) 
sin (B )- R(A1) +Aı f sin (8 + Ri (A1) 0/0 0 V2rcos (8 
—f sin 8) - R(u) + (1— X1) 0|0 0 —V2fcos 3) 
0/0 0 0 
completely | symmetric 


Bz & BX 
which follows from Bauer and Böhlke (2022a, Equation (22)) with 
R(A1) = (Aı - 2/2 (6.12) 


and is represented in Kelvin-Mandel notation (Thomson, 1856; Mandel, 
1965; Mehrabadi and Cowin, 1990). This notation is explained in detail 
in Appendix B.1 and the basis B? & B? is spanned in the orientation 
coordinate system {v;}, i.e., for example the fifth basis vector in Equation 
(6.11) becomes BY = v2 (v1 Q v3 + v3 Q vı). The parameterization in 
Equation (6.11) can be combined with the set of parameter combinations 
which leads to admissible and distinct NP!ara, This set is 


A a 1 A 
planar _ [grim (^. f, 8) |5<MS1L0<F<S RM), (613) 


U {pane (^. f, 8) E 


following Bauer and Bóhlke (2022a, Equation (28)). A visualization of the 
body of admissible and distinct orientation tensors NP!" in Cartesian 
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coordinates 
d, = sin (8) j dg = f cos (8) (6.14) 


is given in Figure 6.1b. The data set VSI. of fiber orientation tensors 
obtained by CT, which second-order information is visualized in Figure 
6.1a, is represented in the parameter space VP!"*" in Figure 6.2. The 
data points are preferably concentrated around the axis of vanishing 
values of d; and dg. This observation supports the statement of Hine 
et al. (2004), that real fiber orientations preferably scatter around orienta- 
tion states with maximum entropy, i.e., those which are approximately 
isotropic. Analysis of additional data sets reveals, that the size of the 
reference volume associated with each discrete fiber orientation tensor 
significantly determines the resulting distribution inside the admissible 
region. However, quantitative investigations on the dependence of fiber 
orientation states on processing and CT settings is moved to a follow-up 
publication. Based on the parameterization Nplanar( X. f. B ) and the set 
AfPlarer a visualization method, motivated by Bauer and Bóhlke (2022a), 
can be used to study the influence of planar fourth-order fiber orientation 
tensors on the mechanical properties predicted by orientation-averaging 
mean field homogenization techniques. 

Although the planar fiber orientation tensors given in Equations (6.9) 
and (6.11) fit into a two-dimensional framework, e.g., used in Bauer 
and Bóhlke (2022a), a three-dimensional tensor framework is deployed 
within this work, as the mechanical behavior of real materials is defined 
by the laws of physics of the three-dimensional world. Boundary condi- 
tion into the out-of-plane direction influence the mechanical response of 
planar structures, see Nordmann et al. (2020). 
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6.4 Orientation averages 


Orientation-averaging mean field homogenization demands analytical 
or numerical schemes which yield averages of direction-dependent 
tensorial quantities over orientations. Advani and Tucker III (1987) 
define the orientation average of a tensorial quantity A by 


(A) ar = B Y(n) A(n) dS, (6.15) 
i.e., as a weighted summation with weights defined by the FODF. Sev- 


eral formulations and approxmiations of Equation (6.15) exist in the 
literature. An exact formulation given in the original paper Advani 


0 —— 3/6 


1/8 
dg / 


Figure 6.2: Fourth-order orientation tensors ASI obtained by computer tomography 
(CT) scans at several positions of a SMC specimen, represented by points inside the 
admissible and distinct planar parameter space VP!*"*", The size of the reference volume 


associated with each discrete fiber orientation tensor significantly determines the resulting 
distribution inside the admissible region. 
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and Tucker III (1987, Equation (28)), frequently used in, e.g., Lielens 
et al. (1998); Jack and Smith (2008); Camacho et al. (1990); Brylka (2017); 
Kehrer et al. (2020); Hessman et al. (2021) is labeled (-) rn in this 
work. This formulation is given directly in terms of fiber orientation 
tensors of second and fourth-order Kanatani first kind (Kanatani, 1984) 
, and in its original formulation is restricted to transversely isotropic 
elasticity tensors H = H;jrı ei 8 e; ® ej; & ej, i.e., stiffness or compliance- 
like quantities with major and both minor symmetries which may be 
specified in index notation by H;jj = Hji = Hiju. However, the 
original formulation in Advani and Tucker III (1987), which is repeated in 
slightly varying notation in, e.g., Kehrer et al. (2020); Brylka (2017); Heller 
et al. (2016), can be extended to transversely isotropic tensors which 
lack the majory symmetry of an elasticity tensor. Such a generalization 
is given in Appendix C.4 based on the formulation introduced in the 
next section. The explicit formulation of (- ) ary in orientation tensors is 
simple and efficient, but complicates adaption of empirical fiber dam- 
age (Schemmann et al., 2018b) or incorporation of direction-dependent 
fiber length distributions (Brylka, 2017). A reformulation based on the 
harmonic decomposition Rychlewski (2000); Forte and Vianello (1996); 
Olive et al. (2018); Bóhlke and Brüggemann (2001) explicitly revealing 
the structure of the scheme is given in the next section. 


6.4.1 Reformulation of the explicit Advani-Tucker orien- 
tation average 


A transversely isotropic elasticity tensor is completely defined by five 
independent parameters and a direction, i.e., a normal vector q, and can 
be written as C™"SY (hy, ha, ha, ha, hs, q) following Equation (C.3) in 
C.2. The orientation average (- ) ary of a transversely isotropic elasticity 
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tensor is independent of q and given by 


(coe (hi, ha, ha, ha, hs, q))ATN (N) (6.16) 
= hı Pı + hg Po + hg J3 [dev (IN)] + ha J4 [dev (N)] + hs dev(N) 


with the irreducible parts of N 


dev(N) = N— N*? (6.17) 


dev(N) = N — sym (dev(N) & I) — N^ (6.18) 


and the isotropic parts of the second- and fourth-order fiber orientation 
tensors of Kanatani first kind (Kanatani, 1984) 

Ni» — li N = n (IG I) (6.19) 
37 35 
following Bauer and Bóhlke (2022c). The isotropic sixth order tensors 
Js and J; are defined in Equations (C.6) and (C.7). The operator sym(-) 
extracts the totally symmetric part with respect to symmetry following 
Spencer (1970) and the operator dev(-) extracts the deviatoric part also 
discussed in Spencer (1970). Linearity of the harmonic decomposition in 
Equation (6.16) implies that the orientation average following Advani 
and Tucker III (1987) is linear in both C™™sY and N. This implication 
might be obvious to some authors, as the only source of orientation 
dependent structural information is N. However, the original formu- 
lation in Advani and Tucker III (1987) does not directly reveal, that 
the isotropic part of the averaged elasticity tensor is unaffected, as one 
would expect from a physical point of view. In addition, the first and 
second deviatoric parts of the orientation average in Equation (6.16), i.e., 
hg dev(N) and h4 dev (N), only differ by a scaling factor and share the 
tensor structure of dev(N). As N is completely symmetric, its first and 
second deviatoric parts coincide, see Bauer and Bóhlke (2022c). As the 
manual derivation of Equation (6.16) is lengthy, listing C.1 combined 
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with code in Bauer (2022) is used to validate the representation of the 
orientation average. Any material model which contains the Advani- 
Tucker orientation average, implicitly contains a linearity assumption 
of at least an intermediate quantity in the fiber orientation or fabric 
tensor. For reference, some authors Schemmann et al. (2018a;b); Karl 
et al. (2021b) explicitly postulate linearity of an effective stiffness or 
viscosity in fiber orientation or fabric tensors. In contrast, a material 
model of a porous, anisotropic, linear elastic material being non-linear 
in a fabric tensor is given in Cowin (1985). 


6.4.2 Direct numerical integration and the adaptive 
scheme based on angular central Gaussian distribu- 
tions 


Direct numerical integration of Equation (6.15) requires both a given 
FODF and a large number of integration points on the two-dimensional 
area of integration 5?, but leads to insufficient accuracy despite high 
numerical effort, especially in the case of strongly localized FODF. This 
performance and accuracy issue motivates a recently developed nu- 
merical scheme (Goldberg et al., 2017) denoted by (-) argos which is 
based on a special class of FODF basis functions called angular central 
Gaussian. Although this scheme is not used within this work, it is 
shortly introduced as a comparison. The scheme (-) \rcos leads to an 
approximation resulting in a weighted summation of the quantity of 
interest pointing into a number of N©° discrete directions 


NGOS 


(A) av © (A) araos (N) = 3 wy? (N) A(n;) (6.20) 


with the weights being a function of the second-order fiber orientation 
tensor. The weights can be pre-calculated and stored in efficient look-up 
structures. Implementation details are given for the three-dimensional 
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case in Goldberg et al. (2017); Hessman et al. (2021) and for the two- 
dimensional case in Górthofer et al. (2020). This formulation yields 
reasonable accuracy even for localized FODF with a limited number of 
integration points (Goldberg et al., 2017), e.g., N GOS — 434 in Hessman 
et al. (2021). Due to the structure in Equation (6.20), any direction- 
dependent quantity can be averaged and the weights may be modeled 
for empirical simulation of damage. However, the orientation average 
(-) argos is solely based on the directional information contained within 
a second-order orientation tensor, which is limited to two scalars de- 
scribing the half axes of an ellipsoid combined with the definition of a 
coordinate system. 


6.4.3 Orientation average by reconstructed planar FODF 
based on a maximum entropy method 


A maximum entropy FODF reconstruction for planar fiber orientations 
developed by Bauer and Bohlke (2022a) based on Miiller and Bohlke 
(2016), can be combined with any numerical integration on S' leading 
to an approximation of the Advani-Tucker orientation average, which 
shares the simple structure of (A) argos, uses fourth-order fiber orien- 
tation tensor information and makes use of the reduced dimension of 
planar fiber orientation tensors. However, as this scheme uses direct 
numerical integration on S ! it leads to insufficient accuracy for localized 
fiber orientation states. Fortunately, the localized orientation states are 
known through analysis of the admissible region \’?*"*" in Bauer and 
Bóhlke (2022a) and can be captured in implementations. Following 
Bauer and Böhlke (20222), for any admissible non-localized fiber orienta- 
tion tensor of fourth-order, a planar FODF, m (p, NPlanar) parameterized 
in one spherical angle defined in Equation (6.8) with 0 < y < 27 can be 
identified. This leads to an approximation of the Advani-Tucker average 
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in Equation (6.15) for the planar case by 


(A) ar (NPA) = (A) (NP) 


2n 


E! d. (p, NPP) QC) x A dip 


NME 


„ME 
zn) wih (pi, NP) Q(yi) x A 
w=1 


NME 
= 27), wi (NP) Q(yi) x A (6.21) 


i=l 


with a rotation around the vs axis Q, evaluated at specified angles which 
are defined by any numerical integration scheme on 5! with weights 
w;. Within this work, numerical integration on S! with homogeneous 
weights w; = 1/NVF and equidistant angles w; = 27 i/ NV following 
Krylov and Stroud (2006) is used. The weights of the numerical integra- 
tion scheme might be combined with the direction specific value of the 
FODF in Equation (6.21). The accuracy of the average (-),;, is assessed 
by averaging of moment tensors, as 


(nSt) ar (NP!anar) — Nplanar x ae") up (N) (6.22) 


has to hold. The assessment is done for all fourth-order fiber orientation 
tensors specified in Table C.1 except for the unidirectional case notated 
as NER ‚as this localized orientation state is insufficiently handled by 
the numerical integration. The deviation is quantified by 


@4 NME Nplanar = Nplanar 
AR AIME = max | (n Imel ’ ) | 


NPlanar œ AyPlanar | | Nplanar | | 


selected 


(6.23) 


and given for varying order of the numerical integration NMÉ in Figure 
6.3. The target precision of 1077 is reached. Studying the convergence 
of the numerical averaging scheme, by averaging a specific, randomly 
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Maximum error among N in N! plar A 
1071 Target precision | 
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Figure 6.3: Influence of the number of integration points on the deviation of (n®*) y 
from (n84) yry with A™* defined in Equation (6.23) 


chosen, elasticity tensor instead of the moment tensor in Equation (6.23), 
following Hessman et al. (2021), leads to similar convergence. 


6.5 Orientation-averaging mean field homoge- 
nization 


Sheet molding compound is a two-phase composite with glass fibers 
embedded irregularly in a thermoset polymer matrix. For isothermal 
applications, isotropic linear elasticity states a reasonable assumption 
for both the fiber and the matrix behavior, see Kehrer et al. (2020). In 
consequence, the local stress-strain relation inside the fiber and matrix 


phase at position x is 


C; [e(x)] Vx € Bi (6.24) 
Cm [e(x)] Vx € Bm (6.25) 


a 
x, 
J 


a 
E 
J 
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Young's mod. | Poissons ratio || Bulk mod. | Shear mod. 


Fibers 73.0 GPa 0.22 43.45 GPa | 29.92 GPa 
Matrix 3.4 GPa 0.385 4.93 GPa 1.23 GPa 


Table 6.1: Material parameters of the SMC constituents glass fiber and neat matrix (UPPH, 
i.e., unsaturated polyester-polyurethane hybrid) adopted from Kehrer et al. (2020, Table 2) 
citing additional references. Bulk and shear modulus are derived from Young's modulus 
and Poisson's ratio. The volume fractions of the fiber and matrix phases are cy — 0.256 
and cp — 0.744. 


with parts of the representative volume element (RVE) occupied by 
the fibers B; and the matrix Bm, respectively. The stiffness of a fiber is 
denoted by C; and Cm represents the stiffness of the matrix. Isotropy 
implies that C; and Cm in Equations (6.24) and (6.25) do not depend on 
the spatial orientation and are given by 


C; = 3K; Pı + 2G; P» (6.26) 
Cn = 3K, Py + 2G, P (6.27) 


with the first and second isotropic projectors Pı and P» respectively 
and bulk and shear modulus denoted by K and G. Throughout this 
work, material parameters are adopted from Kehrer et al. (2020, Table 
2) and listed in Table 6.1. Two one-point characteristics describing the 
microstructure of a RVE at a material point inside a SMC-component are 
considered. The first characteristic is the volume fraction of the fibers cz, 
which implies the volume fraction of the matrix cm = 1 — cr. The second 
characteristic is a fourth-order fiber orientation tensor N. In the absence 
of any non-linearities, such as cracks, the effective elastic stiffness C of 
the RVE exists and is of interest. The effective stiffness maps effective 
strains € = (e(x)) onto effective stresses « = (o(x)) by 


g = Cl[g]. (6.28) 


The operator (-) takes the volume average of a spatially varying field 
quantity, e.g., a (x), inside the RVE. Four mean field approximations 
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are briefly reviewed and examined for given material parameters and 
variable fiber orientation tensors. 


6.5.1 Two-step Hashin-Shtrikman 


A two-step Hashin-Shtrikman scheme for fiber reinforced composites 
following Kehrer et al. (2020), based on Walpole (1966a;b; 1969) and 
proposed by Fernández and Bóhlke (2019) is used to study the influence 
of a fourth-order orientation tensor on the effective stiffness. Detailed 
derivations of the Hashin-Shtrikman bounds, based on a variational 
principle, are given in Walpole (1966a); Fernández and Bóhlke (2019). 
In this work only the resulting equations are of interest. Under the 
assumption of phase-wise constant stress polarizations, the Green’s 
function of a material without long-range order and ellipsoidal two- 
point statistics, is constant (Fernandez and Bohlke, 2019, section 2.2). 
This constant Green's function is represented by the Hill's (Hill, 1965a) 
polarization tensor P and reflects the ellipsoidal symmetry of the two- 
point statistics. Under these assumptions, the effective stiffness of a 
material made up of n phases with stiffnesses (C;) for i € [1,..., n] is 
given as a function of a comparison stiffness Co by 


CHSW (Co, P (Co, a), {Ci}, ((-))) = Co - P^! + [((W (Co, P, C)))] * 
(6.29) 


with 
W (Co, P, Cj) = [P-! + C; - Co] ^ (6.30) 


and with P (Co, a) depending on the symmetry of the two point statistics. 
For simplicity, the dependence of P upon the two point statistic is 
represented by an aspect ratio a with 1/2 < a < oo reaching from the 
isotropic case P? (Co) = IP (Co, 1/2) to the limiting unidirectional case 
PUP (Co) = P (Co, oc). This restriction of P from generally ellipsoidal 
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shapes to spheroidal shapes is common in the context of fiber reinforced 
composites. The operator ((-)) takes the average over the RVE. Al- 
ternative formulations of the Hashin-Shtrikman scheme exist, e.g., in 
Willis (1977; 1981) and remarks on the different formulations are given in 
Fernández and Bóhlke (2019). The formulation in Equation (6.29) is not 
directly applicable to singular, i.e., non-invertible, polarization tensors. 
However, Walpole (1969, page 238) proposes a dual scheme which can 
be used to evaluate Equation (6.29) for singular polarization tensors, e.g., 
PVP, If the comparison material with stiffness Co is selected, such that the 
stiffnesses of all phases are smaller (larger) than Co, e.g., using first order 
bounds, Equation (6.29) yields a lower (upper) Hashin-Shtrikman bound. 
As SMC is a two-phase composite, for any choice of the comparison 
stiffness Co in between Cm and C;, Equation (6.29) gives an admissible 
effective stiffness in between the Hashin-Shtrikman bounds. Following 
Kehrer et al. (2020, Equation (24)), a comparison stiffness as a function 
of an interpolation parameter k is introduced by 


Ĉo (k) = [1 — k] Cm + k Ci (6.31) 


with 0 < k < 1. With this notation, the scheme of Kehrer et al. (2020) 
is reformulated as a generic two-step Hashin-Shtrikman scheme. In a 
first step, the effective stiffness of an artificially pseudo grain CUP is 
calculated by 


CUP (ki, eg, Cg, Cm) = CHSW (Co = Ĉo (k1), P = PYP (Co), (632) 
{Ci} = (Cs Cm}, (00) 7 (02). 


This pseudo grain contains unidirectional fibers with a volume con- 
tent which is equal to the volume content of the fibers in the overall 
two-phase composite cy. To be explicit, the operator (-) in Equation 
(6.32) applied to W (Co, P, C;) with i € [m,f] reads as (W(..., C;)) = 
Cm W (..., Cm) + c¢W (..., Ce), with the abbreviation "^..." for additional 
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function arguments. The stiffness CUP obtained by Equation (6.32) is 
transversely isotropic. As the Hill polarization P? in Equations (6.30) and 
(6.32) reflects the shape of the two-point statistics, the specific choice 
of PUP for a single grain implies a laminate-like two-point statistic. It 
should be noted, that the inclusion shape is not specified. However, 
specification of the two-point statistic is not independent of the inclusion 
shape. In the second step, the stiffness obtained in step one is orientation 
averaged based on the Hashin-Shtrikman scheme in Equation (6.29), 
where the average over the RVE ( (-)) is transferred into an average over 
orientations. The operator (-) ‚rn defined in Equation (6.16) represents 
a commonly used Lielens et al. (1998); Jack and Smith (2008); Camacho 
et al. (1990) average over orientations, based on a fourth-order fiber ori- 
entation tensor N, following Advani and Tucker III (1987). As the pseudo 
grain is artificial, the two-point statistics in the second step are assumed 
to be isotropic and are therefore reflected by the Hill polarization tensor 
for spherical inclusions P°. This leads to the approximation C1? of the 
effective overall stiffness with 


com (ki, ka, N, Cf, Ch Cm) = (6.33) 
CHW (Co = Co (ka), P = P^ (Co) (C) = {EP (k, .)). 


(0) = (ann N) ). 


The connection to the formulas in Kehrer et al. (2020) is given in C.5. 
The formulation of (-) yry in Equation (6.16) shows the linearity of the 
orientation average in both arguments. 

For the special case of planar fourth-order orientation information, i.e., 
Nplanar and a given transversely isotropic pseudo grain stiffness CUP, the 
structure of the effective stiffness is given analytically. Starting from a 
given transversely isotropic stiffness CU obtained in a first homogeniza- 
tion step, e.g., following Equation (6.32), the second homogenization 
step in Equation (6.33) for the special case of a planar orientation tensor 


150 


6.5 Orientation-averaging mean field homogenization 


Nplanar, is of interest. Omitting functional dependencies, the expansion 


of Equation (6.29) reads as 
CHSW = Co — [p] +G 


with 


—1 


o = [P] +0- co] uns] 


=[A Yarn] =[(B)arn] = [E]. 


The Hill polarization P? is isotropic and its inverse is given by 
[P^] ^ = ATP, + APP; 


with 


SAS [h + 2h2] 


PO- 0 0 PO- _ 
Au cef ag, ha = MS a Bn 


(6.34) 


(6.35) 


(6.36) 


(6.37) 


and ^? = 3K?, h$ = 2G? and the bulk and shear modulus of the isotropic 
comparison material K? and G?. As both the Hill polarization and the 


stiffness of the comparison material are isotropic, the intermediate quan- 


tity A in Equation (6.35) inherits the transversely isotropic symmetry 
from the unidirectional stiffness CUP. With the short hand notation of 


a transversely isotropic elasticity tensor defined in Equation (C.13) the 


tensor A is given by 


A = [his has has ha; hs q = vi | 


(6.38) 
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with 
hy = hk — h? RATS, ha = hb -h8 + AP, (6.39) 
lig es hk”, ha = hP, (6.40) 
hs = hyP (6.41) 
based on 
CP = |nYP, AYP, AYP, AYP, AYP, — v; |. (6.42) 


As the orientation of the pseudo grain stiffness is arbitrary, without loss 
of generality, q = vı holds in Equations (6.38) and (6.42). The inver- 
sion of a transversely isotropic stiffness is discussed in, e.g., Lubarda 
and Chen (2008, Equation (31)), and leads to a transversely isotropic 
compliance 


=A= |t, h, hB, hë, hE, q = vı] (6.43) 
By, By Biz 0 0 0 
Bə2 Bo 0 0 0 
B 0 0 0 
= 22 B. @ Be 
Boz — B» 0 0 
sym 2 Dss 0 
2B55 


with A? for i € (1, 2, 3, 4, 5} given in C.6 in Equations (C.33) to (C.38). 
The general correspondence of stiffness and compliance is discussed in 
Rychlewski (1984). Simpler representations of a transversely isotropic 
compliance (Lubarda and Chen, 2008; Cowin and Van Buskirk, 1986; 
Vannucci, 2018), e.g., in tensor components or engineering notation, 
exist. However, Equations (C.33) to (C.38) show the interaction of 
the harmonic coefficients h; with i € (1, 2, 3, 4, 5} during the tensor 


inversion. Each coefficient of B is a function of all five coefficients of A. 
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Taking the orientation average (Advani and Tucker III, 1987) of B, e.g., 
using Equation (6.16), leads to 


Ei En Es 0 0 Eis 
E» Es 0 0  —£ig 
i = (B) apy = 28 » : ° Be & Be. (6.44) 
44 
sym Ess 0 
Ess 


The tensor components in Equation (6.44) are specified in Equations 
(C.41) to (C.50) in C.6. The inversion of E leads to 


Gi G1» Gi3 0 0 Gis 
G22 G23 0 0 G26 
G = [E]! = Ga O O Ga Be&B. (6.45) 
G44 0 
sym Gs5 0 
Gee 


with tensor components given in Equations (C.51) to (C.63) in C.6. 
Equations (C.51) to (C.63) combined with Equations (C.41) to (C.50) 
show that the effective stiffness CĦSW is nonlinear in the parameters 
of the fiber orientation tensor A, di, dg. In addition, the structure of 
CHSW and G differs from that of the intermediate quantity E, as the 


tensor component G36 in Equation (6.45) does not vanish, although E36 
in Equation (6.44) does vanish. CHSW inherits the structure from G. 


6.5.2 Orientation-averaging Mori-Tanaka following Ben- 
veniste (1987) 


Various aspects of the orientation-averaging Mori-Tanaka approximation 
are discussed in Brylka (2017); Weng (1990); Qiu and Weng (1990). 
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Nevertheless, a comprehensive summary of the basic equations for the 
special case of homogeneous fiber lengths, i.e., an isotropic fiber length 
distribution, is given hereafter. The RVE is denoted by B and any point 
inside the RVE is identifiable by a position x € B. Introducing a field of 
strain localization tensors A(x) mapping the effective strain of the RVE, 
€, onto the local strain e(x) with 


e(x) = A(x)[é] (6.46) 


C = (C(x)A(x)) (6.47) 


due to 


= (C(x)A(x)[ē]) = (C(X)A(x) ) [e]. (6.48) 


If the exact strain localization field A(x) for an RVE is known, the 
effective stiffness C is given by Equation (6.47). The volume average of 
Equation (6.47) implies (A(x) ) = I? with I? being the identity acting on 
the space of symmetric second-order tensors. An exact decomposition 
of the volume average over the RVE into volume averages over the 
different phases in the special case of a two phase composite yields 


(A) 2 P = cm (A), + cr (A), (6.49) 
with the volume fraction c, of phase a € [f, m] and therefore 


(A), = — [P — c4]. (6.50) 


m Cm 
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For a two-phase material with piece-wise constant stiffnesses, Equation 
(6.47) combined with Equation (6.50) leads to 


€ = om (CG)A G9), + CAR), (651) 
= Cm Con (A (x) ),, + er Cr (A (x); (6.52) 
= Cn [I — (A), + CHAR), (6.53) 
= Cm +c AC (A), (6.54) 


introducing the short hand notation for the material contrast AC = 
C;— Cm. The strain localization tensor of the Mori-Tanaka approximation 
(Mori and Tanaka, 1973) for a two-phase composite reads as 


(A), © (A)MT = A [es IP eg AS!] 
= [em [Af] +e] (6.55) 


with the strain localization tensor of the inclusion (fiber) phase in the 
single inclusion problem Atl (Eshelby, 1957) given by 


AF = [P (Cm, a) [Ce — Cm] + IB] 5 . (6.56) 


The Hill polarization tensor P(Cm, a) for spheroidal inclusion shapes is 
parameterized by an aspect ratio a. Inserting the approximated strain 
localization tensor of Equation (6.55) into Equation (6.54) leads to the 
effective stiffness of the Mori-Tanaka approximation for a two-phase 
composite 


CMT (cr, Ce, Cm, P) = Cm to AC (AMT, (6.57) 
f f 


The Hill polarization P in the single inclusion problem reflects the shape 
of the inclusion. A study on the influence of the inclusion shape onto the 
effective stiffness is given in Müller (2016). The large aspect ratio of SMC 
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fibers motivates the use of PU, e.g., in Kehrer et al. (2020). Orientation- 
averaging Mori-Tanaka for two-phase composites with anisotropic con- 
stituents and an inclusion phase consisting of aligned or randomly 
orientated ellipsoidal particles following Benveniste (1987) combined 
with an orientation-averaging scheme following Advani and Tucker III 
(1987), depicted by (-) yry, reads as 


(A), & (Ale PN, cr Ce, Cm) (6.58) 
—1 


= [em [ACI CAB ame] arl 


leading to the Mori-Tanaka orientation averaged approximation follow- 
ing Benveniste (1987) with 


CMTOAB (N, cj, Cr, Cm) = Cm + cr AC (A) ^ (N, cr, Cr, Cm). (6.59) 


Orientation-averaging Mori-Tanaka following Benveniste (1987) is used 
in Brylka (2017); Schemmann et al. (2018b) and limitations are discussed 
in, e.g., Benveniste et al. (1991). As the orientation-averaging scheme 
(Jary in its original formulation is limited to elasticity tensors, i.e., 
any tensor C = Cj;xi ei Q e; & ej ® ej with both minor and the major 


symmetries Cijk = Cjiy; = Cijtk = Cklij, but strain localization tensors 
may lack the main symmetry, an intermediate elasticity tensor AC A?! = 
[P + AC 1] ~* is averaged in Equation (6.58) before the effect of AC is 
removed by AC~!. This step could be avoided using the generalization 
of the orientation average (-) pm in Appendix C.4. The unity of all 
fibers is treated as one phase. 
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6.5.3 Direct orientation average of a transversely isotropic 
stiffness 


Some authors (lorga et al., 2008; Schjedt-Thomsen and Pyrz, 2001) iden- 
tify the linear orientation average of transversely isotropic stiffness 
(Advani and Tucker III, 1987) as an approximation for the effective 
stiffness with 


CM TlinearStiffness (N, ...) = cm (av (N) (6.60) 


which is based on the Mori-Tanaka approximation in Equation (6.57) 
and is used in commercial software (Smith, 2019). 


6.5.4 Direct orientation average of a transversely isotropic 
compliance 


In analogy to the approximation defined in Equation (6.60), another 
approximation is given by the orientation average of an unidirectional 
compliance obtained by the Mori-Tanaka approximation 


(@MTlinearCompliance (N, ...) = | (jem (...)] —1 lows (N)| l (6.61) 


6.6 Graphical representation of elasticity ten- 
sors 


The last building block for the systematic investigation of the depen- 
dence of linear elastic effective stiffness tensors on planar fiber orienta- 
tion tensors, is a compact, preferably two-dimensional, visualization of 
effective stiffness tensors. Following Böhlke and Brüggemann (2001), 
two scalar functions E (C, n) and K (C, n), equivalently describe an 
anisotropic stiffness C. The direction-dependent Young’s modulus 
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E (C, n) defined by 


iG gy e = fete] (6.62) 


n) 


_ BE ‘Clio ep = [c*. ne]! 


c (n) 


q 


represents the ratio of the tensile stress c (n) and the tensile strain e (n) 
during a tensile test in direction n. The tensile stress caused by a virtual 
unidirectional tensile test into direction n, is related to the stress tensor 
o (n) by e (n) = o n?? = en & n. The tensile strain e (n) is obtained by 
a projection of the infinitesimal strain tensor e onto the direction n, i.e., 
e = € n??, Following He and Curnier (1995); Böhlke and Brüggemann 
(2001), the generalized bulk modulus K (C, n) is defined by 


3 K(C,n)- oT = Ba (6.63) 


ce (n) 


= [r.c pe] 


and measures the relative change of volume I: e = [dV — dV] / dvo 
caused by uniaxial tension in direction n, i.e., o = o n??. For a given 
stiffness C, both E (C, n) and K (C, n) are functions on the unit sphere. 
The influence of NPlanar on the mechanical properties in the vi-v2-plane 
is investigated by 


gp" (C, g) = E (C, 
K Planar (C, p) —K (C 


(e, 6 — 1/2) (6.64) 
(e, 6 — 1/2) (6.65) 


n 
,n 
with the unit vector n (v, 0) specified in Equation (6.8). The two quan- 
tities HP!" (C. p) and KPlanar (C, p) uniquely define the mechanical 


response of the effective stiffness in the plane spanned by v; and v». 
Alternative representations of three-dimensional fourth-order tensors 
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Figure 6.4: Ten points in the plane of planar orthotropic fiber orientation tensors of fourth- 
order used in Figures 6.7a and 6.7b 


are, e.g., given in Vannucci (2018, p. 62). 


— CHW? (kı = 0, k2 = 0, N, ...) 
— CHW (gi = 1, ko = 1, N, ...) 
— MTOA (N, ...) 

—— (CMrTlinearStiffness (N " 


0 " 1 /8 3/ 6 = (MTlinearCompliance (N, si) 
(a) Points in the plane Ay = 3/6 used in 
Figures 6.8a, 6.8b (b) Shared legend for Figures 6.7, 6.8, 6.9, 


6.10, 6.11 and 6.12. 


Figure 6.5: Definition of representative points in the parameter space A/P'?'?' in the plane 
A1 = 3/6 and shared legend for mean field approximations. 


159 


6 On the Dependence of Orientation Averaging Mean Field Homogenization 


(y2) 

(g2) 
PAST 
(x) (m1)(m2) 


(oy) 
(c2) 


(02) 


(c) Point labels for 
(a) Points in the plane Aı = 4/6 (b) Points in the plane 4; = 5/6 overview Figures 
used in Figures 6.9, 6.11 used in Figures 6.10, 6.12 6.6a and 6.6b 


Figure 6.6: Definition of representative points in the parameter space N’Plarar of planar 
fourth-order fiber orientation tensors. Parameters of each point are listed in Table C.1. 


6.7 Effective stiffnesses: Polar plots and the 
dependence on planar fourth-order fiber 
orientations 


In this section, the building blocks developed in the previous sections are 
combined to study the dependence of linear elastic effective stiffnesses 
obtained by mean field homogenization on fourth-order fiber orientation 
tensors. As the dimensionality of the study is already high due to the 
flexible orientation, the material parameters of the constituents are fixed 
and specified in Table 6.1. 


6.7.1 Visualization setup 


The space of admissible and distinct planar fourth-order fiber orienta- 
tions N Planar is discretized by a small number of representative fiber 
orientations. Slices through the three-dimensional body \P!*"*" are 
defined in Figures 6.4 to 6.6b. On each slice, a small number of points, 
each representing a specific fourth-order fiber orientation tensor NPlanar, 
are selected. For each point, effective stiffnesses are calculated using the 


160 


6.7 Effective stiffnesses: Polar plots 


E M can 


(a) Young's modulus EPlanar (C (N), e) 


LE 


(b) Generalized bulk modulus «Planar (C (N), e) 


Figure 6.7: Polar plots for selected N specified in Figure 6.4 and Table C.1 and mean field 
approximations specified in the legend in Figure 6.5b 
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orientation-averaging homogenization schemes listed in Figure 6.5b and 
described in Section 6.5. The planar properties of each effective stiffness 
can be visualized by the planar projection of the direction dependent 
Young’s modulus and generalized bulk modulus, introduced in Section 
6.6. In consequence, for each point, two sets of polar plots representing 
Young’s modulus on the one hand and generalized bulk modulus on 
the other hand, are obtained for a selection of homogenization schemes. 
Each set of polar plots is combined into one sub-Figure which is arranged 
according to the position of the point inside the slice of APln?r, This way, 
a graphical representation of the influence of the fiber orientation on 
effective stiffnesses is generated. The first slice, visualized in Figure 6.4, 
contains ten points, each representing one planar orthotropic (Bauer and 
Bóhlke, 2022c, Equation (84)) fiber orientation tensor Nplanar(\,, dı, dg) 
with vanishing parameter dg defined in Equation (6.14). For each point, 


- 


(01) X AN 


- - 


"ea “ʻe” 


(b) Generalized bulk modulus KPlanar (C (N), e) 


Figure 6.8: Polar plots for selected N with Aı = 3/6 specified in Figure 6.5a and Table C.1 
and mean field approximations specified in the legend in Figure 6.5b. The order along the 
path in Figure 6.5a is given from left to right. 
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(02) _ : (2) . 


a=" =.. 


Figure 6.9: Young’s modulus EPlanar (c (N), e) for mean field approximations specified 
in Figure 6.5b. The arrangement of the polar plots follows the arrangement of points in the 
parameter space of Nplanar with A, = 4/6 in Figure 6.6a listed in Table C.1. 
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e udis Mr 7" o- - 


(02) _ P (2) _ 2 


Figure 6.10: Young's modulus EPlanar (C (N), e) for mean field approximations specified 
in Figure 6.5b. The arrangement of the polar plots follows the arrangement of points in the 
parameter space of NPlanat with Aı = 5/6 in Figure 6.6b listed in Table C.1. 
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(02). 


-—--— ~.. 


Figure 6.11: Generalized bulk modulus KP!" (€ (N), e) for mean field approximations 
specified in Figure 6.5b. The arrangement of the polar plots follows the arrangement of 
points in the parameter space of NP'@"@ with Aı = 4/6 in Figure 6.6a listed in Table C.1. 
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(02) _ N (2). P 


Figure 6.12: Generalized bulk modulus K'Planar (C (N), e) for mean field approximations 
specified in Figure 6.5b. The arrangement of the polar plots follows the arrangement of 
points in the parameter space of NP?"?' with A1 = 5/6 in Figure 6.6b listed in Table C.1. 
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a polar plot of Eplanar (C(N), p) is given in Figure 6.7a and a polar plot 
of Kplanar (C(N), y) is given in Figure 6.7b. The arrangements of the 
polar plots in Figures 6.7a and 6.7b mimic the position of the points in 
the overview Figure 6.4 and the correspondence is stressed by the color 
of the thick dashed circles around the polar plots. The line colors of the 
polar plots indicate the mean field approximation which are specified in 
the shared legend in Figure 6.5b. Parameters of each point are given in 
Table C.1 in C. 

Similar visualizations based on slices in Figures 6.5a to 6.6b are given in 
Figures 6.8 to 6.12. Each of the slices in Figures 6.5a to 6.6b represents 
the variety of planar fourth-order fiber orientation tensors for a fixed 
second-order fiber orientation tensor, as the parameter A, is constant. 
The limits of all plots of EPlana are homogeneous and given by 0 GPa 
and 22 GPa. The limits of all plots of KP!*"*" are also homogeneous and 
given by 0 GPa and 24 GPa. 


6.7.2 Observations on bounds 


In the unidirectional case, visualized in polar plot (10) in Figure 6.7a, 
all approximations except CĦSW2(1, 1) coincide (Weng, 1990). Young's 
modulus EP'"?'(C. p) obtained for the approximations CMTOAB and 
CM Tlinear$tiffness are between the values obtained for CHSW?(0,0) and 
CHSW?(1. 1), This does not hold for the approximation CMTlinearCompliance, 
The limiting two-step Hashin-Shtrikman homogenizations are bounds 
on the energy of the effective material. Both Youngs’s modulus and 
the generalized bulk modulus are material characteristics related to 
uniaxial tensile tests. Motivated by Figure 6.7b(10), the complementary 
energy density (Bertram and Gliige, 2013) induced by a virtual direction- 
dependent unidirectional tensile test starting from a stress-free initial 
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configuration 
1 
w*(C,o,n) = 5 c(c,n)-C ! [c(o,n)] = 50? C7! . n84 (6.66) 


is plotted in Figure C.1, C.2 and C.3 in C.8 for a unit stress ø = 1 and for 
all views defined in Figures 6.4 to 6.6b. The energies of the approxima- 
tions CMTOAB and CMTlinearStiffness are within the bounds for all inspected 


(CMrTlinearCompliance 


orientation tensors whereas the approximation violates 


the bounds and therefore is non-physical. 


6.7.3 Observations on the shape of the Young's modulus 


The maximum number of extrema of the Young's modulus EP^^' within 
the v;-v»5-plane is four, as this quantity is obtained by contraction of the 
effective compliance with a fourth-order moment. The maximum values 
of the Young's modulus for all mean field approximations point into the 
directions of maximum fiber content, which are visualized in Bauer and 
Böhlke (2022a). For example, polar plots with Ê = —7/2, i.e., da = 0 and 
dı < 0, have their maxima aligned with the axes v; or va, see Figures 
6.9 and 6.10. 


6.7.4 Observations on the shape of the generalized bulk 
modulus 


In contrast to the Young's modulus, the number of maxima of the 
generalized bulk modulus X planar and the resulting shapes are strongly 
limited, as the highest moment of n which enters KP!" is of second 
order, see Böhlke and Brüggemann (2001). The changes in KPlanar for 
fixed second-order fiber orientation tensor contribution, i.e., fixed values 
of A, but different fourth-order fiber orientation tensor contribution, 
ie, values of d; and dg defining dev(NPlanar), are marginal for the 
selected approximations CHW?(k,, k2) and CMTlinearCompliancee These 
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three approximations take the orientation average on compliance-like 


quantities. In contrast, the influence of both dy and dg on the remain- 


MMTOAB rMTlinearStiffness 
C C 


ing approximations and is clearly visible. The 
latter two approximations take the orientation average on stiffness-like 
quantities. For these approximations, increase of da induces a clock-wise 
rotation of the K?!*"*"-body whereas increase of d, leads to a stretch 
of this body. In Figure 6.8b, no stretch is induced by increase of dı. 
It is noted, that EPlanar and KP!" are obtained by contraction of the 
effective compliance. The strong influence of the fourth-order fiber 


orientation contribution dev(NP!"*") on the shape of the bulk modulus 


MMTOAB rMTlinearStiffness 
C C 


for the approximations and is subject of further 


research. 


6.7.5 Implications of closure approximations based on 
second-order fiber orientation tensors 


Closure approximations (Advani and Tucker III, 1990; Han and Im, 1999; 
Cintra Jr and Tucker III, 1995; Chung and Kwon, 2002; Montgomery- 
Smith et al., 2011a) identify a fourth-order fiber orientation tensor which 
corresponds to a given second-order fiber orientation tensor. As any 
second-order orientation tensor is orthotropic or has even stronger 
material symmetry, closure approximations based on a second-order 
orientation tensor lead to orthotropic fourth-order fiber orientation 
tensors. In consequence, any virtual process chain which involves a 
closure approximation is limited to orthotropic effective stiffnesses. The 
planar orthotropic subspace is discretized in Fig. 6.4 and discussed in 
Bauer and Bóhlke (2022c). Among the discrete stiffnesses represented in 
Figures 6.9 to 6.12 only the left column, i.e., those stiffnesses labeled (y2), 
(y1), (x), (o1), (02), are orthotropic. 
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6.8 Summary and conclusions 


Computer tomography scans combined with knowledge on the variety 
of fiber orientation tensors and experience in process simulation, reveal 
that fiber architectures of sheet molding compound specimen are ap- 
proximately planar. The planarity assumption significantly reduces the 
dimensionality of the space of fiber orientation tensors from eleven to 
three in the orientation coordinate system (Bauer and Bohlke, 2022c). An 
invariant and redundancy-free parameter set of structurally differing 
fiber orientation tensors following Bauer and Bohlke (2022a) states 
the main building block for investigations on the influence of fiber 
orientation tensors on effective mechanical properties. 


The orientation average of transversely isotropic elasticity tensors follow- 
ing Advani and Tucker III (1987) formulated directly in fiber orientation 
tensors, is explicitly recast as linear invariant composition in the fiber 
orientation tensors of second- and fourth-order Kanatani third kind 
(Kanatani, 1984) . To the best of the authors knowledge, this essential as- 
sumption of the popular orientation average is not mentioned explicitly 
in the literature. It should be noted, that, e.g., in the field of biomechanics, 
material models being non-linear in fabric tensors (Cowin, 1985; Turner 
and Cowin, 1987; Biegler and Mehrabadi, 1995; Cowin and Cardoso, 
2011), are established. Such models are based on isotropic tensor func- 
tions, i.e., representation theory. A numerical orientation-averaging 
scheme restricted to the special class of planar fiber orientations based 
on a maximum entropy reconstruction of fiber orientation distribution 
functions following Bauer and Bóhlke (2022a) is proposed. The new 
scheme shows fast converges against the exact formulation of Advani 
and Tucker III (1987) for non-localized fiber orientation tensors which 
are dominant in sufficient large fiber arrangements. 


Four mean field homogenization approximations are reviewed and in- 
vestigated. The first two approximations are orientation-averaging Mori- 
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Tanaka following Benveniste (1987) and a two step Hashin-Shtrikman 
homogenization scheme (Kehrer et al., 2020) both in formulations based 
on the Advani-Tucker orientation average (Advani and Tucker III, 1987). 
The remaining two approximations are direct Advani-Tucker orientation 
averages of either a unidirectional stiffness or compliance obtained 
by the Mori-Tanaka approximation (Mori and Tanaka, 1973). Effec- 
tive stiffnesses obtained by the approximations, are visualized by two- 
dimensional polar plots of Young’s and generalized bulk modulus 
(Bohlke and Briiggemann, 2001). Plots are generated for specific points 
in the three-dimensional body of structurally differing planar fiber 
orientation tensors. The developed views establish a new methodology 
for studies and visualizations of the dependence of material models 
on planar fourth-order fiber orientation tensors. Inspection of polar 
plots of the generalized bulk modulus reveals clear structural differ- 
ences between effective stiffnesses obtained by mean field schemes 
taking the orientation average in the stiffness or compliance domain. 
Those stiffnesses averaged in the compliance domain, hardly show a 
dependence of the generalized bulk modulus on dev(NPlanar) whereas 
this dependence is pronounced for stiffnesses averaged in the stiffness 
domain. The reason for this observation and the physical interpretation 
of the pronounced dependence is subject of further research. Lack of 
fourth-order contribution to the generalized bulk modulus (Bóhlke and 
Brüggemann, 2001) restricts the possible shapes of its contour plots. 
Bauer and Bóhlke (2022a, Figure 7) visualize reconstructions of fiber 
orientation distribution functions based on leading fiber orientation 
tensors up to fourth order. The resulting plots visualize the direct 
connection between the causing orientation measure and the effected 
mechanical behavior, e.g., visualized by the effective Young's modulus 
in Figures 6.9 and 6.10. 

The main conclusions of this work are 


* Fiber orientations of sheet molding compound are approximately 
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planar, leading to a reduction of the independent components of 
fourth-order fiber orientation tensors from eleven to three in the 
orientation coordinate system. 

e An invariant and redundancy-free set of structurally differing fiber 
orientation tensors of fourth order can be used to investigate the de- 
pendence of effective material properties on fiber orientation tensors. 

* The orientation average of an elasticity tensor based on an fiber 
orientation tensor following Advani and Tucker III (1987) is linear in 
both the elasticity tensor and the fourth-order fiber orientation tensor. 

e A new numerical formulation of the Advani-Tucker orientation av- 
erage in fiber orientation tensors (Advani and Tucker III, 1987) for 
the special class of planar fiber orientations is proposed based on 
a maximum entropy reconstruction of fiber orientation distribution 
functions following Bauer and Böhlke (20222). 

* The orientation dependence of the generalized bulk modulus differs 
significantly for homogenizations which take the orientation average 
in the stiffness domain and those which take the average in the 
compliance domain. In contrast, the direction of maxima of the 
orientation-dependent Young's modulus is homogeneous for both 
groups of homogenizations and the difference on the size of the 
Young's modulus is barely influenced by the fiber orientation tensor. 

* The orientation dependence of the effective anisotropic material re- 
sponse, e.g., described by orientation-dependent Young's and general- 
ized bulk modulus, is restricted. Restrictions are caused by the limited 
averaged information given by fourth-order fiber orientation tensors 
and due to the constraints of linear elasticity. Possible directional 
dependencies of the elastic response for Advani-Tucker averaged 
two-phase materials of isotropic constituents and planar orientation 
measures are comprehensively presented. This presentation is com- 
plete and it's methodology states the main contribution of this work. 
The methodology can be used to express the orientation dependence 
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of any quantity which is a function of planar fiber orientation tensors 
up to fourth order. 

* Orthotropy of closure approximations, which are based on second- 
order fiber orientation tensors, is shown to be a major restriction on 
the effective anisotropic material response. 

The methodology developed for the low-dimensional subspace of planar 

fourth-order fiber orientation tensors may be applied to sub spaces of 

fiber orientation tensors with less symmetry. 
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Chapter 7 


Conclusions and Outlook 


7.1 Conclusions 


Combination of the detailed summaries in Sections 4.8, 5.5 and 6.8, leads 
to the following insights and conclusions. 

Significant progress has been made in the search for a complete de- 
scription of the phase space of fourth-order fiber orientation tensors. 
A complete description is obtained for selected subspaces. A natural 
parameterization of generic fourth-order fiber orientation tensors is 
given by a linear invariant composition of the isotropic reference state 
and one second- and one fourth-order deviator in an eigensystem of the 
corresponding second-order FOT. The use of the eigensystem separates 
structural properties represented by deviators and superposed rotations. 
In addition, the parameterization naturally separates second- and fourth- 
order contributions and can be refined according to material symmetries, 
thereby reducing the number of independent components significantly. 
A similar parameterization of second-order FOTs based on an isotropic 
central state and deviators increases the knowledge, e.g., on trans- 
versely isotropic states inside the commonly used orientation triangle. 
Demanding positive semi-definiteness of fourth-order FOTs, defines 
admissible parameter ranges. These ranges are given in closed form for 
transversely isotropic as well as planar orientation states. For weaker 
constraints, admissible parameter ranges can be obtained numerically 
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and are visualized for the important orthotropic case. The use of material 
symmetries artificially reduces the complexity of problems and realizes 
solvability. 


A compact and natural description of admissible fiber orientation tensors 
allows systematic investigations of processes and methods which use 
fiber orientation tensors either as input or result. The combination 
of simple and mathematically motivated parameterizations with the 
condition of positive semi-definiteness appears trivial, but holds great 
potential. For example, closure approximations which lead to fourth- 
order FOTs outside the admissible region, e.g., visualized in Fig. 4.3, 
can be easily discarded. Algebraic constraints on admissible parameter 
ranges are connected with intuitive limiting fiber arrangements by the 
following observation. Orientation states at the boundary of admissible 
parameter regions for transversely isotropic and planar states, can be 
represented by a minimal number of Dirac distributions in particu- 
lar directions, which might be called localized. For measured data, 
depending on the resolution of the measurement and the size of the 
reference volume, associated with a single measured FOT, these localized 
states are not obtained. This observation is explained by the localized 
character of such orientation states, which is in contrast to the white 
noise induced by a large number of fibers within the particular reference 
volume. Furthermore, the structure of admissible regions demonstrates 
the pronounced role of the isotropic, unidirectional or planar-isotropic 
orientation states and shows links between these special states. 

A deep understanding of the variety of FOTs is crucial for the assessment 
of closure approximations and demonstrates the ill-posed character of 
closure problems, postulating an artificial one-to-one mapping between 
second- and fourth-order FOTs. 

The transition from analytical investigations on admissible FOTs to 
mechanical properties of SMC within this work is guided by a detailed 
inspection of planar FOTs and corresponding FODF. Therefore, the 
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connection between FODF and FOT for the planar subspace is inves- 
tigated. Due to symmetry-induced ambiguities of the eigensystem, 
admissible parameter sets do not necessarily lead to distinct FOTs. This 
observation leads to a refinement of the admissible region into a set 
of admissible and distinct planar FOT of fourth order by removing 
redundancies. The observed redundancies are two-fold. On the one 
hand, full isotropy or partial isotropy in the case of transversely isotropic 
states, induce a continuous randomness of the eigensystem, which leads 
to a degeneration of the corresponding part of admissible ranges. On the 
other hand, indefiniteness or ambiguity of the eigensystem of positive 
definite second-order tensors (Bro et al., 2008), due to their inherent 
orthotropy, induces a discrete ambiguity which can be handled by a 
sign convention in a suitable parameterization. Three reconstruction 
methods of FODF by leading planar FOT lead to fundamentally different 
results. Deploying a generic truncated Fourier series in a 3D-framework 
leads to non-planar and non-positive semi-definite FODF. Truncated 
Fourier series in a 2D-framework ensure planarity, but not positive 
semi-definiteness. Maximum entropy reconstruction in a 2D-framework 
leads to reasonable reconstructions which fulfill algebraic constraints. 
The proposed explicit and minimal formulation of the maximum entropy 
method represents a valuable tool for non-linear applications which 
might require numerical integration of tensorial quantities weighted by 
FODF based on leading FOTs. The reconstruction methods associate a 
given FOT with the intuitive FODF which is directly associated with 
fiber distributions, and is more general than ensembles of discrete Dirac 
distributions. 

In order to represent reconstructed FODFs based on distinct planar 
FOTs, an innovative visualization method is developed in this work. 
This method can be used to visualize any direction-dependent scalar 
quantity, which is a function of distinct planar FOT. The method deploys 
a discretization strategy and consists of overview plots as well as polar 
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plot ensembles. The overview plots slice the admissible and distinct 
parameter space and define points on these slices. For each slice, an 
polar plot ensemble consisting of one polar plot for each point on the 
slice is given. The views on reconstructed FODFs directly associate FOTs 
with fiber alignments. 


In a subsequent step, the visualization method is applied to effective 
stiffnesses obtained by several mean field homogenizations which are 
popular in both academia and industry. The obtained views on the 
stiffnesses directly associate mechanical properties with FOTs, which 
are associated with fiber alignment through the views on reconstructed 
FODFs. In consequence, a complete study of the influence of fiber 
orientation tensors on effective stiffnesses of SMC, based on the selected 
homogenizations, is given. 

Again, consequences of deploying closure approximations to gain seem- 
ingly fourth-order information can be presented. Any reasonable closure 
approximation is orthotropic and the left most columns in Figures 6.9 to 
6.12 represent orthotropic orientation states. In consequence, only one 
third of the visualized stiffnesses can be obtained, if closure approxi- 
mations are contained within a virtual process chain. Observations on 
the algebraic structure of the selected homogenizations complete these 
investigations. 

In summary, the present work provides comprehensive contributions to 
the understanding of the variety of FOT, being a widely used directional 
measure. Interaction of this variety with fiber orientation density distri- 
butions as well as consequences for derived effective mechanical proper- 
ties are presented. These findings are seminal to the understanding and 
design of fiber reinforced plastics in a wide variety of applications. 
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7.2 Outlook 


Departing from this work’s conclusions, the following statements pro- 


vide starting points for future research questions. 


The parameterization of triclinic fourth-order FOTs implies a classifi- 
cation scheme, which can be applied to study fields of FOTs. 

Closure approximations can be systematically assessed by their abil- 
ity to remain within the admissible parameter region and maintain 
symmetry. 

The parameterizations based on deviators from the isotropic state 
naturally induces an interpolation scheme for fourth-order FOTs. 
Advanced material models, e.g., including damage or failure surfaces, 
can be defined based on the developed parameterizations. 

The planarity assumption with respect to SMC combined with the de- 
rived variety of effective stiffnesses offers great potential for efficient 
clustering techniques for mechanical properties within a virtual pro- 
cess chain. Similar to Goldberg et al. (2017), the admissible parameter 
space could be discretized. 

The visualization method presenting the dependence of effective me- 
chanical properties on planar FOTs offers the possibility to study and 
compare additional analytical as well as numerical homogenization 
schemes. 

The investigations on the dependence of material models on FOTs can 
be extended to thermo-mechanical models. 
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Appendices to Chapter 4 ' 


A.1 Material symmetries of second-order ten- 
sors 


Following Cowin (1985), any symmetric tensor of second order A is 
either isotropic, transversely isotropic or orthotropic. This means, for 
symmetric second-order tensors, three material symmetry classes exist 
and representations of A in its material coordinate system {v;} are given 
for the case of isotropy A*9, transversely isotropy A'?'*Y and orthotropy 


1 This appendix reproduces the appendix of (Bauer and Böhlke, 2022c), i.e., Bauer, J. K., 
Bóhlke, T., 2022. Variety of fiber orientation tensors. Mathematics and Mechanics of 
Solids 27 (7), 1185-1211, 10.1177/10812865211057602. Reproduced with permission. 
€2021 The Authors. Published by SAGE Publications Ltd under CC BY-NC 4.0 
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A ortho by 


1 0 0 

A's a1 0 1 0 vi ®& Vj, 
0o d 
Q1 0 

A Fansv = 0 az 0 v;® MT 
0 ag 
Q1 0 

A ortho = 0 az 0 v;® Vj. 
0 3 


(A.1) 


(A.2) 


(A.3) 


The principle direction of the transversely isotropic case is randomly 


chosen to be the first axis. The number of independent coefficients is 


one, two and three respectively and is indicated by the number of a;. 


A.2 Parameter sets of specific second-order 
orientation tensors 


Table A.1 contains arguments of parameterizations given in Equations 
(4.25), (4.29) and (4.31) leading to special orientation states. 


N, Aa) N(aı, a3) N(à €) N((, wa, 102) 
A Aa Q1 a3 a [e 01 In» W3 
isotropic 1/3|1/3 10 0 0 0 a la la 
unidirectional vı ||1 JO 2/3 |0 2/3 |—-1/3|]a |O 10 
planar isotropic vı | 0 |1/2 | —1/3/0 —1/3|1/6 |O |a Ja 
planar isotropic v3 || 1/2|1/2 ||0 -1/3 ||1/6 |-1/3|a |a |O 
Table A.1: Arguments of parameterizations of N for special orientation states with a € 7, 


i.e., a > 0. Note: The orientation states labeled by planar isotropic vı and planar isotropic 
v3, differ only by ordering convention of the eigenvalues and describe the same physical 


state. 
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A.3 Coefficient-wise extrema of moment ten- 
sors 


Given a triangulation which divides the unit sphere into m surfaces Fy 
with centered normal vectors n., for y € [1, ..., m]. The integration over 
S? in Equation (4.10) can be interpreted as a weighted summation with 
discrete weights w., derived by the FODF by 


Wy = I. (n) dn. (A.4) 


The properties of v (n) imply non-negative and normalized weights, i.e., 
0 € w, and =. wy = 1 holds. Definition of N in Equation (4.10) is the 
limiting case 


m 


; 4 
N= lim X uut (A.5) 
Y 


Defining coefficients of the Y-th moment tensor by n®? = (n2) ec Be ® 
Be, Equation (A.5) reads 


m 


Ngo un a) (A6) 
Y 


The set of weights {w, [0 < wy, $7 wy = 1) describes an m-dimensional 
region W and N is linear in these weights. Due to the linearity, the 
extreme values of N arise at the boundary of W, especially at points 
where only one weight equals one and the remaining weights are zero. 
Such points are corners of the boundary of W. As 


Y wy [(025),, — T 20 (A) 
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is valid it follows from 


m 


ru At (A.8) 
- 
2 wy (n$ ^) = (nus (A.9) 
^Y 


min) ec 


bound for the corresponding coefficient of the fourth-order orientation 


that the minimum of a coefficient of a moment tensor (n is a lower 
tensor NI. Note, that for a fixed coordinate system, there exists a 
singular FODF 7™" (n) of shape 
; 1 
ym (n) = z O(n - n") +ô (n +n*)) (A.10) 


@4 


leading to N, a = (m. ) «c Similar reasoning leads to coefficient-wise 


upper limits of N A defined by maximum coefficients of a moment 


@4 


tensor (n4) .. corresponding to a singular FODF. As a consequence, 


) 

extreme sad of tensor coefficients of N can be derived from the mo- 
ment tensor n?*. However, this coefficient-wise consideration neglects 
the interaction of the coefficients of N and the procedure defining the ori- 
entation coordinate system {v;}. Tensor coefficients of a representation 
in a coordinate system, which depends, e.g., on the eigenvectors of N, 
may not reach the extreme values reachable in a fixed coordinate system. 
Another consequence is that the norm of N is bound by the value one, 
which is reached for a singular FODF. Parameterization of n?* in two 


angles with 


n (8, p) = sin(0) cos(p)eı + sin(0) sin(y)e2 + cos(0)ea (A.11) 
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leads to a representation of the moment tensor by n?! = nee Be @Be 
with 


@4 _ A.12 
Nec (A.12) 
sick 52 stc2, s262 2/25 sic co V288e8cg — y/2s, sd cà 
6C 9559€ $905 9 ge 0^ Q?0wp 
s450 Erle, va säco V252 3c, ce Va 
4 3 3 Da Ad 
Ch C V2sgc.c 23,856, 
completely symmetric 


using the short hand notation 
sg —sin(0), co =cos(#), s,=sin(y), c.=cos(p). (A.13) 


The moment tensor has extreme coefficient-wise values 


0 0 0 —b -a -—a 1 1/4 1/4 b a a 

0 0 -a -b -a 1 14a b a 

-a -a -b 1 b 

— < ng! < = (A44) 
completely | symmetric completely | symmetric 
I m 
—min n94 = Em 
min Nge Es Neg 


with a = /2 avs and b = A, Exemplary application of the coefficient- 
wise limits on ey isotropic fourth-order orientation tensors 
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defined in Equation (4.66), leads to the set 


transv _ J wtransv inn$84 « NOD < marn?’ A.15 
{ | min nec zu dec xi (A.15) 


linear 


1 2 
= transv ee ER 
fn Te. 


1 1 
max (zx (—30a — 7), 105 (15a n) 
<ps 
T nibus w} 
840 


which is visualized in Figure 4.3. 


A.4 Parameterization of admissible N°° with 
Isotropic N 


Explicit parameterization of those fourth-order orientation tensors which 
contract to the isotropic second-order orientation tensor is given by the 


set 
{ee (N, di, de, ds) | (A.16) 

1 1 

BE an 
i5 $5 S 

bg < T 4 5 
15 ^ ?^45 ' 9(45d, — 7)' 
1 4 — 60d; — 60 dz + 675 dı da 

zm Le oe ortho. 
p 9 0 8d re, bow 
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Appendix B 


Appendices to Chapter 5 ' 


B.1 Kelvin-Mandel notation and completely 
symmetric tensors of fourth order 


This appendix directly follows Bauer and Bóhlke (2022c;a) for the current 
work to be self-contained. Kelvin-Mandel notation, explicitly intro- 
duced in Mandel (1965), originating from Thomson (1856), discussed 
in Mehrabadi and Cowin (1990); Cowin and Mehrabadi (1992); Bóhlke 
(2001) and also known as normalized Voigt notation, enables compact 
two-dimensional representations of fourth-order tensors with at least 
minor symmetry. A fourth-order tensor A = A;jj; e; Q e; Q ex ® e; is 
minor symmetric if it has both minor symmetries, i.e., Aij; = Ayıkı = 
Aij holds. Introducing base tensors in an arbitrary Cartesian basis {e; } 


1 This appendix reproduces the appendix of (Bauer and Bóhlke, 2022a), i.e., Bauer, 
J. K., Böhlke, T., 2022. Fiber orientation distributions based on planar fiber orien- 
tation tensors of fourth order. Mathematics and Mechanics of Solids, online first, 
10.1177/10812865221093958. Reproduced with permission. ©2022 The Authors. 
Published by SAGE Publications Ltd under CC BY-NC 4.0 
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by 
2 
Bi =e, ® e,, B1 = Y [e28 es + e3 eel, 
/2 
By = e2 8 e2, B; = cz [e1 Q e3 + e3 G ei], (B.1) 
2 
B3 = ea 8 es, Bs = 2 e ee +61 ©ea], 


any minor symmetric tensor A is represented by a six by six matrix of 
coefficients Ae 


A = Aij e; @ e; Q ey @ e; = Age Be 8 Be (B.2) 


with € and ¢ summing from 1 to 6. Complete index symmetry of a tensor 
N implies the structure 


N= (B.3) 
ND NO NG | VIND VING  V2Nis 

N? NG | VING VINS VƏN 

Nga) | VZN VANS VANS 


Be 9 Be. 

2ND NP on EUM 
major symmetric 2N (9 2N. (9 
an) 


As complete index symmetry implies major index symmetry, the coeffi- 
cient matrix in Equation (B.3) is symmetric. In Equation (B.3) indices of 
redundant tensor coefficients are colored. The redundancy implies that 
six coefficients in the upper left quadrant and nine coefficients in the 
upper right quadrant of the coefficients in Kelvin-Mandel representation 
define a completely symmetric tensor. This motivates a short hand 
notation ,completely symmetric", see, e.g., Equation (5.15). 
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B.2 Parameter sets in polar plots 


Table B.1 contains parameter combinations leading to selected fourth- 
order fiber orientation tensors, which are used to generate polar plots 
in figures 5.7, 5.8, 5.9 and 5.10. Numerical values of the Lagrange 
multipliers in Equation (5.72) are given with limited precision and 
absolute values smaller than 107? = 1E? set to zero. 
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Fig. Aı B Ê dı dg L fi ig gi g2 
5.7(1) 1/2 7/2 9/80 69/560 0 4.4250 0 0 —5.30EU 0 
5.7(2) 2/3 7/2 1/10 61/630 0 6.43E0 —7.69E0 0 —6.12E0 0 
5.7(3) 5/6 7/2 1/16 89/5040 0 1.66E! —2.46El 0 —9.77E0 0 
5.7(4) 1/2 | 0 o 3/280 o 8.88E-! | 0 o o o 
5.7(5) 2/3 o o —1/315 o 9.56E- 1 —6.67E-1 o 1.20E- i [n 
5.7(6) 5/6 o o —113/2520 o 1.45E0 —1.34E0 o 6.49E7 1 0 
5.7(7) 1/2 -r/2 | 9/80 —57/560 0 4.429 0 0 5.30E0 o 
5.7(8) 2/3 —r/2 1/10 —13/126 0 5.01E0 —8.63E- 1 0 5.89E0 0 
5.7(9) 5/6 | -r/2 | 1/16 —541/5040 o 8.26E0 —s8.34E-! | o 9.13E0 o 
5.7(10) 1 [n [n —4/35 [n 1.81E3 —5.51E0 1.03E! 1.81E3 —5.18E9 
5.8(x) 1/2 | 0 o 3/280 o 8.38E~! | 0 o o o 
5.8(01) 1/2 -7/2 1/16 —29/560 0 1.15E0 0 0 1.16E9 o 
5.8(02) 1/2 —7/2 9/80 —57/560 0 4.42E0 0 0 5.30E0 0 
5.8(cl) 1/2 n/a 1/16 3/280 — V2/32 V2/32 1.15E0 o 0 8.20E-1 .20E- 1 
5.8(c2) 1/2 n/a | 9/80 3/280 — 9\/2/160 92/160 4.429 o 0 3.75E9 .75E0 
5.8(m1) 1/2 0 1/16 3/280 1/16 1.15E0 0 0 0 .16E0 
5.8(m2) 1/2 | 0 9/80 3/280 9/80 4.429 o 0 o 30E? 
5.8(g1) 1/2 n/4 1/16 3/280 + V2/32 V2/32 1.15E0 0 0 —8.20E-!1 .20E- 1 
5.8(g2) 1/2 n/4 9/80 3/280 + 92/160 92/160 4.42E0 0 0 —3.75E0 .75E0 
5.8(y1) 1/2 7/2 1/16 41/560 0 1.15E0 0 0 —1.16E9 
5.8(y2) 1/2 7/2 9/80 69/560 0 4.42E0 0 0 —5.30E0 
5.9(x) 2/3 [n [n —1/315 [n 9.56E- 6.6761 [n 1.2061 
5.9(01) 2/3 —1/2 1/18 —37/630 0 1.30E0 —4.54E- 1 0 1.30E0 
5.9(02) 2/3 | -r/2 | 1/10 —13/126 o 5.01E0 3.638871 | o 5.8950 
5.9(cl) 2/3 —n/4 1/18 — 2/36 — 1/315 V2/36 1.3169 —5.86E- 1 3.36E- 1 9.53E- 1 36E-1 
5.9(c2) 2/3 —n/4 1/10 — 2/20 — 1/315 2/20 5.20E0 —1.39E0 2.51E0 4.22E0 15E0 
5.9(ml) 2/3 0 1/18 —1/315 1/18 1.32E0 —9.01E-1 4.71bE- 1 1.20E- i 19E9 
5.9(m2) 2/3 | 0 1/10 —1/315 1/10 5.69E0 —3.93E0 3.600 8.21E-2 —5.99E0 
5.481) 2/3 m/4 1/18 —1/315 + V2/36 V2/36 1.33E0 —1.21E0 3.30E7 1 —7.18E7 1 —8.40E- 1 
5.9(g2) 2/3 n/4 1/10 —1/315 + V2/20 2/20 6.21E0 —6.56E0 2.599 —4.26E0 —4.35E0 
59(y1) 2/3 7/2 1/18 11/210 0 1.3359 —1.33E0 0 —1.07E0 0 
5.9(y2) 2/3 | «/2 1/10 61/630 o 6.43E —7.69E0 o —6.12E0 o 
5.10(x) 5/6 0 [n —113/2520 [n 1.45E9 —1.34E0 [n 6.49E- 1 0 
5.10(01) 5/6 —7/2 5/144 —401/5040 0 2.00E0 —9.88E- 1 0 1.96E9 0 
5.10(02) 5/6 —n/2 1/16 —541/5040 o 8.26E0 —8.34E- 1 o 9.13E0 0 
5.10(c1) 5/6 —mn/A 5/144 —113/2520 — 5,/2/288 5V/2/288 2.0469 —1.29E0 9.33E- 1 1.5959 —9.71E- 1 
5.10(c2) 5/6 —r/4 1/16 —113/2520 — V2/32 V2/32 9.379 —4.00E0 7.92E0 6.69E9 —6.27E0 
5.10(ml) 5/6 0 5/144 —113/2520 5/144 2.07E0 —1.93E0 1.2560 6.77E- 1 —1.36E0 
5.10(m2) 5/6 0 1/16 —113/2520 1/16 1.23E1 —1.20E1 1.17E1 5.68E- 1 —9.41E0 
5.10(g1) 5/6 | 7/4 5/144 —113/2520 + 52/288 5V2/288 1.9569 —2.37E0 7.97E~ 1 —1.91E7 1 —9.20E- 1 
5.10(g2) 5/6 | 7/4 1/16 —113/2520 + V2/32 v2/32 1.55El —2.09E! 8.8750 —6.48E0 —7.31E0 
5.10(y1) 5/6 7/2 5/144 —17/1680 0 1.86E0 —2.47E0 0 —5.16E-1 0 
5.10(y2) 5/6 7/2 1/16 89/5040 0 1.66El —2.46El 0 —9.77E0 0 

Table B.1: See Appendix B.2 
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Appendices to Chapter 6 ' 


C.1 Kelvin-Mandel notation and completely 
symmetric tensors of fourth order 


This appendix is shared with Chapter 5 and contained in Appendix B.1. 


C.2 Harmonic decomposition of transversely 
isotropic elasticity tensors 


Following Spencer (1982); Lubarda and Chen (2008); Walpole (1969), 
any transversely isotropic elasticity tensor can be parameterized by five 
scalars and a direction q. For any choice of a non-unique orthonormal 
coordinate system m; with m, — q, there exists a mapping 


Q (q) =m; 6 e; (C.1) 


1 This appendix reproduces the appendix of (Bauer and Bóhlke, 2022b), i.e., Bauer, 
J. K, Bóhlke, T., 2022. On the dependence of orientation-averaging mean field 
homogenization on planar fourth-order fiber orientation tensors. Mechanics of Materials 
170, 104307, 10.1016 /j.mechmat.2022.104307. Reproduced with permission. ©2022 The 
Authors. Published by Elsevier Ltd under CC BY 4.0 
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and enables the following representation of a transversely isotropic 
stiffness 


gum (hi, ha, hg, ha, hs, q) 


Cii Cio Cia 0 0 0 
C22 Coa 0 0 0 
A C 0 0 0 
= Q(q)* m Co2— Cos 0 o | Be@Bc] (C2 
sym 2C55 0 
2055 
= ha Py + he Pa + Q (q) * [hs Js [FF] (C.3) 


+ ha Ja Er] + hs transvx 


with transversely isotropic structure tensors (Olive et al., 2018; Müller 
and Böhlke, 2016; Bauer and Böhlke, 2022c) 


1 0 0 
2 
Fitansvx _ dev (NUP) = 8 —1/2 0 ei de; (C.4) 
sym =1/2 
ptransvx _ dev (NUP) (C.5) 
8 —4 —4|0 0 0 
3 1/0 0 0 
1 3/0 0 0 
= B: e Be 


completely | symmetric 


and isotropic sixth order tensors 


Ja = [9,581 m Ó1n F ÖkıdimÖjn] e; & ej Q ey Q €] © Cm 8 en (C.6) 
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1 
J4 = 8 [din OjxO1m F ÜinÓjlÓkm Tr ÜimÓjkÓIn (C.7) 


Eg ÜimÓjiÓkn Ar ÖikdjmÖln st dikdjnOlm 
+ bi15jmOkn + O10 jndkmle; Q ej Q ex 8 ei Q em 8 En. 


The scaling factors 2/3 and 1/35 of the structure tensors in Equations 
(C.4) and (C.5) do not lead to a norm of value one. This is in contrast to 
formulations in, e.g., Fernández and Bóhlke (2019, Equation (75)), but 
leads to a direct connection to unidirectional fiber orientation tensors 
NUP = v, & v; and NUP = BY & BY and compact expressions for h; 
with i € (1, 2, 3, 4, 5} in Equation (C.3) which are related to the tensor 
components in Equation (C.2) by 


E E : OEEO EEO RET (C.8) 
pcdes = [DC — 4C yo + 7022 — Cog + 12055] CY 
jys : Catit O= T6 dog] (C.10) 
nes L [2C11 — 4C ya — 5022 + 7C23 + 6Cs5] (C.11) 
hg = = [C11 — 2C12 + C22 — 4055] (C.12) 


with the bulk modulus X and the shear modulus G. Equations (C.1) to 
(C.3) motivate a short hand notation for a transversely isotropic elasticity 
tensor using | - | with 


(ctransv (hi, ha, hs, ha, hs, q) = Jin, ha, ha, ha, hs, al * (C.13) 


Following Rychlewski (2000), Equation (C.3) represents the classical lin- 
ear invariant decomposition of a transversely isotropic elasticity tensor 
into irreducible parts. The classical linear invariant decomposition is one 
out of an infinite number of possible linear invariant decompositions 
(Forte and Vianello, 1996; Rychlewski, 2000). Any non-classical linear 
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invariant decomposition leads to a representation of C*?^** which differs 
from Equation (C.3) by the set of isotropic tensors J3, J4 and the values 
of the corresponding coefficients ha and h4. 


C.3 Reformulation of the Advani-Tucker orien- 
tation average 


Listing C.1 contains symbolic tensor manipulation code to validate 
the representation of the Advani-Tucker orientation average in the 
formulation given in Equation (6.16) and which is based on the classical 
linear invariant decomposition (Rychlewski, 2000). 


Listing C.1: Validation of the orientation average following Advani and Tucker III (1987) 
given in Equation (6.16) 


import sympy as sp 
import symbolic as sb 


D2 = sb.fabric tensor.dev2 by 1a0 lal() 
D4 = sb.fabric tensor.dev4 triclinic by d( 
N4 = sb.special.combine to N4(D2-D2, D4=sb.tensorr (D4) 


# Equation (C.3) Transvsely isotropic stiffness 
stiffness = sb.material_symmetry. 
transversely_isotropic_x_minimal_harmonic_normalized ( 


normalization-"structure tensor multiplicity" 


print ("Orientation tensor fourth-order Kanatani first kind =\n", N4) 
print ("Transversely isotropic stiffness =\n", stiffness) 


averager = sb.orientation_average.AdvaniTucker (N4=N4) 
average = averager.average (stiffness) 


decomposer Rych = sb.elasticity.Rychlewski2000 () 

decomposition = decomposer_Rych.decompose_classical_harmonic ( 
stiffness=sb.tensorr (average) 

)["parts"] 
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def test_equality(A, B): 
return sp.simplify(A) == sp.simplify (B) 


assert test equality (decomposition["h1"], sb.abc.hl) 

assert test equality (decomposition["h2"], sb.abc.h2) 

assert test equality(decomposition["H1"], D2 " sb.abc.h3) 

assert test equality (decomposition["H2"], D2 " sb.abc.h4) 

assert test equality (sb.mandel(decomposition["H"]), DA " sb.abc.h5) 


C.4 Advani-Tucker Orientation Average of Mi- 
nor Symmetric Tensors 


A generic transversely isotropic fourth-order tensor, which has both 
minor symmetries but lacks major symmetry, i.e., Hijri = Hijik = Hjikl 
but Hijki £ Ag; with Hransv — Hijkxiei & e; 6 ek 8 E, depends on 
a direction q and six scalars, see Schróder and Gross (2004); Brannon 
(2018). For any choice of a non-unique orthonormal coordinate system 
m; with m, = q, there exists a mapping 


Q (q) = m;Ge;, (C.14) 


which enables the following representation of a transversely isotropic 
minor symmetric fourth-order tensor lacking the major symmetry 


H" (Hii, H22, Hi2, Haı, Has, H55, q) = (C.15) 

Hıı Hıa Hıa 0 0 0 
Hz, H22 Has 0 0 0 
A Az, Hos | Hae 0 0 0 

* B:®B 

ala) Ü. 0 0 Ha-Ha 0o 3e $2 
0 0 0 0 2 H55 0 

0 0 0 0 0 2 Hss 
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In other words, if the transversely isotropic axis is aligned with the 
eı-axis of the coordinate system underlying the Kelvin-Mandel basis 
B; & B, in Equation (C.15), the matrix representation of Kelvin-Mandel 
components has the simplified structure shown in Equation (C.15). 
Compared to the representation of C*?"* in Equation (C.2), lack of the 
major symmetry of H" results in one additional independent tensor 
component in Equation (C.15). In order to extend the reformulation 
of the Advani-Tucker orientation average in Equation (6.16) to tensors 
which lack the major symmetry, harmonic decomposition (Forte and 
Vianello, 1996; Rychlewski, 2000) of such tensors is investigated follow- 
ing Fernández and Bóhlke (2019, Section 2.4) and Lobos et al. (2017, 
Appendix A). Lack of major symmetry of a minor symmetric fourth- 
order tensor, potentially adds three additional sub spaces of orders one, 
two and three, to the sub spaces of its harmonic decomposition. If the 
decomposed quantity is major symmetric, the contributions of these 
three additional sub spaces are zero (Lobos et al., 2017, Section 2.4). It is 
easy to show, that no non-vanishing first-order tensor a fulfills 


Qxa=a V Qest™ (C.16) 


with the transversely isotropic symmetry group S"?"°V. Following Bran- 
non (2018), the same holds for harmonic tensors of third order. In 
consequence, the additional degree of freedom of H'?"** in Equation 
(C.15) compared to C*?"** in Equation (C.2) corresponds to the additional 
subspace of second order in the harmonic decomposition of minor 
symmetric fourth-order tensors which lack major symmetry. As the 
symmetry group of a tensor is the intersection of the symmetry groups 
of its harmonic parts (Forte and Vianello, 1996), transversely isotropy of 
"sv demands that this subspace is also transversely isotropic, making 
it one-dimensional and being a multiple of a structure tensor F'?rsv* 
defined in Equation (C.4). Consequently, H'*?"*" can be represented by 
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extending Equation (C.3) to 


]transv (hi, ho, ha, ha, hs, he, q) (C.17) 
= hy Py + ha P» + Q (q) de [hs Js [pres din lia Ja [pras] 


+ hs rtransvx 4 he Is | 


with F'nsvx, Firansvx I^ and J4 defined in Equations (C.4) to (C.7) and 
with the isotropic sixth-order tensor 


Je = [—5ij5kn dim + Öindzmöri] €i Q e; Dek Ge; Ges Ben — (C.18) 
following Fernández and Bóhlke (2019), acting as 
Js. [H] = -IG H + Heol. (C.19) 


The connection between h; for i € [1, 2, 3, 4, 5] and components H;; in 
Equation (C.15) follows the pattern in Equations (C.8) to (C.12) combined 
with 


1 
he = 2 (Hi m H»1) " (C.20) 


Based on the representation of H*?"** in Equation (C.17), the Advani- 
Tucker orientation average (Advani and Tucker III, 1987) given in equa- 
tion (6.16) can be extended to minor symmetric fourth-order tensors 
which lack the major symmetry, leading to 


(HE®S (hy, ha, ha, ha, hs, he, q)) apn (N) = hi Pi + ho P2 (C.21) 
+ ha Ja [dev(N)] + ha Ja [dev(N)] + ha dev(N) + he Je [dev (N)] 


with dev(N), dev(N) defined in Equations (6.17) and (6.18). An alterna- 
tive representation closely following the original formulation of Advani 
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and Tucker III (1987) is given by 


(EEES mw) iia = by NS (C.22) 
+ bə Nij Ont + be Ni Oi; 

+ bg (Nij ji + Nu jy + Nji dir + Nyx Sir) 

+ b4 di; Örı + bs (dir O51 + Sir Ok.) 


with component representations in an arbitrary coordinate system 


N= NA e; Be; Pez G e (C.23) 


The coefficients b; for i € [1, 2, 3, 4, 5, 6] are defined by 


b, = Hiii + H2222 — H1122 — H2211 — 4 H1212 (C.25) 

bə = H1122 — H2233 (C.26) 
1 

b3 = Ay212 + P] (H2233 — H2222) (C.27) 

bs = H2222 — H2233 (C.29) 

bs = H2211 — H2233 (C.30) 


based on tensor components of H" = H,jki ei Q ej Q ex Q & ina 
coordinate system with e; aligned along the transversely isotropic axis 
of H'*"s’, i.e., the tensor components H;;,; are directly related to the 
Kelvin-Mandel components H;; in Equation (C.15). If the quantity to 
be averaged is major symmetric, the coefficients in Equations (C.25) 
to (C.30) coincide with those for the original formulation Advani and 
Tucker III (1987), e.g., given in Brylka (2017, Equation (2.89)) or Kehrer 
et al. (2020, Equation (27)). 
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C.5 Connection to Notation in Kehrer et al. 
(2020) 


The two-step scheme in Equations (6.32) and (6.33) simplifies to the 
one proposed in Kehrer et al. (2020), if kı = 0 or kı = land k = kə. 
Equations (C.31) and (C.32) connect the notation of this work on the left 
hand side and the notation of Kehrer et al. (2020) on the right hand by 


CHW? (kı = 0, ko =k, ...) = CH (k), (C.31) 
CHSW2 (kı = 1, ka = k, ...) = CB (k). (C.32) 


C.6 Component representations of tensor in- 
versions 


Explicit representations of components after matrix inversion in Equa- 
tion (6.43) are given by 


h? = 3 (105h5 + 140h4 + 36hs) /b (C.33) 

hb = 105(33075h4 h2 + 22050hıhaha + 5670hıhahs (C.34) 
— 17640hı h3 — 3780hıhahs — 864hı h? — 52920h3h2 
— 141120hahzhy — 94080h5h2 + 17640h3h4 + A536h2 hs 
+ AT040ha 4 + 12096hah4hs + 313608; + 806A hs) / (512) 

hË = 45(102900hhah4 + 98000h1h2 — 2520hih4hs — 3600h1h2 (C.35) 
— 77175h3h3 — 102900h2h4 — 441003 h2 — 66150h2h3h4 

+ 13230h2h3hs — 9800h3h2 + 17640hahahs — 117600h3h4 

+ 22680h2hs — 245000h3h2 + 34020h3h4hs + 1008h3h2 

— 117600h} + 5040h2hs + 1344h4h2)/(bib2) 


E 
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hb = —225(15435hihah4 + 14700h4 h2 — 378hıhahs — 540hıh? (C.36) 
— 6615h2h3 — 17640h2h3h4 — 11760h3h2 — 17640h2h4 
+ 3402h2 hs — AT040hah2 + 9072h3h4hs 
— 313601} + 6048h2 hs) / (bi 63) 

hÈ = —1575(2205h;hahs — 5880hih2 — 2730hıhahs + 216hıh? (C.37) 
— 6615h2h3 — 17640hahaha — 11760h3h2 + 13230h2h4 
— 1008h3hs + 35280h3h2 — 2688h3h4hs + 23520h7 
— 1792h2hs) /(b1b2) , (C.38) 


with 
by = 315hiha + 420hıha + 108h4 hs — 630h3 — 1680h3h4 (C.39) 


— 1120h3 
bə = (105h2 — 140h4 + 6h) (105h2 + 70h4 — 24ha) , (C.40) 


after orientation-averaging in Equation (6.44) by 


1 
Eu = ^35 [35 Budi = 35B11Aı + ABıı m 70Bıadı ad 8D1» (C.41) 
+ 35Bood, + 35B33A1 — 31 Bəə — 140Bs5d, — 16Bs5 
1 
Eis = 2 [35 Budi +4Bıı — 70Bıadı + 27By + 35 Baadı (C.42) 


+ 4B»5 — 140 B55dı — 16Bss 


E13 = B13À1 — Ba3\ı + B23 (C.43) 

Eis = V2dg (Bi — 2Bi2 + B22 — 4Bss) (C.44) 
1 

E2 = ^35 [35B11dı +35B11Aı — 31B11 — 70Bı2dı — 8Bı2  (C.45) 


+35 Baadı — 35B33À1 + ABaa — 140Bgsd, — 16Bss 
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E53 = — B121 + Big + B23A1 

E33 = B22 

E44 = B22å1 — Ba3Aı — 2.B55A1 + 2Bss 

Ess = -BaaAı + B22 + B231 — B23 + 2B55A1 


Ess = 


2 
35 
+ AB55 — 140Bs5d1 + 19Bs5 


[35 Budi + ABy — Brod; — 8B12 + 35 Baadı 


(C.46) 
(C.47) 
(C.48) 
(C.49) 


(C.50) 


with dı = d, + R(À1) — & (see Bauer and Böhlke (2022a)) and after 


35 


matrix inversion in Equation (6.45) by 


2 Es E33 — E22 E33 E66 + E33 E66 


b 


_ Er2E33E66 — E13 E23 E66 + Eig Esa 


b 
E12E23 Egg + Ei5 E] — E13 E22 Ego + Efe E23 


b 
Eis (E12 Es3 — E13 Eo3 + E22E33 — E34) 
b 
E11 E33 Eoo — Ets Eoo — Eig E33 


b 


_ Ej Eos Eg — E12E13 Egg — E13 Ej — Eig E23 


b 
Eis (Eti Ess + E12E33 — E24 — E13 E23) 


b 


E Ey, EQ — Ey, Ej Egg + E]; Egg + 2E19 Ej + E25 E22 


b 


= Eıs (E11 E23 — Ej Eia + Ey2E23 — Eig E22) 


b 
1 


Eu’ 


(C.51) 
(C.52) 
(C.53) 
(C.54) 
(C.55) 
(C.56) 
(C.57) 
(C.58) 


(C.59) 


(C.60) 
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= —_ C.61 
G55 Ess ( ) 
Ges = — [Bus B22 Ens — Eu E34 — Ea E33 + 2E12 E13 Eo (C.62) 

z Ej En] n 


b = Eu E? E33 — E11 Es; E33 Egg + E11 E24 E66 + Ez E33 Egg (C.63) 
— 2F 19 E13 E23 E66 + 2E 12 E25 E33 — E24 E28 + E24 E25 Eco 
— 2E; 4 E2, E23 + Ed, E23 E33 — E, E2,. 


C.7 Parameter sets in polar plots 


Table C.1 contains parameter combinations leading to selected fourth- 
order fiber orientation tensors, which are used to generate polar plots in 
Figures 6.7, 6.8, 6.9, 6.10, 6.11 and 6.12. The corresponding parameteriza- 
tion NPlanar(A, f. B ) is defined in Equation (6.11). Numerical values of 
the Lagrange multipliers (L, fi, fo, 91, 92) specifying j^ (gi, ye ener) 
following Bauer and Bóhlke (2022a, Equation (71)) are given with abso- 
lute values smaller than 1075 = 1E? set to zero. 
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C.8 Effective complementary elastic energy 
density 


Effective complementary elastic energy density due to uniaxial tension 
with unit stress c = 1 MPa defined in Equation (6.66) for mean field 
homogenizations listed in Figure 6.5b are given in Figures C.1a to C.3. 
The bounds of the contour plots in Figures C.1a to C.3 are homogeneous 
and given by 0 MJ/m? and 0.1 MJ/m?. 
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nur 


(b) Polar plots for selected N specified in Figure 6.5a and Table C.1. The order along the 
path in Figure 6.5a is given from left to right. 


Figure C.1: Effective complementary elastic energy w* (C,o = 1MPa, n) due to uniaxial 
tension in varying direction defined in Equation (6.66) for mean field approximations 
specified in Figure 6.5b. 
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Sanpa” “anm” 


Figure C.2: Effective complementary elastic energy w* (C, o = 1MPa, n) due to uniaxial 
tension in varying direction defined in Equation (6.66) for mean field approximations 
specified in Figure 6.5b. Polar plots for selected N specified in Figure 6.6a and Table C.1 
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(02) 4 (c2) 


RII i |I 


Figure C.3: Effective complementary elastic energy w* (C, 0 = 1MPa, n) due to uniaxial 
tension in varying direction defined in Equation (6.66) for mean field approximations 
specified in Figure 6.5b. Polar plots for selected N specified in Figure 6.6b and Table C.1 
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